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ABSTRACT 

The postbuckling failure of square composite plates with central holes is analyzed nu- 
merically and experimentally. The particular plates studied have stacking sequences of 
[ + 45/0/90] zs , [ + 45/0 2 ] 2S , [ ± 45/0„]s, and [ ± 45] 4S A simple plate geometry, one with a hole 
diameter to plate width ratio of 0.3 , is considered. Failure load, failure mode, and failure lo- 
cations are predicted numerically by using the finite-element method. The predictions are 
compared with experimental results. In the experiments in order to be accomodated by the 
test fixture it is necessary for the plates to be slightly larger than the analysis region, ex- 
tending somewhat beyond the supports. The region outside the supports is included in the 
numerical study. It is shown that not considering these regions can lead to erroneous nu- 
merical predictions. In numerical failure analysis the interlaminar shear stresses, as well as 
the inplane stresses, are taken into account. By comparing the interlaminar shear stress 
calculations from the finite-element method with analytical results for simple cases, a solid 
foundation for interlaminar shear stress calculation is established. As a failure criterion, the 
maximum stress criterion is used. A special test fixture was designed for loading the plates. 
In the experiments, strain gauges, linear variable displacement transducers (LVDT s), and the 
shadow moire and acoustic emission tests are used to monitor plate response. An issue ad- 
dressed in the study is the possible mode shape change of the plate during loading To ac- 
count for the fact that plates can experience mode shape change, two deformation 
configurations are considered. One configuration assumes the plate responds with one half- 
wave in the loading direction and the other configuration assumes the plate responds with two 


II 



half-waves in the loading direction. Failure predictions are made for each configuration. No 
attempt is made to predict the configuration change. However, it is predicted that the first 
three laminates fail due to excessive stresses in the fiber direction and, more importantly, that 
the load level is independent of whether the laminate is deformed in the one half-wave con- 
figuration or the two half-wave configuration. The location of failure does depend on the de- 
formed configuration. It is predicted that the fourth laminate fails due to excessive inplane 
shear stress. Interlaminar shear failure is not predicted for any laminate For the Hrst two 
laminates the experimental observations correlated well with the predictions. In the exper- 
iments the third laminate failed along the side support due to interlaminar shear. The 12 
clustered 0° layers resulted in very low interlaminar shear strength S 23 The fourth laminate 
failed due to inplane shear in the location predicted However, material softening resulted in 
an actual failure load quite different from the predicted value. 
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1.0 Introduction 


1.1 Overview 


Since plates with circular holes and other openings are extensively used as 
structural members in aircraft design, the buckling and postbuckling behavior of such 
plates has received much attention. Much of the work done in this area has occurred 
since the mid-1940's. Owing to the abundant accumulation of information and 
understanding of the mechanical behavior of isotropic plates, including plates with 
holes, isotropic plates are often designed to operate in the postbuckling range. As 
an example, fuselage panels are frequently designed to sustain loads of three or four 
times the buckling load [1]. This results in a great benefit of weight savings. 

With the development of advanced composite materials, the potential application 
of composite materials to aircraft design has increased rapidly. Consequently, the 
mechanical behavior of composite plates in the buckling and postbuckling range of 
loading has become a topic of great interest. As the postbuckling response of com- 


Introduction 


1 


posite plates is not yet fully understood, the initial buckling load is usually considered 
as the ultimate load in composite material applications. For increased efficiency, 
however, operation in the postbuckling range is important and, as a result, many in- 
vestigators are conducting research regarding the postbuckling behavior of compos- 
ite plates for various geometries, loadings, and support conditions. Recent 
developments [2] show that even composite plates with large holes have load bear- 
ing capability beyond the buckling load. The key issue is determining failure load 
levels in the postbuckling range. Failure in these situations involves a coupling of 
structural response and material failure. It is a complex topic and because of it, the 
failure of postbuckled composite plates with holes is not yet fully understood. The 
present study is devoted to increasing the understanding of the failure of composite 
plates with holes loaded in the postbuckling range of response. To follow is a review 
of some of the literature felt to be relevant to the topic studied. Cited are buckling 
and postbuckling studies of plates with holes, and postbuckling failure studies of 
composite plates. After the review of the literature, the specific problem addressed 
in this study are described. 


1.2 Literature Review 


Buckling analysis of isotropic plates began with the work of Bryan [3] a century 
ago. However, the history of buckling analysis for composite materials does not go 
back that far. It began with the analysis for plywood by March [41 in 1942. With the 
development of advanced composite materials and increasing need for light-weight 
structural members for flight vehicles, scientists and engineers began to study the 
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mechanical behavior of composite materials. The history and nature of buckling an- 
alyses of composite plates to 1987 are well explained in Leissa s review paper [5]. 


1,2.1 Review of Buckling and Postbuckling of Plates with Holes 

Inevitably some plates used as structural members will contain holes. These 
holes can be access holes, holes for hardware to pass through, or in the case of 
fuselages, windows and doors. In some cases holes are used to reduce the weight 
of the structure. When a plate contains a hole, it is well known that the tensile 
strength is reduced due to the stress concentration around hole. However, the re- 
duction of tensile strength may have little bearing on the buckling behavior of a plate 
with a hole. The two problems are quite different. For the buckling problem the 
prebuckling stresses must first be found, and then the stability of the plate in the 
presence of these prebuckling stresses be studied. The former problem is an inplane 
problem, like the tensile problem, and involves inplane stiffnesses. The stability 
problem, however, involves out-of-plane effects, in particular, bending stiffnesses. 
Thus the mechanics of the buckling problem are much more complicated than the 
mechanics of the tensile problem. The study of the buckling of an isotropic flat plate 
with a central circular hole began with the work of Levy, Wooley, and Kroll [6], Fol- 
lowing Levy, et al. many investigators studied the stability of plates with central holes 
using various methods to find the prebuckling and buckling solutions. The history of 
solution techniques and the findings for the stability of plates with holes to 1983 is 
explained in detail in the literature review of Nemeth's thesis[7]. A brief summary 
of Nemeth's literature review is presented in Table 1-1 [6, 8-20]. In the present review 
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of buckling and postbuckling of plates with holes only the papers and reports after 
1983 will be explained in detail. Herein, the length of the plate will be denoted as L, 
the width W, and the hole diameter D. An important parameter when discussing 
plates with holes is the hole diameter to plate width ratio D/W. Another important 
parameter is the plate aspect ratio L/W. These notations will be used in this study. 

Nemeth [7] investigated the buckling behavior of orthotropic composite plates 
with central holes by using Kantorovich's method in both the prebuckling and 
buckling analyses. The loading conditions considered were both uniform stress and 
uniform displacement on opposite edges. The adjacent unloaded edges were simply 
supported and the loaded edges were simply supported or clamped. In addition to 
isotropic and orthotropic plates, quasi-isotropic and symmetric angle-ply plates with 
very weak anisotropy were analyzed by neglecting D 1( andD 2 e . From the buckling 
equation anisotropic parameters were derived. These parameters were used to de- 
termine the degree of anisotropy of the plate. The details of this argument were 
published in ref. 21. The results from this approach were compared with finite- 
element results. For isotropic plates good correlation with the previous results was 
obtained. Laminates with all 90° layers appeared to be the least sensitive to changes 
in hole size, support conditions, and loading conditions. As the number of 0° layers 
increased, the sensitivity to changes in these variables increased. Quasi-isotropic 
plates exhibited buckling behavior similiar to that of isotropic plates. For symmetric 
angle-ply plates fiber orientation and type of the loading substantially affected the 
buckling behavior. The [( ± 60) 8 ] s and [(±45),] s laminates showed high buckling re- 
sistance. No experimental results were reported. 
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Table 1-1 Buckling and postbuckllng analysis for plates with holes [7]. 

(continued) 


year 

Investigator 

Num./Exp. 

Remarks or Findings 

1947 

Levy, 

Wooley, 

Kroll 

[6] 

Num. 

- First study on buckling for a plate with hole 

- Kirsh solution for prebuckling and 
energy method for buckling 

- Monotonic decrease with hole size for simple 
support and stress loading 

1949 

Kroll 

[8] 

Num. 

- Similar approach to Levy, et al. 

- Shear loading is considered 

- Buckling load is more sensitive to hole 
size than uniaxial loading 

1951 

Kumai 

[9] 

Num./Exp. 

- Same approach as above 

- w satisfies stress free condition on hole 
boundary 

- For uniform stress loading simple support and 
clamped considered 

- Simple support : monotonic decrease 

- Clamped : reaches minimum at about D/W = 0.2 
then increases 

1964 

Schlack 

[10] 

Num. /Exp. 

- First study for edge displacement loading 

- Rayleigh-Ritz method for prebuckling 
and buckling 

- Biaxial loading test apparatus designed 

- Concluded singular term is not necessary 

- Studied up to D/W = 0.3 

1967 

Yoshiki, 
Fujita, 
Kawamura, 
Arai[1 1] 

Num./Exp. 

- In addition to size effects of the hole 
the effect of hole position was examined 

- Substantial loss of buckling load is possible 
for eccentric hole 

1968 

Schlack 

[12] 

Num. /Exp. 

- Studied up to D/W = 0.7 

- Rayleigh-Ritz solution agrees well up 
to D/W = 0.35 

1968 

Kawai, 

Ohtsubo 

[13] 

Num. 

- FEM for prebuckling and Rayleigh-Ritz for 
buckling 

- Good agreement with previous experiments 

- The role of prebuckling stress distribution 
explored 

1971 

Yu & Davis 
[14] 


- The prebuckling stress distribution for 
plate with small hole and large hole is 
sustantially different 

1972 

Martin 

[15] 

Num./Exp. 

- First study on composite plate with hole 

- Rayleigh-Ritz method used 
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Table 1-1 Buckling and postbuckling analysis for plates with holes [7]. 

(concluded) 


year 

Investigator 

Num./Exp. 

Remarks or Findings 

1972 

Vann & Vos 
[16] 

Num, 

- FEM used for isotropic plate 

- The interaction between prebuckling stress 
distribution and reduction of bending stiffness 
due to the presence of hole addressed 

- Loss of bending stiffness dominates for 
small hole and prebuckling stress 
dominates for large hole 

- When the aspect ratio is two or greater, the 
difference between stress and displacement 
loading diminishes 

1973 

Yu & Davis 
[17] 

Exp. 

- Postbuckling, shear buckling, web crippling 

1975 

Ritchie, 

Rhodes 

[18] 

Num./Exp. 

- Buckling and postbuckling, isotropic plate 
FEM for prebuckling and Rayleigh-Ritz 
for buckling 

- Up to D/W = 0.4 postbuckling prediction 
agreed with experiment 

1978 

Knauss, 

Starnes, 

Henneke 

[19] 

Exp. 

- Buckling , postbuckling, failure study 

- Aspect ratio 2, D/W up to 0.3 

- At D/W = 0.15-0.2 mode shape changes from 
full sine wave to half sine wave 

1980 

Proebrazh- 

enskii[20] 


- Survey paper for mostly Eastern Europe and 
some of America, Western Europe, and Asia 
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In 1983 Azizian and Roberts [22] investigated the buckling and elasto-plastic col- 
lapse of an isotropic square plate with a hole. The finite-element method was used 
for this study. The shape of the hole was either square or circular. The support con- 
dition was either simple support or clamped on all four sides. The loading was a 
uniform displacement on one edge. There was no difference in the buckling load 
between circular and square holes up to D/W = 0.3 for the simply supported case. For 
D/W > 0.3 the plate with a square hole exhibited a slightly higher buckling load. For 
the clamped case the plate with a square hole showed a higher buckling load re- 
gardless of the hole size. In both support cases the decrease in buckling load due 
to the holes was very small, and for a large hole the buckling load was even higher 
than for the plate without a hole. Even though the buckling load of a plate with large 
hole increased, the elasto-plastic collapse load decreased monotonically as the hole 
size increased. In addition, there was no difference in collapse load between the 
circular and square hole. When the plate was thin, the reduction in collapse load due 
to the hole was small. When the plate was thick, the reduction was pronounced. 

In 1985 Vandenbrink and Kamat [23] studied the buckling and postbuckling be- 
havior of isotropic and laminated composite square plates with central circular holes. 
They used three-dimensional finite elements in their analysis. The laminates ana- 
lyzed were symmetric angle-ply laminates, ([±0]s) . where the fiber angles were 
6 = 0°, 30°, 45°, and 60°. They studied initial imperfection sensitivity in the 
postbuckling behavior of plates with a central hole. They concluded that these plates 
are insensitive to initial imperfection and the postbuckling strength of a plate with a 
hole is reduced compared to a corresponding plate without a hole. 
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Marshall, Little, and El-Tayeby [24] studied the stability of composite panels with 
holes by using an expression which had a one-term double cosine and a one-term 
exponential for the out-of-plane displacement. They followed a formal energy method 
solution procedure. They compared their results with their experimental results and 
previously reported results. They obtained good agreement. 

In 1986 Nemeth, Stein, and Johnson [25] performed experiments for orthotropic 
[Oi 0 ]s. [90io] s , and [(0/90) 5 ] s plates and compared the buckling loads with those from 
Nemeth [7], In the experimental apparatus the plates were clamped on the loaded 
edges and simply supported on the unloaded edges. In all cases the difference be- 
tween experiment and analysis was less than 13 %. 

Sabir and Chow [26] obtained the solutions for the elastic buckling load of square 
plates with eccentrically located holes. Uniaxial, biaxial, and pure shear loading 
cases for simply supported and clamped plates were considered. Generally, the 
buckling loads decreased with eccentricity for holes larger than D/W = 0.2. The re- 
duction in the buckling load was more appreciable for clamped plates. The excep- 
tional case occurred when the circular hole was within the compression zone of a 
plate loaded by a uniform shear. The buckling load in this case increased appreciably 
with eccentricity for both simply supported and clamped plates. 

Yettram and Brown [27] studied the effect of biaxial loading on the buckling be- 
havior of isotropic plates with central holes. Uniform edge stress was applied on all 
four edges. The buckling behavior was studied by changing the ratio of the load in the 
x direction to the load in the y direction, and the size of the hole. Generally, for all 
support conditions the buckling load was sensitive to the hole size when the loads in 
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the x direction were opposite in sign to the loads in the y direction. When the loads 
were both compressive, the dependence of buckling load on the hole size was not 
pronounced. 

Marshall, Little, and El-Tayeby [28] investigated the membrane stress distribution 
in the postbuckling range for composite plates with holes. They showed that when 
the hole size was small, D/W = 0.2 , the stress concentration effect near the hole 
boundary diminished and the areas near the supported edges were subjected to 
higher stresses. When the hole size was large, D/W = 0.4 , the stress concentration 
effect still dominated. 

Larsson [29] analyzed the buckling and postbuckling behavior of square plates 
having central circular holes by the finite-element method. Both uniaxial and biaxial 
compression were considered for various types of boundary conditions. They also 
studied an annular plate with a hole. 

For a single layer the dependence of buckling load on fiber orientation was in- 
vestigated by Turvey and Sadeghpour [30]. They found that the buckling loads of 
square carbon-fiber reinforced and glass-fiber reinforced laminates without a hole 
decreased monotonically as the fiber angle changed from 0°to 90° with respect to the 
loading direction. The loading was a uniform stress loading and the support condi- 
tions were both simple and clamped support. The inplane movement on the unloaded 
edge was not constrained. When the hole size was D/W = 0.3, the buckling load of the 
simply supported plate was not sensitive to the fiber angle change. However, the 
buckling load was sensitive to the fiber angle change for the clamped support. For 
clamped supports the buckling load also decreased monotonically with fiber angle. 
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In 1988 Nemeth [31] reported experiments focused on the buckling behavior of 
symmetric angle-ply plates with holes. The results were compared with his numerical 
results from ref. 7. He conducted a parametric study by changing the hole size, 
loading and support conditions, and the fiber angle of the laminate. For [( ± 15) m ] s 
laminates, m ^ 6 , the dependence of buckling load on the aspect ratio, L/W, was 
studied for simply supported and clamped conditions. The results indicated that the 
buckling load was independent of the support condition when the aspect ratio was 
greater than 2.5 

In 1989 Lin and Kuo [32] performed several parametric studies for the buckling 
behavior of laminated plates with a central hole with a shear deformable finite- 
element formulation. They studied the dependence of buckling load on hole size for 
uniaxial and biaxial loading, for simple and clamped supports, and for different com- 
binations of anti-symmetric cross-ply laminates. They also studied the dependence 
of buckling load on the thickness of the plate for different hole sizes. The dependence 
on the longitudinal-to-transverse stiffness ratio, and the dependence on the biaxial 
loading ratio, were also studied. For anti-symmetric angle-ply plates the dependence 
on the fiber orientation was studied for different hole sizes. 

Nemeth [2] investigated experimentally the postbuckling behavior of square 
composite plates with central holes. For [Oi 0 ]s, [9Oi 0 ]s. [(0/90) 5 ] s , [( ± 30) e ]s, 
[( ± 45) e ] s , and [( ± 60) e ] s laminates with hole sizes from D/W = 0 to 0.6 postbuckling 
behavior was investigated. Results for aluminium plates were reported also. It was 
observed that some plates failed at low load levels. As an example a [9Oi 0 ]s plate 
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failed at low loads when the hole sizes were large. However, no plate failed before 
the buckling load was reached. This means that even plates with large holes exhibit 
some postbuckling strength. Depending on the hole size the postbuckling stiffness 
changed. Most of the plates with holes showed smaller postbuckling stiffness than 
the plates with no hole. However, some of the highly orthotropic plates with holes 
showed higher postbuckling stiffness than the corresponding plate without a hole. 
Failure predictions were not part of this study. 

Hyer and Lee [33] studied the buckling load using the curvilinear fiber format in 
a composite plate with a circular hole. By using the finite-element method and opti- 
mization technique they tried to find the curvilinear fiber path which produced the 
maximum buckling load. They showed that the buckling load of a composite plate 
can be increased by using the curvilinear fiber format instead of the conventional 
straight fiber format. They also showed that by combining curvilinear fiber format 
layers with the straight fiber layers both the buckling load and tensile failure load can 
be increased. 


1.2.2 Review of Postbuckling Failure 

There are few references in the literature dealing with the postbuckling failure 
of composite plates. Starnes and Rouse [34] presented extensive experimental re- 
sults for composite plates. All of the plates were rectangular and their aspect ratios, 
L/W, ranged from 2.1 to 6.7. Most of the plates were solid, without holes. Some 
plates, however, had small circular holes, and some had been subjected to impact 


Introduction 


11 




damage. For a few plates experimental results were compared with numerical re- 
sults from a general-purpose FEM code. They found that most of the specimens 
failed along a nodal line of the buckling mode in an interlaminar shear failure mode. 
Failure predictions were not included. 

Gurdal, Haftka, and Starnes [35] studied the effect of slots on the buckling, 
postbuckling, and postbuckling failure of laminated plates. Centrally located slots did 
not affect the failure loads or the failure modes of the plates. Slots located close to 
an unloaded edge of a plate reduced the postbuckling strength of the plate and af- 
fected the failure mode. For one slightly offset slot location, interaction of out-of-plane 
deformations with the slot tip was observed to affect the failure mode significantly 
without affecting the failure load. 

Buskell, Davies, and Stevens [36] showed other evidence that in postbuckling the 
rectangular plates which had aspect ratios, L/W, of 3 - 6.9 failed at the nodal lines and 
interlaminar shear stresses were the governing parameters. They used accoustic 
emission tests to monitor the failure development, and with ultrasonic techniques 
they observed the development of delamination growth from the unloaded edges. 

Noor, Starnes, and Waters [37] published numerical failure analysis results which 
were compared with limited experimental data for composite plates of aspect ratio 
2.0. All of the plates were solid plates without holes. As a numerical tool a mixed 
finite-element formulation was used. The formulation included shear deformation and 
each node had 13 degrees of freedom. Up to the buckling load the correlation be- 
tween numerical and experimental results was good. For loads beyond that the nu- 
merical results overestimated the postbuckling stiffness. The plates of Noor, et al. 
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[37], like the plates of Starnes and Rouse [34], failed at the nodal lines in a shear 
failure mode. From the numerical results Noor, et. al. [37] found that the shear strain 
energy density was maximum at the nodal lines. In their analysis interlaminar shear 
stresses were calculated from the three-dimensional equilibrium equations of 
elasticity by using the inplane stresses obtained from the finite-element analysis. By 
integrating the equilibrium equation in the thickness direction, the interlaminar shear 
stresses were calculated. 

Shin [38] studied postbuckling failure of solid plates and attempted to find a 
weight-saving design by changing the stacking sequence so that the postbuckling 
failure load was maximized. He used the maximum strain failure criterion by con- 
sidering inplane strains only. The phenomenon of mode shape change has been ob- 
served by several investigators [39-42]. Some investigators [43,44] followed Stein s 
approach [39] in considering the mode shape change. Stein's approach [39] was to 
find the load-end shortening relation for successive buckling mode shapes and 
superpose them. Then, the broken continuous curve that forms the lower bound en- 
velope is assumed to be the load-end shortening relation of the plate. In the 
postbuckling analysis Shin [38] tried to predict mode shape change by examining the 
energy state for several successive mode shapes. Shin claimed that mode shape 
change occurs at the point where the energy level crosses, instead of at the crossing 
point of the load-end shortening relation. 

Engelstad, Reddy, and Knight [45] studied postbuckling failure numerically. They 
used a first-order shear deformable degenerated shell finite-element theory to predict 
the postbuckling response and failure. By applying a progressive damage failure 
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mechanism and comparing with Starnes and Rouse's experimental results [34], they 
obtained good agreement in failure load, failure mode, and failure location. 


1.3 Problem Studied 


1.3.1 Description of Problem 

The plate geometry, nomenclature, and a description of the problem studied are 
shown in Fig 1-1. The coordinate system used in the study is illustrated in the figure. 
The origin of the coordinates is at the center of the hole and the x axis is vertical. The 
positive z axis is out of the plane of the figure. The top edge and the bottom edge 
are clamped and the two side edges are simply supported. The top edge is loaded 
by a known uniform inplane displacement denoted as u. The bottom edge is re- 
strained to have zero displacement in the x direction. In numerical simulation the in- 
plane displacement in transverse direction is not restrained on both the simply 
supported and clamped edges. Figure 1-1 illustrates an important aspect of the 
problem. Specifically, in the experimental phase that accompanies the numerical 
analysis, the actual width of the plate, b, is wider than the distance between the 
simple supports, and the actual length of the plate, a, is longer than the distance be- 
tween the clamped supports. The width and length of the plate between the supports 
are denoted as W and L, respectively. The portions of the plate outside the dimen- 
sions W and L are necessary from the experimental point of view. To clamp the top 
and bottom edges part of the plate must be within the clamps. The clamps extend 
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inward some distance, yet the known inplane displacement, u, is applied at the outer 
most edge of the upper clamp, and the zero inplane displacement is enforced at the 
outermost edge of the lower clamp. Likewise, when the plate deflects out-of-plane 
due to buckling the simply supported edges move inward. Therefore, the plate must 
be wider than the distance between supports to accommodate this narrowing, or the 
plate will not stay in the fixture. In the experimental set-up knife edges are used to 
simulate simple supports. The dotted lines in Fig. 1-1 indicate the lines of support 
provided by the knife edges and clamped supports. The hole diameter is given by D 
and the hole is centrally located. The ratio of hole diameter to plate width is again 
given by D/W and attention must be given to the definition of W, as opposed to b, in 
Fig. 1-1. 


With the problem described, the statement of the goal of the work is simple: 
Determine what level of inplane load (or inplane displacement u! at the clamped 
boundaries causes the plate to fail. Inherent in this is : (1) the definition of failure and 
(2) the determination of the location of the failure and the mode of failure. 

1.3.2 Approach 


The approach used to study postbuckling failure involves both analysis in the 
form of numerical predictions, and experiments. To put the postbuckling analysis into 
context, the buckling problem is studied first. Then the postbuckling response, or 
more accurately stated, the geometrically nonlinear response, is computed 
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The numerical analysis for determining the buckling loads uses a mixed- 
formulation finite-element method following Reddy [46]. The analysis for studying 
response in the postbuckling range relies on a shear-deformable finite-element for- 
mulation following Reddy and Chao [47]. Imperfections in the plates in the form of 
initial out-of-plane deformations are being included. The inclusion of imperfections 
actually precludes the plate from exhibiting ideal postbuckling response. Thus, the 
study of postbuckling is done using a geometrically nonlinear analysis which predicts 
out-of-plane deformation response for even small inplane loads. Since the emphasis 
is on material failure, the calculation of stresses is included in the geometrically 
nonlinear analysis. Inplane stresses are computed directly from the analysis. The 
interlaminar stresses are calculated by integrating the three-dimensional equilibrium 
equations of elasticity in the thickness direction. Failure is predicted by using the 
maximum stress failure criterion. 

The actual responses of the plates are studied by using a specially built fixture 
designed to uniformly displace the top edge of the plate. The fixture is used within 
a displacement-controlled Tinius-Olsen load frame. Two linear variable displacement 
transducers (LVDT's) are used to measure the downward motion of the top edge. The 
load required to do this is measured with the load cell installed in the load frame. 
Hence load vs. end shortening relations from the experiment can be compared with 
similar data from the numerical analysis. The out-of-plane displacement near the 
center of the plates is measured with a third LVDT. This response can be compared 
with numerical prediction. Strain gauges at specific locations are used to measure 
the uniformity of the response and bending effects. The shadow moire method is 
used to provide a whole-field picture of the out-of-plane displacements. These 
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whole-field measurements are compared with contour plots generated by the nu- 
merical analysis. 

The issue of failure, the prime focus of the work, is studied primarily by visual 
inspection, and by considering the response of the plate as measured by the strain 
gauges and LVDT's. However, there was an attempt to study failure using acoustic 
emission. This approach was partially successful in detecting the onset of failure. 
The failure load, mode, and location, as determined by the instrumentation and by 
visual inspection, respectively, are correlated with predictions. 

The plates actually considered are fabricated from AS4/3502, a graphite-epoxy 
material. The plates are a = 11 in. by b = 11 in. square with 3 in. diameter hole. The 
distance between the simple supports is 10 in. and the distance between the inner 
edges of the clamped supports is 10 in.. Thus one-half inch on all four sides is used 
for enforcing the support conditions, and D/W = 0.3. Four laminates are considered. 
The laminates are [ ± 45/0/90]», [ ± 45/0*]*$, [ ± 45/0 6 ] s , and [ ± 45]«. All laminates 
are 16 layers. The first laminate represents a conventional design. The second lami- 
nate is also common and represents a departure from the quasi-isotropic nature of 
the first laminate. There is a preferred direction with the second laminate. The third 
laminate has a preferred direction, and represents, perhaps, an extreme case. The 
fourth laminate has no 0° layers and represents a laminate often studied. 

To follow are descriptions of each chapter of this study. The next chapter ad- 
dresses the issue of buckling. The finite-element formulation used is briefly over- 
viewed and the predicted buckling characteristics of the laminates considered in the 
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experiments are discussed. For contrast, the case of plates with no holes is exam- 
ined. Buckling loads and accompanying mode shapes are presented. 

Chapter 3 addresses the issue of postbuckling response. Since plates generally 
exhibit imperfections in the form of initial out-of-plane displacements, the chapter 
actually addresses the geometrically nonlinear response of plates which exhibit an 
initial out-of-plane deflection. The finite-element formulation to solve the governing 
equations is explained. The results from the formulation are compared with results 
from the general purpose code ABAQUS [48]. The postbuckling response of the four 
laminates with and without holes are studied using the formulation developed. The 
influence of the plate extensions outside the analysis region is considered in this 
chapter. 

Interlaminar stress calculations are addressed in chapter 4. The formulation for 
computing the stresses is presented. This formulation is based on integrating the 
equations of equilibrium from the theory of elasticity through the thickness. Since the 
finite-element analysis is an approximation, the interlaminar stress calculation is 
compared with other analyses to validate the numbers. After the validation 
interlaminar stresses in the plates in this study are examined. The stresses are 
computed at several locations along the net-section for the case of the plate 
postbuckling with two half-waves in the loading direction. Stresses at the hole edge, 
partway between the hole edge and support, and inside and outside the simple sup- 
ports are considered. 

Chapter 5 addresses the primary focus of this study, namely, the analysis of fail- 
ure. The maximum stress failure criterion is used on a layer-by-layer basis to study 
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failure. The failure load, failure mode, and failure location for each of the plates 
studied, both with and without holes, are predicted. A so-called damage accumulation 
scenario is depicted for each plate. Because of the aforementioned mode shape 
change, the failure analysis considers failure scenarios for both the one half-wave 
and two half-wave configurations of deformation in the loaded direction. 

Chapter 6 describes the test fixture and experimental set-up. Detail drawings of 
each part of the fixture are provided. The test procedures and the layout of the 
equipment are explained in this chapter. 

Experimental results and the comparison of experimental results with numerical 
results are presented in chapter 7. For each laminate the load vs. deflection relations, 
load vs. strain relations, shadow moire results, failure, and acoustic emission results 
are discussed. For the [ ± 45/0 6 ] s laminate, which showed a difference between ex- 
perimental results and numerical predictions, a re-examination of postbuckling failure 
is attempted, by revising the interlaminar shear strength to take into account the layer 
clustering effect. 

Chapter 8 provides concluding comments and recommendations for future re- 
search. 
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2.0 


Buckling Finite-Element Formulation 


2.1 Governing Equations 


In this chapter the buckling analysis using a finite-element method based on a 
mixed formulation will be explained. The formulation used in this chapter follows 
reference [46]. The basic equations used in this mixed finite-element formulation are 
the stability equation, 
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and the moment-curvature relations, 
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Recall, the coordinate system used was illustrated in Fig. 1-1. Here M x , M y , and M ^ 
are moment resultants defined by 
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where h is the laminate thickness. The symbol k is the buckling parameter which is 
defined as 
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where N x , N y , and are the inplane force resultants. These are defined by 
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The quantities n 10 , and n M are the inplane force resultants due to a unit load applied 
on the edge of the plate. Interpreted differently, ni 0 , n m , and n M are the prebuckling 
force resultants due to a unit applied load. The meaning of unit load is different for 
a uniform stress loading than it is for a uniform displacement loading. This difference 
in meaning will be explained in the numerical procedure section of this chapter. The 
quantity w is the out-of-plane displacement. The quantities D , are bending stiffness 
coefficients which are defined as 
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where Qii is the transformed reduced stiffness matrix. To facilitate the derivation of 
the variational form, the three moment-curvature relations are expressed in the al- 
ternative form as 
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(2.7) 


a 2 w 

dxdy 


(D^ M x + D 26 My 4- Dee M X y), 


where 


, 

pi 

0 12 

1 

<0 

IQ 


1 

p 

D 

fO 

Di6 

D i2 

d 22 

0 2 6 

= 

D i2 

D 22 

0 2 6 

Oi6 

d 26 

0 6 6 


Ol6 

0 2 6 

^66 


( 2 . 8 ) 


The variational form of eqs. (2.1) and (2.7) can be established by using the inverse 
procedure of deriving the governing equations. This procedure appears as 


6J 


-JH 


d 2 M, 


+ 2 


d 2 M 


xy 


dx * dx d y 
+ A. (n 10 


+ ■ 


d 2 w 


dx ‘ 


d 2 M y 

~d?~ 
+ n 20 


d 2 w 


~~T- + 2 Hqq . 

dy 2 5 y 


) ^ dxdy 
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Using integration by parts and rearranging the terms leads to 

r 


SJ = 5 


Y J {D„M 2 x + D 22 My + D 66 a4 + 2D, 2 M x M y 

* 

+ 2 D :6 M x M xy + 2D 26 MyM xy ) dxdy ] 



( 2 . 10 ) 


4- 8 


QnW + M * ^ 7 v * + My ^7 r,y + M ^( **)}* 


where tj x and r/ y are direction cosines, and Q n is defined as 


dw 


dw 


Q «-[ N * dx +Nxy dy 


+ ■ 


dM* 

dx 


dM 


+ • 


xy 


dy 


Mx 


( 2 . 11 ) 


+ 



+ Ny 


dw 

dy 


dMxy 

dx 


dMy 

dy 


n y . 


The direction cosines, force resultants, and moment resultants are shown in Fig. 2-1. 


Therefore, the variational functional can be expressed as 
J(W , M X ' ft/fy j M x y) 



{D^M X + D 22 My -I- D 6e M xy + 2 D, 2 M x My + 2D, 6 M x M xy + 2D 26 M y M xy }dxdy 
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Fig. 2-1 Force resultants, moment resultants, and direction cosines. 
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i" e^ + 3M 2 _ ( ^L + ^p_ )l rfxrfy 

dx ' 5x dy ' 5y ' dx dy f 


+ 


U 


dw , dw . dw dw 

n «17 + n ’°17 +2 "“ 1717 r xrfy 


( 2 . 12 ) 


» 

{", f - «. + **,-«=■ », + <V ^ + -f=- *> + Q„*}*, 


Interpolating the variables w, M„, M y , and M „ in terms of nodal values of these vari- 
ables, and constructing algebraic equations, leads to an eigenvalue problem which 
will be explained in the following section. From the above procedure to derive a 
functional it can be observed that the essential variables are 


w, M x , M y , M xy , 


( 2 . 13 ) 


and the natural variables are 


0 


n, 


dw 

dx Vx ’ 


dw 

W' 1 * 


, dw 

( 17 , >' + 


dw 


dy 


Vx)- 


( 2 . 14 ) 


Note that the natural variables involving the bending moments in the displacement 
formulation become the essential variables in the mixed formulation. 
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2.2 Finite-Element Formulation 


The quantities w, M„, M y , M „ can be interpolated by expressions of the form 



(2.15) 


where <f>i are appropriate interpolation functions, and w, , M xt , M yt , and are nodal 
values of out-of-plane displacement and moment resultants to be solved for. These 
constants are the nodal values of w, M x , M y , M respectively. By substituting eq. 

(2.15) into eq. (2.12) and equating the variations of J with respect to 
Wj, Mu, M yi , to zero separately, a matrix equation of the following form is ob- 


tained : 


k’ 

K- 12 

K ij 

fc' 13 

K n 

1 

w 



0 

0 

0 

w 

CM 

> T— I^i 

CM 

. CM 1^. 

* 

ts 23 

K ij 

K 24 

Ktj 

M x 

= A 

0 

0 

0 

0 

M x 

13 

K ii 

U 23 
K 'j 

v-33 

*ij 

U' 34 

Kij 

My 

/I 

0 

0 

0 

0 

My 

1 

* 

C: 

K 24 

K ij 

K 34 

K 'J 

K 44 

Kij 

_M X y 


0 

0 

0 

0 

_M X y_ 


where 


K,; 1 - 0. 


L f 12_ 

N - 


d(j>i d<f)j 
dx dx 


dx dy 


/% 

i/13 

Kfj = 

4 


50/ d<t>j 
dy dy 


dx dy, 


C= l( 


d<t>i d<t>j dtfij d<t>j 
dx dy + dy dx 


)dx dy 


(2.16) 
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K 


24 

I] 


K 


34 


D 16 0/ 4>jdx dy. 

K/f = - Jd 22 0/ 4>jdx dy 

D 2 6 0 1 <t>jdx dy. 

K“ = - ju^^^jdxdy 


(2.17) 


M 


ii 

i I 


d<t>, 
ni ° dx 


d(f>j d<t>i d(f)j 50 / d<t>j 

5x /?2 ° 5y 5y + neo dx dy 


d<t>i 

dy 



dx dy. 


As mentioned in the previous section, moments are essential boundary conditions in 
this mixed formulation. Therefore these essential boundary conditions must be sat- 
isfied exactly. On the hole boundary the moments must satisfy the conditions 


M n = 0, M ns = 0, 


(2.18) 


where, as shown in Fig. 2-1, M n is the bending moment in the normal direction on the 
hole boundary, and M„, is the twisting moment on the hole boundary. However, the 
nodal values being solved for are the moments in the x-y coordinate system. In order 
to satisfy eq. (2.18) on the hole boundary, M„ and M n , need to be expressed in terms 
of the moments in the x-y coordinate system. Such an expression can be obtained 
from the moment equilibrium equation applied to the differential triangular element 
shown in Fig. 2-1. The relations are expressed as follows: 

M n = M X t]x + My t]y + 2 M X y t] tj y 

M ns = (My ~ M x ) r\ x rjy + /W xy (^ - »/J). (2.19) 

Since M„ and M n , are zero on the hole boundary, from eq. (2.19) the following relations 
are valid on the hole boundary : 
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( 2 . 20 ) 


M, 


Vy 


M X y=~(^)M X . 


Substituting these relations for the nodal values on the hole boundary in eq. (2.15) 
satisfies the free-edge condition on the hole boundary. 


2.3 Numerical Procedure and Results 


As stated earlier, the ni 0 . n 20 . and n w in the stability equation, eq. (2.1), are the 
prebuckling force resultants which must be calculated before the buckling analysis. 
These force resultants are calculated at Gauss points so that they can be used in the 
numerical integration procedure used in conjunction with the last of eq. (2.17). The 
prebuckling force resultant calculation is a two-dimensional planar analysis. A dis- 
cussion of this step is not included here because it can be found in any one of a 
number of books discussing finite elements. For example, refer to [48] for more de- 
tails on the analysis procedure and the computer program "FEM2D". 

For a uniform stress loading if a unit force resultant is applied on the edge, i.e., 
N x = *\ on the edge, then the resulting eigenvalue X represents the critical value of 
N x . For a uniform displacement loading, an arbitrary but uniform displacement is ap- 
plied on the boundary. Then the force resultant distributions in the plate and average 
force resultant on the displaced edge are calculated. The average force resultant on 
the displaced edge is defined as 
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( 2 . 21 ) 


w_ 

*> = irL N ’ dy ' 

2 

where W represents the width of the plate. It needs to be mentioned that if the anal- 
ysis includes the region outside the supports, b is used instead of W (refer to Fig. 1-1). 
The force resultants within the plate are divided by this average force resultant on the 
displaced edge so distributions of the force resultants due to a unit average force 
resultant are obtained. These force resultant distributions due to a unit average force 
resultant are the n 10 , and n w in the notation used here and these resultants are 
substituted into the last of eq. (2.17). The resulting eigenvalue equation is eq. (2.16). 
Again, the resulting eigenvalue k represents the critical applied average force re- 
sultant that results from a uniformly displaced edge. 

For verification of the computer program the buckling coefficients were com- 
pared with those of Nemeth [31] for symmetric [±0 ]k angle-ply laminates under a 
displacement loaded, clamped boundary condition. When the plate has the Du and 
D w components in the bending stiffness matrix, even though the geometry and load- 
ing are symmetric with respect to the x and y axes, the out-of-plane deflection is not 
symmetric with respect to these axes. In this case a full plate analysis must be used. 
Strictly speaking, a quarter-plate analysis is not correct. However, Nemeth [21] 
showed that if the anisotropic coefficients defined as 

y <5 = r (2-22) 

(Di 3 iD 22 )4 (D 11 D 2 3 2 )4 
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are less than 0.2, a result very close to a full-plate analysis can be obtained with a 
quarter-plate analysis with Du and D 28 set to zero. Nemeth used this approach and 
obtained results from Kantorovich's method. A comparison between Nemeth's re- 
sults and the present analysis is shown in Fig. 2-2. The buckling coefficient in this 
comparison is defined as 


N? W 2 
*VDiiD 22 


(2.23) 


For fiber angles of 0° , 30° , and 45° and changing D/W from 0.1 to 0.5 the results are 
compared. Fig. 2-2 shows that the present analysis is in good agreement with 
Nemeth's result. The present analysis considered one-quarter of the plate and as- 
sumed D 16 = Dj« = 0. Nemeth [31] also assumed D,« = D» = 0. 


Considering the plates in the postbuckling study, for a buckling analysis with the 
finite-element method large computer memory capacity is necessary. In the buckling 
analysis of a plate with a hole the memory size can be a limitation. The lowest 
buckling load, and the mode shape corresponding to this load, can generally be ob- 
tained with a small number of finite elements. However, it is possible that some of the 
other low buckling modes would not be detected with a full-plate analysis when using 
a small number of finite elements. In order to avoid this possibility of missing a 
buckling mode, yet staying within the limits of the computer capacity, based on 
Nemeth's argument a quarter-plate analysis is used in this study. However, a variety 
of boundary conditions, discussed below, must be used to account for the higher 
modes. Material properties for AS4/3502 used in the present study are tabulated in 
Table 2-1. With these material properties the anisotropic coefficients of the laminates 
selected for study, namely [ ± 45/0/90] 2S , [ ± 45 /OJk , [ ± 45/0 8 ] s , and [ ± 45]4s lam- 
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Nemeth Present Nemeth Present Nemeth Present 
Odeg. 0 (teg. 30 deg. 30 deg 45 _^®9 45 jjeg. 



Fig. 2-2 Comparison with Nemeth's angle-ply results [31] for displacement loaded, clamped 
[ ± *]« laminate. 
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Table 2*1 Layer material properties of AS4/3502. 


E 1 

18.500 

Msi 

E 2 

1.600 

Msi 

G 12 

0.832 

Msi 

G 23 

0.332 

Msi 

V 12 

0.350 

thickness 

0.00531 

in. 
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inates, are shown in Table 2-2. Since the coefficients are considerably less than 0.2, 
a quarter-plate analysis can be used with sufficient accuracy. 

For a complete study with a quarter-plate analysis, to examine more than one of 
the lower modes, four kinds of boundary conditions must be considered on the x and 
y axes. These boundary conditions are explained in Fig. 2-3(a)-(d). The first case is 
free on the x axis and free on the y axis. This case produces mode shapes that are 
symmetric with respect to both the x and y axes. The second case is free on the x 
axis and simply supported on the y axis. This case produces mode shapes that are 
symmetric with respect to the x axis and anti-symmetric with respect to the y axis. 
The third case is opposite to the second case. The fourth case is simply supported 
on both the x and y axes, and this case produces mode shapes that are anti- 
symmetric with respect to both the x and y axes. Of course, in all of these cases the 
clamped and simply supported edges remain that way. 

From all of the buckling loads calculated with the four sets of boundary condi- 
tions, the smallest three values are selected for each plate and they are listed in Ta- 
ble 2-3 for plates with holes. As can be seen as the number of 0° layers is increased, 
the buckling load increases. The [ ± 45]«s lamimate shows the smallest buckling load. 
For future use the critical (buckling) force resultants are nondimensionalized in the 

p 

form of , where P cr is the buckling load ( N„ cr multiplied by load applied length), 
E is the equivalent Young's modulus in the x-direction of the corresponding laminate, 
and A is the area of the edge of the laminate where the load is applied. The equiv- 
alent Young's modulus of each laminate is tabulated in Table 2-4. 
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Table 2-2 Anisotropic coefficients for laminates in this study. 


Laminate 

y 

5 

[ ± 45/0/90] 2S 

0.067 

0.073 

[ ± 45/0 2 ] 2S 

0.062 

0.092 

C ± 45/0 6 ] s 

0.041 

0.068 

[ ± 45] 4S 

0.134 

0.134 


Buckling Finite Element Formulation 


35 






Table 2-3 Buckling loads of plates with holes. 


Laminate 

Parameter 

First mode 

Second mode 

Third mode 

[ ± 45/0/90] 2S 

Nx r 

258 

424 

470 

Pci 

EA 

0.000411 

0.000676 

0.000749 

[ ± 45/0 2 ] 2S 

M Cr 

Nx 

294 

496 

497 

Prr 

EA 

0.000321 

0.000541 

0.000542 

[ ± 45/0 s ]s 

Nx cr 

318 

508 

545 

Per 

EA 

0.000245 

0.000407 

0.000436 

[ ± 45]« 

Nx* 

240 

373 

427 

-Eg 

EA 

0.000980 

0.001523 

0.001 744 


(unit of FTf* = lb/in) 
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Table 2-4 Equivalent Young's modulus In x direction for each laminate. 


[ ± 45/0/90] 2S 

7.381 E+06 psi 

[ ± 45/0 2 ] 2S 

1 .078E+07 psi 

[ ± 45/0 6 ] s 

1 .469E+07 psi 

[ ± 45] 4S 

2.881 E+06 psi 
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In conjunction with the buckling loads, the buckling mode shape must be under- 
stood. From Fig. 2-4 and Fig. 2-5 it can be seen that for every plate the first mode is 
one half-wave in both the x and y directions. For the second mode the [ ± 45/0/90] 2S 
and [ ± 45]« plates have two half-waves in the x direction and one half-wave in the 
y direction. However, for the second mode the [ ± 45/0*]*$ and [ ± 45/0 8 ] s plates have 
one half-wave in the x direction and two half-waves in the y direction. In the third 
mode the [ ± 45/0/90] 2S and [ ± 45]^ plates have two half-waves in the y direction and 
one half-wave in the x direction, but the [ + 45/0 2 ] 2S and [±45/0 6 ] s plates have one 
half-wave in the y direction and two half-waves in the x direction. It is important to 
note from Table 2-3 how close together the second and third mode buckling loads can 
be. In an actual experiment it might be impossible to distinguish between these two 
modes. 

For comparison, the buckling loads of solid plates of the same lamination as the 
plates with holes were calculated. The smallest three buckling loads for each plate 
are tabulated in Table 2-5, and the buckling mode shapes are shown in Fig. 2-6 and 
Fig. 2-7. The buckling loads of solid plates are higher than those of the plates with 
holes. All the solid plates have the same mode shapes in the first and second 
modes. The first mode shape is one half-wave in both the x and y directions, and the 
second mode shape is two half-waves in the x direction and one half-wave in the y 
direction. In the third mode the [±45/0/90]^ and [±45]« laminates showed two 
half-waves in the x direction and one half-wave in the y direction. The [ ± 45/0 2 ] 2S and 
[±45/0»]s plates showed one half-wave in the x direction and two half-waves in the 
y direction. For the case of plates without holes the buckling loads of the second and 
third modes are far apart. With the buckling analysis established, and numerical re- 
sults computed, the next chapter addresses the issue of postbuckling. 
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Second mode 



Third mode 


Fig. 2-4 Buckling mode shapes of [±45/0/90]j, and [±45]«j plates with holes. 
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First mode 



Second mode 



Third mode 


Fig. 2-5 Buckling mode shapes of [ ± 45/0 ? ] >s and [ ± 45/Ot], plates with holes. 
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Table 2-5 Buckling loads of solid plates. 


Laminate 

Parameter 

First mode 

Second mode 

Third mode 

[ ± 45/0/90] 2S 

R? 

293 

444 

762 

_£a 

EA 

0.000467 

0.000708 

0.001215 

[ i 45/0 2 J 2 s 

r;or 

Nx 

327 

528 

888 

Per 

EA 

0.000357 

0.000576 

0.000969 

[ ± 45/0 6 ] s 

Nx 

340 

571 

841 

Per 

EA 

0.000272 

0.000457 

0.000673 

[ ± 45] <s 

Nx 

291 

405 

650 

Per 

EA 

0.001188 

0.001654 

0.002654 


(unit of Nx cr = Ib/in) 
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First mode 



Second mode 



Third mode 

Fig. 2-6 Buckling mode shapes of solid [ ± 45/0/90 ]m ®nd [ ± 45]4s plates. 
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3.0 Nonlinear Plate Bending With Initial 

Imperfection 


For analysis in the postbuckling range of deflection, a plate bending finite- 
element program which includes the influence of geometric nonlinearity is devel- 
oped. For an ideal plate the stiffness matrix of the algebraic equations becomes 
singular at the buckling load. Therefore, in postbuckling studies numerical difficulties 
can occur at load levels near the buckling load. One of the ways to overcome this 
singularity problem at the buckling load is to impose a small initial out-of-plane de- 
flection on the plate. As plates fabricated to be flat are rarely without an initial out- 
of-plane deflection, the inclusion of these effects in the analysis is quite realistic. In 
addition, numerical difficulties are avoided. This is the approach taken here. 

The basic kinematic equations for plate bending with initial imperfection can be 
found in reference [50]. The finite-element formulation procedure for this develop- 
ment, in the context of a first-order shear deformation theory, will be explained in the 
following sections. The formulation results in a set of nonlinear algebraic equations. 
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The resulting nonlinear algebraic equations are solved by the Newton-Raphson 
method. The general features of the Newton-Raphson method will also be explained 
in a following section. 


3.1 Governing Equations 


For a first-order shear deformation theory, the displacement field is assumed to be 
of the form 


u,(x, y, z) = u{x, y) + z ^ X (x, y) 

u 2 (x, y, z) = v{x, y) + z \f/ y (x ,y) (3.1) 


u 3 (x, y, z) = w(x, y). 


Here u 1f t/ 2 , u 3 are the displacements in the x, y, z directions, respectively, u, v, w are 
the associated midplane displacements, and and \j/y are cross-section rotations of 
the line elements perpendicular to the midplane in the x-z and y-z planes, respec- 
tively. The strains in the plate are expressed as 


o . „ o 

fix = £ x + Z K x 


Cy = + Z 


0,0 

Vxy — yxy "I" z K xy 


(3.2) 
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, . dw 

y « = ^ + 17 


/ I UVV 

Vyr ~'l'y+ Qy • 

where the eS , ej , and yjy are the midplane strains, and the kS , kJ , and kS,, are the 
midplane curvatures. By one of the assumptions of plate theory, lengths normal to 
the midplane remains unchanged. Therefore the strain in the thickness direction is 
zero, i.e., c z = 0. The midplane strain-displacement relations and midplane curva- 
tures for a plate with an initial imperfection can be written explicitly as 

o du 1 / dw \ 2 . dw 

C *~ dx 2 \ dx ) dx dx 

2 

o _ dv _1_ / dw \ dw 

c y~ dy 2 \ dy ) dy dy 


o du dy dw dw dw dw 

V* y dy dx dx dy dx dy dy dx 



(3-3) 



0 Wx d jtl 

K * y dy + dx 


where the influence of geometric nonlinearities has been included and £ is the initial 
out-of-plane imperfection. 
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The governing equilibrium equations can be derived by the principle of virtual 


work. That principle can be expressed as 


0= 5U-8W, 


(3.4) 


where 0 is the internal strain energy and IV is the work done by external forces. The 
body force terms are neglected in this equation. The variation of the work done by 
the external forces is expressed as 


6VJ = J<7 <5 w dx dy + J ( N n du n - M n 


dSw ±N ns du s -M ns -^+Q n 5w)dS, 


dn ' " ns '" s '" s ds 

where N n , Nn,. Mn, M„,. and Q„ are the specified applied resultants defined by 


(3.5) 


JL -£L 

A f 2 A A f 2 A 

N n = J N ns = J ^ t ns dz, 


(3.6) 


h JL -IL 

A f2 A A f 2 A A f 2 J. 

A4 n = J ■ zt n dz , M ns = J ^ zt ns dz , Q n — j ^t nz dz. 


In the above the effect of a transverse load, q, is included for completeness. Here 
t„. L. and L are specified normal and tangential stress components, and u„ and u, 
are normal and tangential components of the midplane displacements on the edge. 
Now, including the strain energy term, the principle of virtual work is expressed as 

0 = f (<r x Se x + < 7 y 8c y + t xy 5y xy + t y Z 5 y yz + t xz Sy xz ) dV (3.7) 


J 


q 5w dx dy — 


N n 5u n dS — 


ns S u s dS 


A 

M n 

% 


ddw 

dn 


dS + 


A 

M n 


ddw 

ds 


dS 


Q n 8w dS. 
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To derive the equilibrium equations, in addition to the moment and inplane force re- 
sultants defined in eq (2.3) and eq (2.5), respectively, the out-of-plane force resultants 
Q x and Q y defined as 




dz. 


(3.8) 


are used. By substituting the strain-displacement relation, eq. (3.3), into the virtual 
work statement, eq. (3.7), integrating in the thickness direction, applying integration 
by parts, and using the relations 

N n — N x Yf X + Ny T}y + 2 N X y Tj % t] y 

N ns = (Ny N x ) y\ x t]y + N X y ( t] x f]y) 

M n = M X Tjl+ MyTI 2 y + 2 M X y Y\ y ( 3 ‘ 9 ) 

M ns = ( My — M X ) T] X t]y + M X y ( X ~ y) 


and 


U n = Ut} x + Vtly, U S = -Ut]y + VT] X , 

u = u n ri x -u s riy, v = u n t]y+ u s rj x , 

•An = iMx + <My- <As = ~ 'I'xVy + ^yVx> 
<Ax= i'nVx- 'I'sVr 'l / y = 'l / nVy+ 'P$Vx< 


_d d d 

dn ~ * tx dx ^ y dy 

_d d d 

8x ~ r,x dn Vy ds ’ 


d 


d 

d 

ds 

— Vx 

dy ‘ 

Vy dx 

d 


d 


dy 

= *7y 

dn 

ds 


( 3 . 10 ) 


where ti» and rj y are direction cosines as shown in Fig. 2-1, the following form of the 
principle of virtual work is obtained: 
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(3.11) 


+ dw 


dQx dQ y _d_ 

dx dy dx 



dw_ 

dx 


+ 



H 

dy 


+ 


d_ 

dy 



dw 

dx 


H 

dx 



H 

dy 



dx dy 



+ J<K (N„ - N n )dS + j"<5u s (N„ - N ns )dS + JW„ (M n - M n )dS 
+ JW, K s - MnsjdS + J<5w (<?„ - Q n )dS. 


Since the variations du, dv, dw, di// x , d\f/ y , du„, du, t d\f/„, and dip, 
equilibrium equations are as follows : 



The boundary conditions are : 
either du n = 0 , or N n — N n — 0 
either du s = 0 , or N ns -N ns = 0 

A 

either d^ n = 0 , or M n — M n = 0 

A 

either d\j/ s = 0 , or M ns — M ns = 0 



are arbitrary, the 


(3.12) 


(3.13) 
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either <5w = 0 , or Q n — Q n = 0 


3.2 Finite-Element Formulation 

In this work a displacement based finite-element formulation for the governing 
equations derived in the previous section is used. The principle of virtual work 
statement in the previous section , eq. (3.11), can be expressed in terms of displace- 
ments u, v, w, \p K , and using Hooke's law and the strain-displacement equations. 
Over a typical element the generalized displacement U, where U represents u, v, w, 
\j/ x , and \l/ y , is interpolated spatially by an expression of the form 

r 

v = ^ u i<i>ix<y)< ( 3 - 14 ) 

i 

where U, is the value of U at node i, 4> , is the finite-element interpolation function at 
node i, and r is the number of nodes in the element. For simplicity in this formulation 
the same interpolation function is used for each of the five deformations u, v, w, 
and i j/ y . Substituting the interpolated expression for the displacements into the prin- 
ciple of virtual work statement, eq. (3.11), and collecting the terms of 3 U< leads to the 
following equation for a typical element: 

[K e (A)] {A} 8 = {F}°. (3.15) 

Here {A}* is the column vector of the nodal values of the generalized displacements, 
[K e ( A)] is the matrix of the stiffness coefficients, and {F}* is the column vector con- 
taining the boundary and transverse force contributions. The notation [K e (A)] indi- 
cates that the stiffness matrix is a function of the nodal values. Thus the problem is 
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nonlinear in the nodal displacements. The assembled finite-element equations are 
of the form 


[K(A)] (A) = {F}, (3.16) 

where [K(A)] is the direct stiffness matrix and {A} is the global solution vector. This 
nonlinear set of equations is solved iteratively. Here the Newton-Raphson method is 
used, and it is illustrated in the next section. 


3.3 Newton-Raphson Method 


To solve the system of nonlinear equations, eq. (3.16), they are transformed into 
the form 


[K r (A r )] {<5A r } = {*}, 


(3.17) 


where {A r } is the intermediate solution at the r-th iteration and [K T (A')] is the tangent 
stiffness matrix which is obtained from the relation 


[K T ( A r )] = 


d{A r } 


(3.18) 


{ R } is the residual given by 

{R} = [K(A r )] {A r > - {F}, (3.19) 

and (<5A'} is the increment of the solution. The tangent stiffness matrix in eq. (3.18) 
is symmetric. The total displacement at the (r + 1)-st iteration becomes 


Nonlinear Plate Bending With Initial Imperfection 


52 



{A r+1 } = {A r } + {M f }. 


(3.20) 


The true solution makes the residual zero, i. e. 


{R} = o. 


(3.21) 


In this implementation iteration is stopped when the following equation is satisfied: 




/ = i 


^EPS, 


(3.22) 




i — 1 


where the value of EPS is 0.01 in most of the work. 


3.4 Numerical Results 

3.4.1 Verification of Present Analysis 

Based on the argument associated with Table 2-2, a quarter-plate model is again 
employed in the postbuckling analysis and Di 6 and D w are assumed to be zero. A 
uniform inplane displacement is applied on one edge. Recall that the displaced edge 
and the edge opposite to it are clamped, and the two other edges are simply sup- 
ported. For the initial out-of-plane deflection a one-term double cosine function is 
used which is of the form 
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i = w o h cos ( ) cos ( )> 


(3.23) 


where h is the plate thickness, and L and W are the plate length and width inside the 
supports, respectively, as shown in Fig. 1-1. In the above equation, w 0 is the fraction 
of initial out-of-plane deflection with respect to the plate thickness at the center of the 
plate. In these results w 0 was taken to be 0.1. As mentioned previously, the nonlin- 
ear simultaneous equations are solved by the Newton-Raphson method, an incre- 
mental scheme using the solution of the previous load step as the initial guess of the 
next load step. 

For the verification of the computer program, postbuckling results were com- 
pared with the widely used commercial computer code ABAQUS [48], The case of a 
solid quasi-isotropic plate, a [±45/0/90] K laminate, was considered. The material 
properties used in the comparison are given in Table 2-1. The same boundary con- 
ditions as in the present study, i.e., clamped on the loaded edge and simply sup- 
ported on the other two edges, were applied. The dimensions of the plate were also 
10" by 10". In ABAQUS [48] the initial imperfect out-of-plane deflection was assumed 
to be the same as the first buckling mode. This is not quite the same situation as eq. 
(3.23), but it is close enough. Fig. 3-1 shows the comparison for the load vs. end 
shortening relation, and Fig. 3-2 shows the comparison for the load vs. out-of-plane 
deflection relation. The out-of-plane deflection being considered in Fig. 3-2 is the 
deflection at the center of the plate. Note that the nondimensionalized load parame- 

p 

ter, -=-r- , is the same as in chapter 2. For both the ABAQUS results and the present 
results 16 elements were used for a quarter-plate. As can be seen, the comparison 
of the present formulation with the commercially available code is quite good. Note 
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Fig. 3-1 
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Comparison with ABAQUS [48] for load vs. end shortening relation of solid 
[ ± 45/0/90] JS plate under displacement loaded, clamped condition. 
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Fig. 3-2 
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Comparison with ABAQUS [48] for load vs. out-of-plane deflection relation of solid 
[ ± 45/0/90] 2 s plate under displacement loaded, clamped condition. 
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that the load vs. end shortening relation is almost linear until the buckling load is 
reached. After that load level the slope of the load-end shortening relation decreases 
sharply. The secondary portion of the relation is associated with the postbuckling 
response of the plate. The load vs. out-of-plane deflection is also nearly linear until 
the buckling load is reached. The deflection then increases rapidly with increasing 
load, the region of high deflection again being associated with the postbuckling re- 
sponse. 


3.4.2 The Effect of the Portion of the Plate Outside the Supports 


As mentioned earlier, in any real experimental set-up the plate extends outside 
the supports. It is physically impossible to perform the experiment any other way. 
Thus it is important to know the influence on the plate response of the portions of the 
plate beyond the supports, and to determine if they must be included in the analysis 
of the plate. Here the effects of the extensions on the postbuckling analysis will be 
studied. To that end, the three cases shown in Fig. 3-3 are studied. In Fig. 3-3 the 
'no extension' case represents the usual numerical analysis where only the region 
inside the supports is considered. The 'simple extension' case represents the case 
where the portion of the plate beyond the simple support is taken into account. The 
analysis region is the region of 'no extension' plus the region beyond the simple 
support. The 'both extensions' case represents the condition where the entire region 
of the plate is considered. In this case the analysis region is the same as the speci- 
men in the test set-up. For each case in Fig. 3-3 the analysis region is enclosed with 
a solid line. The postbuckling response for the quasi-isotropic plate with central hole 
was examined for the above three cases. Figure 3-4 shows the result of the exten- 
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(c) Both extensions 

Fig. 3-3 Plate models to account for the influence of the region outside the supports. 
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Fig. 3-4 Influence of the region outside the supports on the load vs. end shortening relation. 
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sion study. The load vs. end shortening relations for the three cases are considered. 
Note the horizontal axis reflects whether the length within the clamps, L, or the actual 
length of the plate, a, is being considered. As can be seen from the figure, there is a 
significant difference in postbuckling response between the case of 'no extension' 
and the other two cases. By taking into account the extensions the postbuckling re- 
sponse becomes stiffer. There is some difference between the case of 'simple ex- 
tension' and 'both extensions'. However, there is much more difference between the 
case of 'simple extension' and 'no extension'. The case of 'both extensions' is 
somewhat stiffer than the case of 'simple extension'. It can be concluded that the 
extensions cannot be neglected in postbuckling analysis. In the subsequent dis- 
cussions in this study the actual dimensions of the test specimen will be analyzed, 
i.e., the case of 'both extensions'. 

3.4.3 Results for Plates With and Without Holes 

To provide insight into the types of response that might be expected from the 
plates to be tested, the load vs. end shortening relations for the four laminates are 
illustrated in Fig. 3-5 - Fig. 3-6. Figure 3-5 shows results for the case of plates with 
a hole, D/W = 0.3, while Fig. 3-6 shows results for plates with no hole. In the figures 
the load is again nondimensionalized to be in the form of P/EA. The value of E is 
different for each plate, the values having been listed in Table 2-4. The results in Fig. 
3-5 - Fig. 3-6 were obtained for the imperfection given in eq. (3.23). Though not ob- 
vious from Fig. 3-5, the [±45/0 e ] s laminate is the stiffest, followed by the 
[± 45/0 2 ] 2S , the [ ± 45/0/90] 2S , and then the [±45]4s laminate. To be expected, a 
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Fig. 3-6 Load vs. end shortening relations for solid plates. 
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comparison between Figs. 3-5 and 3-6 shows the plate without a hole to be slightly 
stiffer than the counterpart plates with a hole. 

The load vs. midplane deflection relations for the same cases are illustrated in 
Figs. 3-7 and 3-8. The resistance to out-of-plane deflections as a function of lami- 
nation sequence follows the same pattern as for the end shortening deflections, the 
[±45/0 e ]s being the stiffest. Though the load vs. deflection reiations in Figs. 3-7 and 
3-8 seem to indicate that certain laminates will reach higher loads than others, the 
issue of failure must be superposed on those relations. The key question is at what 
load level on the postbuckling portion of the response will failure occur. The calcu- 
lation of stresses and the prediction of failure will be discussed in the next two 
chapters. However, before those topics are addressed, another issue regarding the 
postbuckling response will be discussed. This follows. 


3.4.4 Mode Changes During Loading 


As was discussed in the introduction, it has often been observed that as lami- 
nates are loaded into the postbuckling range, there can be an unexpected and sud- 
den change of the postbuckled configuration [38-41]. If the laminate is responding in 
the postbuckling range by deflecting in a shape similiar to the first buckling mode, 
with a single half-wave in the loading direction, the laminate may suddenly change 
to a shape associated with the second or higher buckling mode characterized by 
multiple half-waves in the x direction. Whether or not such a mode shape change 
occurs is a function of the laminate, the geometry, the boundary conditions, and 
though hard to quantify, possibly any eccentricities or imperfections in the loading 
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or plate itself. Whether the plate is responding in a shape associated with the first 
mode or a shape associated with a higher mode can have a profound effect on the 
failure load and failure mode of the plate. Typically, if the plate responds with multi- 
ple wavelengths in the loading direction, interlaminar shear may dominate the stress 
state. On the other hand, if the plate responds with a single half-wave, inplane 
stresses may be important. Figures 3-9 through 3-12 show the load vs. end shorten- 
ing responses of the plates in the study assuming the imperfection in the plate is in 
the form of two half-waves in the loading direction, i.e., 

{ - w 0 /jsin(-^)cos(-^-). (3-24) 

For comparision the load vs. end shortening response for each plate with the initial 
one half-wave imperfection, eq. (3.23), is also shown. Plates with and without holes 
are considered for each laminate. As can be seen, in all cases the load vs. end 
shortening relation for the two half-wave case coincides with the one half-wave case 
until the load reaches the first buckling load. At loads greater than this level the load 
vs. end shortening relation for the one half-wave case turns toward a postbuckling 
path, while the two half-wave case continues as an extension of the prebuckling load 
vs. end shortening relation. At a load level corresponding to the buckling load of two 
half-waves in loading direction, the load vs. end shortening relation of the two half- 
wave case then turns toward another postbuckling path. The two postbuckling paths 
cross and usually at loads greater than the crossing load a sudden change in config- 
uration occurs. It is important to note that the figures with both the one half-wave 
imperfection and two half-wave imperfection do not represent the complicated 
phenomenom of the sudden change in postbuckled shape. That phenomenon will not 
be studied analytically here. There will be no attempts at predicting the load path that 
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leads to the plate being in the two half-wave configuration. However, failure in the 
postbuckling range due to high interlaminar shear stresses will be studied. That will 
be addressed by examining the interlaminar shear stresses when the plate is 
postbuckled in the shape associated with the two half-wave configuration in the 
loading direction. It should be mentioned that the deflection at the center of the plate 
is zero when the plate is responding in the two half-wave configuration. Thus there 
are no load vs. out-of-plane deflection relations shown for that case. 

With the analysis to predict the geometrically nonlinear response in hand, atten- 
tion shifts to the calculation of stresses associated with this response. The calculation 
of inplane stresses is trivial and is not discussed explicitly. The calculation of the 
interlaminar stresses, however, is more difficult. That is discussed in the next chap- 
ter. 
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4.0 Calculation of Interlaminar Shear Stresses 


4.1 Formulation 


Though understanding the deflection response in the postbuckling range of de- 
formation is important for the study of failure, it is absolutely essential to understand 
the state of stress within the plate. The calculation of the inplane stresses within the 
laminate is relatively straightforward. None of these stresses is generally zero. How- 
ever, in classical plate theory the out-of-plane shear stresses are zero by the as- 
sumption that the line element perpendicular to the midplane surface remains 
normal. The out-of-plane normal stress is also assumed to be zero. In a first-order 
shear deformation theory the out-of-plane shear stresses are constant through the 
thickness, but the out-of-plane normal stress is still assumed to be zero. In reality 
there exist out-of-plane stresses that vary from zero to some maximum value as the 
thickness of the laminate is traversed. This has been shown by Pagano [51]. Using 
an example of cylindrical bending Pagano obtained the distribution of interlaminar 
stresses through the thickness by using an elasticity solution. Pagano compared this 
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elasticity solution, which involves rigorous mathematical calculations, with another 
solution which was obtained from a simpler approach. This simple approach used 
the two-dimensional solution obtained from plate theory. By substituting the inplane 
stress distribution from the two-dimensional solution into the three-dimensional 
elasticity equilibrium equations and integrating these equations through the thick- 
ness, the interlaminar stresses were obtained. This simpler approach showed very 
good agreement with the elasticity solution. As a result of the calculations these 
interlaminar, or out-of-plane, stresses were shown to be small. However, even 
though the out-of-plane stresses may be small, the weakness of the material makes 
the knowledge of the character of these stresses important. As is well known a 
composite material is quite poor in resisting all stresses except for the fiber direction 
normal stress. In this chapter general expressions for the interlaminar shear 
stresses are developed, the interlaminar stresses being calculated from the inplane 
solutions of plate theory. The approach follows that of Pagano. With this approach 
the calculation of interlaminar normal stresses involves derivatives of displacements 
and rotations that are of a higher order than they are in the calculation of the 
interlaminar shear stresses. With the discretized results obtained from the finite- 
element solution, it is very difficult to obtain a converged solution for the higher-order 
derivatives, and hence a converged solution for the interlaminar normal stress. For- 
tunately, in the postbuckling problem there is no transverse load. Therefore it can 
be assumed that the interlaminar normal stress is not important. Because of this as- 
sumption the interlaminar normal stress is not included in this study. However, there 
can still be convergence problems with interlaminar shear stress calculations and 
hence convergence must be studied. This will be discussed shortly. 
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Excluding the interlaminar normal stress, o z , the two equilibrium equations of 


importance are 


d °x 

dx 


+ 


dT xy 

dy 


fcxz 

dz 


0 


(4.1) 


dx + dy 



0. 


Since the distribution of a y , and t,* are determined from the two-dimensional sol- 
ution to the problem, the unknowns are t„ and -zy, . They can be found from the in- 
tegration of eq. (4.1) with respect to z, i.e., 


l xz 


~J< 


do, 


xy 


dx 


dy 


) dz + f(x,y ) 


(4.2) 


yz ” 



-^)dz+g(x,y), 


where f(x,y) and g(x,y) are unknown functions of x and y, and the definitions of z and 
the layer interfaces are illustrated in Fig. 4-1. The key constitutive equation for clas- 
sical lamination theory is 




Ql2 ©16 
Q22 Q26 
Q26 ^66 


£° + ZK° 

L; + ZK° y 1. 
Vxy + Z K° xy 


(4.3) 


By using the constitutive equation and the strain-displacement relations in the above 
equation, performing integration with respect to z on a layer-by-layer basis, and ap- 
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plying boundary conditions on the bottom surface and continuity conditions at the 
layer interfaces, explicit expressions for interlaminar stresses can be obtained. The 
interlaminar shear stresses within the k-th layer are expressed as follows : 


l XZ 
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(4.4) 


(4.5) 


Note that for the purpose of the simplification of the equations, the notations for 
strains and curvatures are changed in above equations. The quantities £?, e8, eg, /c?, 
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kI and kJ correspond to e$, c°, yS*. k°, k°, and k%, respectively. As mentioned earlier, 
the traction-free boundary conditions on the bottom surface and the continuity con- 
ditions at the layer interfaces were used in the process of deriving the above 
equations. The traction-free boundary condition on the top surface remains redun- 
dant. However, the traction-free condition there is automatically satisfied. This can 
be shown as follows. The shear stress r„ on the top surface can be obtained by 
substituting z = z N+1 and k = N into the above equations. The result of sustitution 
appears as 


- A 


11 


K 

bx 


— I A 


16 


dy 


+ A 


12 


fry 

dx 


+ A 


dy, 


5k, 


8k: 


5k, 


16 


+ A 


26 


5iEy 

dy 


+ A, 


66 


■ ^ J 

dx 

■ + Bn 

A 

dx 

■f b 12 

dx ® 16 

dx 

8y% 

+ Big ' 

_ 

F Bo r ■ 

8k° 

+ Sfifi 

bK°y 


(4.6) 


by 


by 


by 


5AL dN, 
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Following the same procedure, it can be shown that 

/ 8N ru 8N U 


r yz ~ 


xy 


bx 


+ ' 


(4.7) 


where A,, and B /y are extension stiffness matrix and bending stiffness matrix, respec- 
tively, defined as 


JL JL 

Ay = [ 2 qU dz, Bij = f 2 z Q* dz. (4.8) 

J —h J —h 

2 2 

From the first two of eq. (3.12) it can be seen that the conditions r yI = 0 and r„ = 0on 
the top surface of the plate are satisfied. One indication of the lack of convergence 
of the interlaminar stress calculations using numerical solutions is the lack of satis- 
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faction of the redundant traction-free boundary condition at z = z N + 1 . This will be 
discussed more shortly. 


4.2 Comparison with Other Soiutions 


From eqs.(4.4) and (4.5) it can be seen that the expressions for interlaminar shear 
stresses involve first-order derivatives of strains and curvatures. Since the strains 
and curvatures are expressed as first-order derivatives of displacements and ro- 
tations in shear deformation theory, the interlaminar shear stresses involve the 
second-order derivatives of the displacements and rotations. If the solution for the 
displacements and rotations is expressed in a functional form which varies contin- 
uously in x and y, the second derivatives can be calculated exactly. However, the 
finite-element solution is a descrete solution. Therefore, in order to obtain the second 
derivatives of the displacements and rotations numerical derivatives have to be used. 
The second derivative obtained numerically can be inaccurate. Thus, in order to ex- 
amine the accuracy of interlaminar shear stress based on the finite-element method, 
simple cases are studied where the exact solutions for interlaminar stresses can be 
obtained in functional form. The numerical values from eq.(4-4) and eq. (4-5) based 
on finite-element calculations are compared with exact solutions which utilize these 
same equations but which can take advantage of closed-form expressions for the 
derivatives. The first case compared is the comparison in the geometrically linear 
range. The second case is a comparison in the geometrically nonlinear range of de- 
flections. 
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4.2.1 Simply Supported Cross-Ply Plate Under Uniform Loading. 


The governing equation within the framework of linear classical lamination theory 
for the out-of-plane deflection of a symmetric cross-ply laminate is 

D 1 ,-^ + 2(D, 2 + 2D 66 )-4tT + 0 22 -^--(7. (4.9) 

dx dx dy dy 

When the four edges of the plate are simply supported, the solution for w can be ex- 
pressed in the form 


M N 

w= X X lv - cos ( ir f L ) cos (-^r L ). ( 4 - 10 ) 

m = 1,3 n = 1,3 

where L and W are the plate length in the x direction and plate width in the y direc- 
tion, respectively, and the origin of the coordinate system is at the geometric center 
of the plate. The W mn are the unknown coefficients to be solved as a function of the 
applied load. From the governing equation and assumed solution form it can be ob- 
served that if the transverse load is expanded into a double cosine series, the sol- 
ution can be obtained easily. The transverse load can be expanded into a double 
cosine series as 


M N 

<7= X X Q m nCos(-^f-)cos(^-). (4.11) 

m = 1,3 n =1,3 

When q is a constant (uniform loading), the coefficients Q mn are obtained as 
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Substituting the expanded form for q and assumed form for w into the governing 
equation, eq (4.9), the expression for the unknown coefficients W mn can be obtained 
as 


W = 

wv mn 
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The solution developed here is exact as it satisfies all displacement and moment 
boundary conditions. The only issue is using a sufficient number of terms in the se- 
ries. 


The midplane curvatures are obtained by the following equations : 
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m = 1 ,3 n = 1,3 
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( 4 . 15 ) 
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- * £ i m r ) (T) s K nr w ^ )• (4 - ,6) 


m = 1,3 n= 1,3 


For this problem the inplane displacements, and thus the midplane strains, are zero. 
Substituting the above expressions into eqs. (4.4) and (4.5), interlaminar shear 
stresses are obtained. For a square plate with 15 terms in m and n (M — 15, N — 15), 
a sufficiently converged solution can be obtained. The stresses from this solution are 
compared with those calculated by the finite-element method. In the finite-element 
method the shape functions are used for the calculation of the curvatures and their 
derivatives. Figure 4-2 shows the comparison between the series and finite-element 
solutions at several locations within a [0/90] s plate. The variables appearing in Fig. 
4-2 and the ensuing figures are defined as 


x = 




(4.17) 


Sixteen elements for one-quarter of the plate (4 elements in the x direction and 4 el- 
ements in the y direction) were used for the finite-element calculation. The 
interlaminar shear stresses were compared at the centroid of each element. The 
agreement between the exact and finite-element solution is very good. 


4.2.2 Interlaminar Shear Stresses in Nonlinear Cylindrical Bending 

In the postbuckling problem geometrically nonlinear effects are important and it is 
not clear that conclusions obtained regarding the accuracy of the interlaminar 
stresses for the geometrically linear problem can be applied. As in the previous 
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1>ja 

(a) Txz along 7=0 0625 



Tyz 

(b) Xyz along x=0.0625 

Fig. 4*2 Comparison of interlaminar shear stresses between exact solution and FEM solution 
for a simply supported square [0/90]* plates under uniform loading. 
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sections for geometrically linear problems, a geometrically nonlinear problem for 
which there is a closed-form solution available, and for which interlaminar stresses 
can be computed, should be used for checking the finite-element calculations. In the 
course of this study, for the purpose of comparing interlaminar shear stresses, a ge- 
ometrically nonlinear problem was posed and solved in closed form. Interlaminar 
stresses calculated from this solution were compared with the stresses computed 
using the finite-element approach for the same problem. The problem is essentially 
the geometrically nonlinear analog of the cylindrical bending studied by Pagano [50]. 
To follow is the exact solution for the nonlinear cylindrical bending problem. 


The governing equations for classical plate bending theory including geometric 
nonlinearities are 


dx 
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xy 
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(4.20) 


where Af*. N y , and N*, are force resultants, M „ , M y , and M „ are moment resultants, 
w is the out-of-plane displacement of the plate's reference surface, and q is the 
transverse loading. The basic assumptions for cylindrical bending are 


v = 0, 



(4.21) 
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where v is the reference surface displacement in the y direction. With these as- 
sumptions the above governing equations are simplified to 


dN x 

dx 
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(4.22) 


d 2 M x 

dx 2 
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(4.23) 


Equation (4.19) is not involved in the solution procedure. The terms of importance in 
the constitutive equation for a symmetric laminate are expressed as 


N x ~ A 11 c x 



du 1 
dx + 2 



(4.24) 
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As usual, u is the displacement in the x direction of the plate's reference surface. 
Due to the basic assumptions for cylindrical bending, eq. (4.21), cj , y% , k? , and k% 
become zero. The boundary conditions studied are shown in Fig. 4-3. Note that the 
location x = 0 is at the plate's midspan. The use of symmetry is employed in the ap- 
plication of the boundary conditions and only one half of the plate is analyzed. 


The above equations are nondimensionalized as 



(4.26) 


By substituting eq. (4.26) into eqs. (4.24) and (4.25), nondimensionalized equations for 
N x and M x are obtained as 
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Fig. 4-3 Boundary conditions for nonlinear cylindrical bending. 
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(4.27) 


fT- du ( L v 1 , dw 2 

N * = W ( 7T ) + T ( tfx >■ 


— _ d 2 w 

M x 2 ’ 

dx 


where 


A/ x =-^-rA/ x , 

*11* 




M x = -pr~r M x . 
D^A 


(4.28) 


(4.29) 


(4.30) 


By substituting eq. (4.26), (4.29), and (4.30) into eq. (4.23), a nondimensionalized 
equation for eq. (4.23) is obtained as 


d 2 M x . A„H 2 _ 
dx 2 Du dx 2 


+ q = 0 , 


(4.31) 


where q is given by 


L 4 

q D 11 /? * 


(4.32) 


From eq. (4.22) it can be seen that /V„ therefore /V„, is a constant. This result that A/, 
is constant makes the problem quite simple. By substituting eq. (4.30) into eq. (4.31) 
a fourth-order ordinary equation for vv results, namely, 


dV 

dx 4 


2 d 2 w - 
a — r = d. 
dx 2 


(4.33) 


where 
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(4.34) 


and q can be a general function of x . For simplicity a uniform transverse loading 
case is considered. By substituting the solution for w into eq. (4.27), u can be de- 
termined. Considering that w is an even function of x , the solutions for u and vv are 
expressed as 


u = a^x — a 2 1 


sinh(2ax) * 

4a T 


+ a 3 < 


x cosh(ax) sinh(ax) 


}-*4 


(4.35) 


w = C 1 cosh(ax) — 



x 2 + C 2> 


(4.36) 


where 



C 2 a 2 /? 
a 2= 2 L 




(4.37) 


In these equations JT X , C,, and C* are unknown coefficients. The above solution satis- 
fies the symmetry conditions at the center of the plate. Applying the boundary con- 
ditions at the right end, the following three equations are derived : 


a ! f sinh(a) 

2 32 | 4a 4 



cosh( -|- ) 
2a 


sinh(-|-) 


24 


0 , 


(4.38) 


C, cosh( ) + C 2 = 0, 

* 8a 


(4.39) 


C^ 2 cosh( ■ = 0. 

* a 


(4.40) 
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In the above equations a is a function of N„. From eq. (4.37) a, , a 2 , a 3 , and a 4 are 
therefore functions of N x and C : . In the above equations an expression for C, in terms 
ofw; can be obtained from eq. (4.40). Substitution of this expression for C, into eq. 
(4.38) gives one nonlinear equation for N x as a function of the uniform load q. The 
final equation is solved numerically for N„ as a function of q , and C, and C 2 are de- 
termined by back sustitution into eqs. (4.39) and (4.40). 

With the displacement response of the plate known, the strains and curvatures, 
and their derivatives, can be used in eqs. (4.4) and (4.5) and the interlaminar stresses 
computed. The interlaminar shear stresses calculated by using the above exact sol- 
ution were compared with finite-element calculations. A [0/90] s cross-ply laminate 
with a load level of q = 1000 is studied. At this load level the maximum deflection at 
the center of the plate is vv = 1.80 , i.e., the plate is deflected 1.80 times its thickness. 
With the specific boundary conditions this deflection is in the geometrically nonlinear 
range. With 2, 5, and 10 elements representing one half of a strip across the width 
of the plate, the convergence of the finite-element solution is studied at the 
centroidal and other Gauss points in the element at the supported edge. This ele- 
ment is studied because the interlaminar shear stress is the highest at the supported 
edge. Figure 4-4 shows the comparison of the finite-element solution and the exact 
solution at the center Gauss point of the element. Even with 2 elements the shear 
stress on the top surface is zero. As was shown earlier, in the exact sense this has 
to be the case. However, with numerically obtained second derivatives there is no 
guarantee that this will happen. Despite the shear stress returning to zero at the top 
surface, with two elements there is a difference between the finite-element calcu- 
lation and the exact solution for shear stress within the thickness of the laminate. 
As the number of elements is increased the finite-element solution approaches the 
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(a) two elements mesh 
(r=0.375) 



(b) Five elements mesh 
(x=0.45) 



(iT-0.475) 

Fig. 4-4 Comparison of interlaminar shear stresses between exact solution and FEM solution 
in nonlinear range at center Gauss point In end element for a [0/90] s laminate in 
cylindrical bending, g = 1000. 
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exact solution. (Note that as the number of elements increases, the location x where 
the comparison is taking place moves closer to the edge of the plate, x = 0.50 , and 
the shear stress actually increases. ) 

Figure 4-5 shows the convergence study at a Gauss point other than the center 
Gauss point in the end element. Even with 10 elements the agreement is not good. 
From this it can be seen that the convergence of interlaminar stress at the Gauss 
point in the center of the element is better than at other Gauss points. Chaudhuri and 
Seide discussed this in ref. 51. 

Fig. 4-6 shows the comparison of the finite-element solution with the exact sol- 
ution for a [±45/0/90]* laminate with a load level of q = 1000. In this comparison 
Du and Da are set to zero. At the load level of q = 1000, w at the center of the strip 
was 1.58. With 10 elements representing the half-width of the plate, the agreement 
at the centroid of the element at the supported edge is good. 


4.3 Results for Plates with Holes 


With the above examples it has been shown that the interlaminar shear stress 
calculation scheme using the finite-element method can be accurate. Some 
investigators[34, 36, 37] have mentioned that the interlaminar shear stresses along a 
node line can be an important parameter for determining the failure of solid plates 
which deform into more than one half-wave in the loading direction during 
postbuckling. In this study specific interest centers on the interlaminar stresses in 
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(a) Two elements mesh 
(x^0.472) 



(5T=0.489) 



Fig. 4-5 Comparison of Interlaminar shear stresses between exact solution and FEM solution 
In nonlinear range at a Gauss point other thanthe center Gauss point in end element 
for a [0/90] s laminate in cylindrical bending, q = 1000. 
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0.5 



Fig. 4-6 Comparison of interlaminar shear stresses between exact solution and FEM solution 
in nonlinear range at a center Gauss point In end element for [ ± 45/0/90] 2S laminate 
in cylindrical bending, q = 1000. 
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plates with holes loaded in a postbuckled range. In this section for the plates with 
holes which have deformed into two half-waves in the loading direction, interlaminar 
shear stresses will be studied. First, convergence of interlaminar stress calculations 
as a function of the mesh refinement will be studied. Then, insight into the 
interlaminar shear stress distribution on the midplane will be discussed. Finally, 
through-the-thickness distributions of interlaminar shear stresses for each laminate 
in this study will be presented. Shown in Fig. 4-7 are several regions of interest 
within the plate, namely the net-section hole edge, half way between the hole edge 
and simple support along the net-section, the region just inside the simple supports 
on the net-section, and the portion of the plate extending beyond the simple supports 
on the net-section. 

For the convergence study of interlaminar shear stresses three levels of mesh 
refinement are used. They are a 57 element, a 121 element, and a 209 element mesh 
representing a quarter plate, as shown in Fig. 4-7. With a [±45/0/90] 2S plate 
interlaminar shear stresses on the elements along the net-section were studied for 
the above three levels of mesh refinements. The load level (P/EA = 0.00152) was se- 
lected so that the maximum out-of-plane deflection of the plate was about two plate 
thicknesses. As shown in Fig. 4-4 and Fig. 4-5, interlaminar shear stress calculations 
with the present method are most accurate at the center of the element. Thus, cal- 
culations are made at the centers of the elements. 

One location of interest is the net-section hole edge, the region marked by the 
double-crossed hatch in Fig. 4-7. The results from three mesh refinements for the 
through-the-thickness distribution of t„ at the center of the element at the net-section 
hole edge are shown in Fig. 4-8. In this figure t„ is normalized with respect to the 
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net-section extension 

net-section simple 
support 

net-section between 
hole edge and simple 
support 

HH net-section hole edge 



(a) 57 element mesh 



(b) 121 element mesh 



Fig. 4*7 Finite-element meshes and regions of interest used to study interlaminar stresses. 


Calculation of Interlaminar Shear Stresses 


94 




X Y7 

X1.0E-04 


Observation point 


No. of elements 

X 

y 

57 

0.176 

1.778 

121 

0.111 

1.688 

209 

0.077 

1.570 


Fig. 4-8 Dependence of the distribution of t„ through the thickness on mesh size for a 
[±45/0/90] 2S laminate at P/EA = 0.00152 at net-section hole edge (Two half-waves 
in the loading direction). 
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equivalent Young's modulus of the laminate in the x direction. This modulus is listed 
in Table 2-4. In this comparison, since the centers of the elements at the net-section 
hole edge for each mesh refinement do not coincide with each other, the stresses are 
not computed at exactly the same point. However, some aspects of the stress calcu- 
lations can be observed from this figure. As discussed with eqs. (4.6)-(4.8), the re- 
dundant traction free condition on the top surface in the derivation of eq. (4.4) and eq. 
(4.5) must be satisfied. Figure 4-8 shows that this condition is satisfied at the net- 
section hole edge with just 57 elements in the mesh. Also, it is clear that calcu- 
lated on the midplane with the 57 element mesh is quite different than t„ calculated 
with the 209 element mesh. As the mesh is refined, the observation point approaches 
the net-section hole edge. With a small movement toward the hole edge t„ on the 
midplane calculated with the 209 element mesh is about 4 times bigger than that 
calculated with the 57 element mesh. This shows that t„ at this region changes very 
rapidly. Using only a small number of elements could result in missing this rapid 
change of interlaminar shear stresses near the hole boundary. 

As a next observation location, the region between the hole edge and the simple 
support along the net-section was selected. The through-the-thickness distribution 
of r„ at this location is shown in Fig. 4-9 for each mesh refinement. This figure shows 
that as the number of elements is increased, r„ on the top surface approaches zero. 
Also, the figure shows that at this point t„ on the midplane is about one-third of the 
value at the net-section hole edge. 

One more location of interest is the extension region outside the simple support. 
The element at this region is marked by dots in Fig. 4-7 and is labeled as the net- 
section extension region. Through the thickness distribution of t xi at this location is 
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Observation point 


No. of elements 

X 

y 

57 

0.351 

3.536 

121 

0.214 

3.636 

209 

0.160 

3.540 


Fig. 4-9 Dependence of the distribution of t„ through the thickness on mesh size for a 
[ ± 45/0/90] 2S laminate at P/EA = 0.00152 between the net-section hole edge and the 
simple support (Two half-waves In the loading direction). 
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Observation point 


No. of elements 

X 

y 

57 

0.497 

5.25 

121 

0.329 

5.25 

209 

0.246 

5.25 


Fig. 4-10 Dependence of the distribution of x„ through the thickness on mesh size for a 
[±45/0/90]2 s laminate at P/EA = 0.001 52 at the net-section extension region outside 
the simple support (Two half-waves in the loading direction). 
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shown in Fig. 4-10. The value of r„ at this location is relatively small, and the sign is 
opposite to the sign at the other two locations. 

As a convergence test for T n , the center of the element just inside the simple 
support, a region marked by a single hatch in Fig. 4-7, was selected. The distribution 
of Tyi in the thickness direction at this location is shown in Fig. 4-11. The 57 element 
mesh shows poor convergence. As the number of elements is increased, the con- 
vergence is improved as indicated by the shear stress going to zero at the top sur- 
face. The value of at this location is also relatively small compared to r„ at the 
net-section hole edge. 

From this convergence study, it is shown that the convergence of interlaminar 
shear stresses is good at the locations where the interlaminar shear stresses are 
large. The convergence is not good at the location where the shear stress is small. 
However, the locations of interest for failure study are the areas where the stresses 
are indeed large. Thus, with the above results, it can be expected that by using a 209 
element mesh an accurate and realistic study of failure can be achieved. 

In order to learn more about the interlaminar shear stresses within the plate, the 
interlaminar shear stress on the geometric midplane is computed as a function of 
location along several loci. In particular, the variation of t„ and r ^ along the centers 
of the elements at the net-section and along the centers of the elements on the hole 
boundary is examined. In Fig. 4-12 the variation of interlaminar shear stresses along 
the centroids of the elements at the net-section is shown. The y coordinate of the 
center of the elements is normalized by the plate width W, and the interlaminar shear 
stresses are again normalized by the equivalent Young's modulus of a [ ± 45/0/90] 2S 
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Fig. 4*11 Dependence of the distribution of t n through the thickness on mesh size for a 
[±45/0/90 ]j 4 laminate at P/EA = 0.00152 at the net-section just inside the simple 
support (Two half-waves in the loading direction). 
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Fig. 4-12 Distribution of interlaminar shear stresses on the midplane of [±45/0/90] 2S along 
the center of the elements at net-section. 
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plate in the x direction. In Fig. 4-12 it is clearly shown that r„ increases very rapidly 
as the hole boundary is approached, as discussed at the beginning of this section. 
On the other hand, r n is negligibly small. As expected, v is largest at the simple 
support. In Fig. 4-13 the variation of interlaminar shear stresses along the centers 
of the elements on the hole boundary is shown. The stresses are shown as a function 
of 0, where 0 = 90° is the net-section location. It can be seen that up to 0 = 45° the 
value of t„ is small, but for 0 > 45° t„ increases very rapidly and reaches its highest 
value at net section hole edge. Along the hole boundary r n has relatively small val- 
ues compared to t « at the net-section hole edge. The distribution of jyz around the 
hole boundary is roughly symmetric with respect to 0 = 45° . 

In Fig. 4-14 the through-the-thickness distribution of interlaminar shear stress 
r„ for the different laminates in this study is shown. The stress was calculated at the 
net-section hole edge. The 209 element mesh was used in this calculation. Note that 
the equivalent Young's moduli used to normalize the stresses are different from plate 
to plate. Also, the load levels are different from plate to plate. However, the load level 
was selected so that it represents the postbuckling range where the maximum out- 
of-deflection is approximately two plate thicknesses. From Fig. 4-14 it can be ob- 
served that t „ changes rapidly in the 0° and ±45° layers, and less rapidly in the 90° 
layers. Each figure shows symmetry with respect to the midplane. In the [ ± 45/0 6 ]s 
laminate the inner 12 layers have the same fiber orientation. In Fig. 4-14 (c) it can be 
seen that the distribution of r„ in these grouped 0° layers results in almost a 
parabolic distribution, as can be observed in an isotropic material. In Fig. 4-14 (d) the 
distribution of in the [ ± 45]« laminate does exhibit a parabolic shape. 
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X 1.0E-04 


(c) [ ±45/0 els plate at P/EA=0.001 17. 


-xz 


X 1.0E-04 


(d) [ +45] 4S plate at P/EA=0.00240. 


Fig. 4-14 Through the thickness distribution of i„ tor each plate at the net-section hole edge. 
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The through-the-thickness distribution of x n was observed at the center of the 
element at the net-section simple support for the laminates in this study. Their dis- 
tribution with thickness is shown in Fig. 4-15. In general on the top surface is not 
zero. This means that x n is not fully converged. However, the value of x n at this lo- 
cation is an order of magnitude smaller than t„ at the net-section hole edge. It can 
also be observed that the distribution of r „ is not symmetric about the midplane. In 
keeping with the results of Fig. 4-14, it is seen that for the 0° layers t,* does not vary 
significantly with z. In terms of actual value of interlaminar shear stress, the 
[ ± 45/0*]« and [ ± 45/0 6 ]s laminates experience the largest values of r„ and the 
[±45/0/90]a and [±45]« laminates experience about 25% less than these two 
laminates. As regards interlaminar stress x „ , the [ + 45/0/90] 2S , [ ± 45/0 2 ] 2S , and 
[±45/0 6 ] laminates experience the same level of x n and the [±45]« laminate ex- 
periences less. 

With the stresses being computed and the calculations being validated, attention 
turns to the prediction of failure using the stress calculation. This is the subject of the 
next chapter. 
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5.0 Numerical Failure Analysis 


As has been frequently stated, the primary goal of this study is to predict failure. 
The load level, the location, and the mode of the failure are of interest. The first step 
in the analysis of failure requires an accurate determination of the stresses. The cal- 
culations for the stresses have been discussed in the previous chapter. The second 
step in a failure analysis is the selection of a failure criterion and the development 
of failure scenarios using the stresses and the failure criterion. Here the maximum 
stress failure criterion is used. This criterion is selected because it is easy to under- 
stand and easy to use. Complicating data regarding failure due to stress interaction, 
as is sometimes required in more complicated failure theories, such as the Tsai-Wu 
criterion, are not needed. 

Because of the numerical nature of this study the failure analysis in this study is 
carried out layer-by-layer by examining the stresses at the centroid of each element. 
Recall from the last chapter that the interlaminar stress calculations at the centroid 
are more accurate than at other Gauss points. With the stresses in each element and 
each layer computed, the following failure parameters are defined and calculated : 
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(5.1) 


Oy Oy 

Ry = —r , when Oy ^ 0, or Ry = -j— , when a 1 < 0, 
x t x c 


On ^2 

R 2 = — r , when ct 2 > 0, or R 2 = -rr~ , when ct 2 < 0, 


Rio — 


r 12 I 


12- S 


Rl3- 5 


12 
t 13 I 


13 


*23 - S 


t 23 I 


23 


(5.2) 

(5.3) 

(5.4) 

(5.5) 


In eqs. (5.1)-(5.5) < 7 , is the normal stress in the fiber direction, ct 2 is the normal stress 
transverse to the fiber direction in the plane of the layer, i.e., the 1-2 plane, ti 2 is the 
shear stress in the 1-2 plane, ti 3 is the interlaminar shear stress in the 1-3 plane, and 
t 23 is the interlaminar shear stress in the 2-3 plane. In addition, X, is the tensile 
strength in the fiber direction, X c is the compression strength in the fiber direction, 
Y t is the tension strength transverse to the fiber direction in the plane of the layer, 
Y e is the compression strength transverse to the fiber direction in the plane of the 
layer, and S 12 is the inplane shear strength. These strengths are defined in the mate- 
rial coordinate system, i.e., the 1-2-3 coordinate system. Also, S 13 is the interlaminar 
shear strength in the 1-3 plane, and S 23 is the interlaminar shear strength in the 2-3 
plane. The values of these strength parameters for AS4/3502 used in this study are 
listed in Table 5-1. These values are taken from reference [53, 54]. Note that S 12 and 
S 13 are assumed equal. This is a reasonable assumption. From the numerical analy- 
sis formulated in earlier chapters the stresses are calculated in x-y-z coordinates. 
These stresses are then transformed into the stresses in the material coordinate 
system for failure analysis with the maximum stress criterion. The stress transfor- 
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Table 5*1 Strengths of AS4/3502 used in this study. 


Strength parameter 

Strength 

(Ksi) 

tensile strength in fiber direction ( Xt) 

200 

compressive strength in fiber direction ( X c ) 

-200 

tensile strength 

transverse to fiber direction ( Y t ) 

7 

compressive strength 
transverse to fiber direction ( Y c ) 

-35 

shear strength in 1 -2 or 1 -3 plane ( So, So ) 

9.4 

shear strength in 2-3 plane ( S23) 

8 
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mation rule can be found in any book dealing with mechanics of composite materials, 
as an example, in reference [55, 56] 

The ratios R, , R 2 , R u , , and R 23 are used as indicators of failure. When Ri is 

equal to 1.0, it is judged that fiber failure occurs due to normal stress in the fiber di- 
rection. In this study it is called fiber failure. In the same fashion, when R z and R n 
are equal to 1.0, it is called transverse failure and inplane shear failure, respectively. 
There are three stress components that can contribute to fiber failure. They are a , , 
a n , and ct 13 . Even though fiber failure can be enhanced by the stresses a n and ct 13 , fi- 
ber failure in this study means the failure of a layer in the fiber direction due to <7,. 

To follow is a considerable amount of information regarding predictions of failure 
in the plates. As mentioned in chapter 3, a plate loaded in the postbuckling range 
can experience unexpected mode shape changes. For the laminates in this study 
failure loads, failure modes, and failure locations are predicted for both the one half- 
wave configuration and the two half-wave configuration. For plates with holes the 
failure analysis results are presented in Table 5-2 for the one half-wave configuration 
and Table 5-3 for the two half-wave configuration. For plates without holes, i.e., solid 
plates, failure analysis results are presented in Table 5-4 for the one half-wave con- 
figuration and Table 5-5 for the two half-wave configuration. The tables are based on 
a first-ply failure analysis. In these tables the element number and the layer number 
of the failed layer are given, where layer number is counted from the bottom concave 
surface of the plate. The element numbering scheme can be found in Fig. 5-1 for the 
plates with holes and Fig. 5-2 for the solid plates. Since the initial deflection due to 
the imperfection is assumed to be in the positive z direction, the plate deflects furthur 
in the positive z direction when loaded. Thus here the bottom surface is on the op- 
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Table 5-2 Failure loads, failure modes, and failure locations for plates with holes for the one 
half-wave configuration. 


Laminate 

Failure 

mode 

N x (lb/in) 

P 

EA 

Failure location * 

[ ± 45/0/90] 2S 

R i 

1605 

0.00256 

element 197 
3-rd layer ( 0°) 

r 2 

1239 

0.00197 

element 189 
16-th layer (45°) 

R 12 

1140 

0.00182 

element 197 
1 -st layer (45° ) 

[ ± 45/0 2 ]2s 

R 1 

1989 

0.00217 

element 169 
2 -nd layer (- 45 °) 

R 2 

1231 

0.00134 

element 205 
14-th layer ( q°) 

R 12 

1342 

0.00146 

element 197 
1 -st layer ( 45 ° ) 

[ ± 45/0 6 ]2s 

R 1 

2030 

0.00163 

element 169 
2 -nd layer (- 45 °) 

R 2 

1247 

0.00100 

element 205 
14-th layer ( q°) 

R 12 

1613 

0.00129 

element 197 
1 -st layer (45°) 

[ ± 45^43 

R 1 

1447 

0.00592 

element 129 
2 -nd layer ( 45 0 ) 

R 2 

1132 

0.00463 

element 140 
1 6-th layer ( 45 ° ) 

R 12 

461 

0.00189 

element 197 
1 -st layer ( 45 0 ) 


* Layer no. is counted from concave side. 
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Table 5-3 Failure loads, failure modes, and failure locations for plates with holes for the two 
half-wave configuration. 


Laminate 

Failure 

mode 

Nx (Ib/in) 

P 

EA 

Failure location * 


R i 

1532 

0.00244 

element 208 
1 -st layer ( 45 0 ) 

j 

[ ± 45/0/90] 2S 

r 2 

1168 

0.00186 

element 79 
16-th layer ( 45 °) 


Rl2 

1536 

0.00245 

element 170 
16-th layer ( 45 °) 


R 1 

1892 

0.00206 

element 208 
1 -st layer ( 45 °) 

[ ± 45/0 2 ] 2S 

R 2 

1232 

0.00134 

element 1 54 
14-th layer ( o°) 


R 12 

1681 

0.00183 

element 170 
16-th layer (45°) 


R 1 

1953 

0.00156 

element 208 
1-st layer ( 45 0 ) 

[ ± 45/0 e ] s 

R 2 

1318 

0.00106 

element 154 
1 4-th layer ( q° ) 


m 

1894 

0.00152 

element 208 
3-rd layer ( o° ) 


R i 

1260 

0.00515 

element 207 
1-st layer ( 45 0 ) 

[ ± 45] 4S 

R 2 

716 

0.00293 

element 79 
16-th layer(45°) 


R 12 

631 

0.00258 

element 158 
1 -st layer ( 45 °) 


* Layer no. is counted from concave side. 
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Table 5-4 Failure loads, failure modes, and failure locations for solid plates for the one 
half-wave configuration. 


Laminate 

Failure 

mode 

N x (lb/in) 

P 

EA 

Failure location 


R 1 

1758 

0.00280 

element 69 
2 -nd layer (. 45 °) 

[ ± 45/0/90] 2S 

r 2 

1544 

0.00246 

element 50 
1 6-th layer ( 45 0 ) 


Rl2 

1710 

0.00273 

element 62 
1 6-th layer ( 45 0 ) 


R 1 

2060 

0.00225 

element 69 
2-nd layer (-45°) 

[ ± 45/0 2 ] 2S 

r 2 

1468 

0.0016 

element 48 
1 4-th layer ( q° ) 


Rl2 

1859 

0.00203 

element 73 
1 -st layer ( 45 0 ) 


R 1 

2076 

0.00166 

element 68 
2 -nd layer (- 45 °) 

[±45/0 6 ] s 

R 2 

1425 

0.00114 

element 46 
14-th layer ( 0 °) 


R 12 

1919 

0.00154 

element 68 
3-rd layer ( 0°) 


R 1 

1443 

0.00590 

element 60 
2 -nd layer (- 45 0 ) 

[ ± 45]„ 

R 2 

1080 

0.00442 

element 50 
1 6-th layer ( 45 ° ) 


R 12 

738 

0.00302 

element 73 
1 -st layer ( 45 ° ) 


* Layer no. is counted from concave side. 
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Table 5-5 Failure loads, failure modes, and failure locations for solid plates for the two 
half-wave configuration. 


Laminate 

Failure 

mode 

N"x (Ib/in) 

P 

EA 

Failure location 


fl 

1606 

0.00256 

element 64 
1 -st layer ( 45 °) 

[ ± 45/0/90] 2S 

r 2 

1546 

0.00246 

element 48 
15-th layer (- 45 °) 


□ 

1568 

0.00250 

element 76 
16-th layer ( 45 °) 


R 1 

1990 

0.00217 

element 64 
1-st layer (45°) 

[±45/0 2 ] 2S 

r 2 

1340 

0.00146 

element 58 
14-th layer ( o°) 


R 12 

1748 

0.00191 

element 76 
1 6-th layer ( 45 ° ) 


R 1 

2096 

0.00168 

element 64 
1 -st layer ( 45 0 ) 

[ ± 45/0 6 ] 2S 

r 2 

1426 

0.00114 

element 58 
14-th layer ( o°) 


R 12 

1990 

0.00159 

element 67 
1 6-th layer ( 45 ° ) 


R 1 

1347 

0.00551 

element 64 
1 -st layer ( 45 ° ) 

[ ± 45] 4S 

R 2 

1117 

0.00457 

element 39 
15-th layer (- 45 °) 


m 

745 

0.00305 

element 66 
1 -st layer (45° ) 


* Layer no. is counted from concave side. 
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Fig. 5-1 Element numbering for a square plate with hole with a 209 element mesh. 
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posite side to the direction of the plate deflection, i.e., on the concave side. These 
tables mean that if the plates remain in the one half-wave configuration without a 
mode shape change, failure occurs in the manner predicted in Table 5-2 for plates 
with holes and Table 5-4 for solid plates. If the plates have a mode shape change, 
failure occurs in the manner predicted in Table 5-3 for plates with holes and Table 5-5 
for solid plates. In Tables 5-2 through 5-5 failure loads are presented in terms of av- 
erage applied end force resultant as defined in eq. (2.21). The nondimensionalized 
values of these failure loads in the form of P/EA are also shown in these tables. As 
an example of the deflected shapes, the three-dimensional surface drawings of the 
deformed shape of the quasi-isotropic plate are presented for plates with and without 
holes for both the one half-wave and two half-wave configurations. 

In figures to follow the predicted failure loads for each plate are marked on the 
load vs. end shortening relations for each deformation configuration. Following the 
load vs. end shortening relations for each plate, figures which show the failure lo- 
cations, failure modes at that location, and the number of failed layers at that location 
are presented. These figures go beyond a first-ply failure analysis and provide in- 
formation on multiple-ply failures, or interpreted differently, damage accumulation. 
These multiple-ply failures are assumed to occur after the first-ply failure and the 
calculations are done for loads greater than the first-ply failure loads. The analysis for 
the multiple-ply failure is not altogether correct, however. If first failure occurs in a 
certain mode and at a certain location in the plate, after this first-ply failure plate 
stiffnesses are no longer the same as before the first-ply failure. The stiffness change 
due to the first failure must be taken into account in the analysis for further ply fail- 
ures. This makes the analysis much more complicated. In the damage accumulation 
scenario here, the stiffness changes due to the first failure are neglected. Except for 
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the [ + 45]« plates, the load level used for the damage accumulation scenarios is 
slightly above the first fiber failure load. For the [ ± 45]4s plates the fiber failure load 
is quite high compared to other failure loads. It is obvious that the final failure of the 
[±45]« plates will occur at much lower load than the fiber failure load. Therefore, 
a load level slightly above the first transverse failure load is selected to show the 
damage scenarios of these plates. 

In chapter 7 the applicability of the first-ply failure theory will be examined by 
comparing the numerical results from the first-ply failure theory with the experimental 
results. Extensive comparisons for load-end shortening relations, load-deflection re- 
lations, strains, failure locations, failure loads, and failure modes will be presented in 
chapter 7. It will be noted in advance that contrary to the previous investigators' ex- 
perimental results for rectangular plates [34, 36, 37], interlaminar shear failure is not 
predicted in this failure analysis for square plates. This will be discussed more in 
chapter 7. 


5.1 Failure of Plates with Holes 

5.1.1 The [±45/0/90] 2S Plate with a Hole 

Figure 5-3 shows the deflected shape of the [ ± 45/0/90]^ plate with a hole in the 
one half-wave configuration. The load level is P/EA = 0.00258. Referring to Table 2-3 
it is seen that this level is about 6 times the first buckling load. The lines where the 
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Si 


line of simple support 


line of clamped support 



Fig. 5-3 Shape of the [ ± 45/0/90],* plate with a hole deflected in the one half-wave config 
uration at P/EA = 0,00258. 
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supports are assumed to act are shown. It can be observed that a curvature opposite 
to the curvature at the center of the plate occurs near the corner of the clamped 
support and simple support. The deflected shape of a quarter-plate in the two half- 
wave configuration is shown in Fig. 5-4. The load level used is about 4 times the 
buckling load corresponding to the two half-wave configuration in the loading direc- 
tion. Figure 5-5 shows the predicted failures superimposed on the load vs. end 
shortening relations for the [ ± 45/0/903es plate with a hole. For the one half-wave 
configuration, the solid line, as the load is increased, inplane shear failure is pre- 
dicted to occur first and it is predicted to occur at the net-section hole edge, element 
197, in the bottom layer. Continuing to load in the one half-wave configuration, 
transverse failure in the top layer of element 189 follows the inplane shear failure. 
Fiber failure finally occurs at the net-section hole edge, in element 197, at a much 
higher load than the loads of the other two failure modes. Table 5-2 summarizes 
these findings and it records that first inplane shear failure and first fiber failure occur 
at the net-section hole edge. From this it can be deduced that the net-section hole 
edge is under severe stress. Figure 5-6, which is one of the aformentioned damage 
accumulation scenarios, shows clearly that the net-section hole edge is severely 
damaged in the postbuckling loading of quasi-isotropic plate responding in the one 
half-wave configuration. The load level used in the scenario is P/EA = 0.00258, 
slightly above the fiber failure load of P/EA = 0.00256. The scenario shows that all 
three failure modes appear in that area, and it shows that a number of layers at that 
location fail. Figure 5-6 also shows that at the hole edge on the x axis 5 or 6 layers 
are damaged by transverse failure. This is because of the high value of N y at this lo- 
cation. Up to the load level of first fiber failure no interlaminar shear failure is pre- 
dicted. At the load level of P/EA = 0.00258 the highest value of R :3 is 0.3, and the 
highest value of is 0.34. 
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Fig. 5-4 Shape of the [ ± 45/0/90]j S plate with a hole deflected in the two half-wave conflg 
u ration at P/EA = 0.00255. 
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Fig. 5-5 Failures and load-end shortening relations for [ ± 45/0/90]2 S plate with a hole. 
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• Fiber failure 
" Trans, failure 
a Shear failure 



Fig. 5-6 Predicted failure locations and failure modes of the [ ± 45/0/90]2s plate with a hole 
in the one half-wave configuration. 
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Referring to Fig. 5-5 and Table 5-3 for the two half-wave configuration, the dashed 
line, as the load is increased, transverse failure is predicted to occur first and it is 
predicted to occur at the hole edge about half way between the net-section and the 
x axis, in the element 79. Inplane shear failure and fiber failure are predicted to occur 
at similar load levels, but at a level much higher than the transverse failure load. The 
first inplane shear failure is predicted to occur in the element 170, in the extension 
area away from the net-section. The first fiber failure occurs at element 208, just in- 
side the simple support along the net-section. This is a significantly different location 
than for the one half-wave case. However, the fiber failure load of the two half-wave 
configuration is only slightly lower than the fiber failure load of the one half-wave 
configuration. 

Considering the damage accumulation scenario of Fig. 5-7, at a load level slightly 
greater than first fiber or first inplane shear failure, transverse failure appears rather 
widespread over an area on the hole edge, near where the first transverse failure was 
predicted (element 79), and also over an area away from the net-section and simple 
support, in the central region of the quarter-plate. Inplane shear failure occurs out- 
side the simple support, slightly away from the net-section. Interlaminar shear failure 
is not predicted in the two half-wave configuration at loads as high as the first fiber 
failure load. At P/EA = 0.00255 the highest value of the failure parameter Ru is 0.44, 
and the highest value of is 0.54. These values are higher than those for the one 
half-wave configuration, but far from failure values. As mentioned at the beginning of 
this chapter, inplane shear failure can contribute to fiber failure. Therefore, inplane 
shear failure in addition to fiber failure can be a major factor of final failure of the 
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• Fiber failure 
“ Trans, failure 
a Shear failure 



Load: —=0.00255 
EA 


Fig. 5-7 Predicted failure locations and failure modes of the [±45/0/90] 2S plate with a hole 
in the two half-wave configuration. 
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plate. It is expected that the final failure of the plate may even occur outside the 
simple support, a location where severe inplane shear stresses, though not the 
maximum values, are predicted to occur. 


5.1.2 The [±45/0 2 ]2s Plate with a Hole 

Figure 5-8 shows the load vs. end shortening relation for the [ ± 45/0 2 ] 2S plate 
with a hole and with the failure loads and modes indicated. As can be seen, for both 
the one half-wave and two half-wave configuration, transverse failure occurs first as 
the plate is loaded from zero load. Figures 5-9 and 5-10 show the damage accumu- 
lation scenarios for the one and two half-wave configurations, respectively. 

Referring to Table 5-2, for the one half-wave configuration transverse failure oc- 
curs first and it occurs in the middle of the net-section between the net-section hole 
edge and simple support in element 205. Shortly after the first transverse failure the 
first inplane shear failure follows at the net-section hole edge in element 197. Fiber 
failure occurs at a load higher than the other two failure loads just inside the simple 
support, slightly away from the net-section. The fiber failure load for this plate in the 
one half-wave configuration is higher than the fiber failure load of the quasi-isotropic 
plate in the one half-wave configuration. 

Considering the damage accumulation scenario at a load of P/EA = 0.00222, just 
above the first fiber failure load of P/EA = 0.00217, as shown in Fig. 5-9 the transverse 
failure develops along the net-section line and also on the hole edge near the x axis. 
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P/EA 



Fig. 5-8 Failures and load-end shortening relations for the [±45/02] 2 s plate with a hole. 
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• Fiber failure 
“ Trans, failure 

▲ Shear failure ix 



12 3 4 

(Number in the figure represents the no. of failed layers.) 


Load : — — = 0.00222 
EA 


Fig. 5-9 Predicted failure locations and failure modes of the [ ± 45/02] 2S plate with a hole In 
the one half-wave configuration. 
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• Fiber failure 
- Trans, failure 

▲ Shear failure 



Fig. 5-10 Predicted failure locations and failure modes of the [ ± 45/Oj ]j S plate with a hole in 
the two half-wave configuration. 
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It is predicted that half of the laminate, i.e., 8 layers, is damaged on the hole edge 
near the x axis at this load just above the first fiber failure load. Inplane shear failure 
develops on the net-section hole edge and also along the simple support. The ex- 
tension area outside the simple support on the net-section, and the area just inside 
the simple support on the net-section are damaged by inplane shear stress. Fiber 
failure occurs just inside the simple support away from the net-section. Compared 
with the quasi-isotropic plate deformed in the one half-wave configuration, the net- 
section hole edge of this [±45/02] a laminate is under a less severe stress state. 
Again, interlaminar shear stresses are low. At slightly above the first fiber failure load 
of P/EA = 0.00222 the highest value of Ri 3 is 0.35, and the highest value of R i} is 0.40. 
It is expected that the final failure will occur along the net-section. 

For the two half-wave configuration, as the load is increased failure occurs in the 
sequence of transverse failure, then inplane shear failure , and then fiber failure, as 
in the one half-wave configuration. However, from Table 5-3 it can be seen that the 
locations of these failure modes are quite different than the locations of failure for the 
one half-wave configuration, and the load level for the transverse failure is different. 
The first transverse failure occurs between the simple support and the hole edge, 
away from the net-section, in element 154. The first inplane shear failure occurs in 
element 170, the extension area away from the net-section. First fiber failure occurs 
just inside the simple support on the net-section, element 208. The location of the 
first fiber failure is the same location of first fiber failure of the [ ± 45/0/90]^ plate in 
the two half-wave configuration. As seen from Fig. 5-8, transverse failure occurs at 
about the same load level as for the one half-wave configuration. In the one half- 
wave configuration, inplane shear failure occurs shortly after transverse failure. 
However, in the two half-wave configuration the inplane shear failure load is much 
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higher than the transverse failure load. The fiber failure load for the [ ± 45/0 2 ]« plate 
in the two half-wave configuration is slightly lower than the fiber failure load of the 
one half-wave configuration, but it is considerably higher than the failure load of the 
quasi-isotropic plate in the two half-wave configuration. Referring to the damage ac- 
cumulation scenario for the two half-wave configuration at the first fiber failure load, 
Fig. 5-10, like the quasi-isopropic plate in the two half-wave configuration, Fig. 5-7, 
transverse failure is widespread in the region between the simple support and hole 
edge, but away from the net-section. Severe inplane shear failure occurs near the 
simple support but slightly away from the net-section. No interlaminar shear failure 
is predicted at the first fiber failure load. At the fiber failure load, P/EA = 0.00206, the 
highest value of R 13 is 0.58, and the highest value of R 23 is 0.44. 

From Fig. 5-10, it is expected that the final failure of the plate will occur at the 
simple support at or slightly away from the net-section. Considering that the inplane 
shear failure load and fiber failure load are higher than those of the quasi-isotropic 
plate in the two half-wave configuration, the final failure load of this plate in the two 
half-wave configuration is expected to be greater than the final failure load of the 
quasi-isotropic plate in the two half-wave configuration. 


5.1.3 The [±45/0 6 ] s Plate with a Hole 

The load vs. end shortening relation for the [±45/0 e ] s plate with hole is shown 
in Fig. 5-11. For the one half-wave configuration the failure sequence is transverse 
failure, followed by inplane shear failure, followed by fiber failure. The inplane shear 
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Fig. 5-11 Failures and load-end shortening relations for the [ ± 45/0 plate with a hole. 
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failure load and fiber failure load are both higher than for the [ ± 45/0 2 ] 2S plate, which, 
as has been seen, are higher than these failure loads for the quasi-isotropic plate. 
Table 5-2 shows that the first transverse failure occurs at the center of the net- 
section, element 205. First inplane shear failure occurs at the net-section hole edge, 
element 197, and first fiber failure occurs at the same location as the first fiber failure 
of the [ + 45/0 2 ]« plate in the one half-wave configuration, element 169, along the 
simple support but away from the net-section. 

Figure 5-12 illustrates the damage accumulation scenario for the [ ± 45/0 8 ] s plate 
in the one half-wave configuration. Like the [±45/0 2 ] 2S plate in the one half-wave 
configuration, Fig. 5-9, transverse failure develops along the net-section line and at 
the hole edge on the x axis. However, by comparing Fig. 5-12, Fig. 5-9, and Fig. 5-6 
it can be seen that as the number of 0° layers are increased, the location of trans- 
verse failure on the net-section moves toward the simple support. Also, since a 
larger number of layers have failed, it can be seen that the transverse failure at the 
hole edge on the x axis is more severe for the [±45/0J s laminate than for the 
[±45/0 2 ] 2S laminate and for the [ ± 45/0/90]«. For the [±45/0 e ] s laminate all 12 0° 
layers on the hole edge on the x axis have failed in the transverse failure mode. For 
the [ + 45/0 2 ] 2S laminate all 8 0° layers on the hole edge at the x axis have failed in 
the transverse failure mode. For the [ ± 45/0 8 ]s plate at the fiber failure load, inplane 
shear failure develops on the net-section hole edge and net-section simple support. 
The net-section hole edge is under concentrated inplane shear failure. However, the 
inplane shear failure on the net-section simple support is spread out. The 
interlaminar shear failure parameters are far lower than the inplane failure parame- 
ters. The highest value of R» is 0.35, and the highest value of is 0.31, both com- 
puted at P/EA = 0.00167. 
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• Fiber failure 
- Trans, failure 
▲ Shear failure 



(Number in the figure represents the no. of failed layers.) 


Load : 


P 

EA 


0.00169 


Fig. 5*12 Predicted failure locations and failure inodes of the [±45/0b]s plate with a hole In 
the one half-wave configuration. 
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From Fig. 5-1 1 and Table 5-3 it can be seen that as the plate is loaded, for the two 
half-wave configuration transverse failure begins to develop at a slightly higher load 
than the transverse failure load of the one half-wave configuration. The location is 
element 154, the same location as the first transverse failure location of the 
[±45/0 2 ] 2S plate in the two half-wave configuration. The inplane shear failure load is 
higher than the transverse failure load, and the inplane shear failure occurs just in- 
side the net-section simple support, element 208. Shortly after the inplane shear 
failure occurs, fiber failure occurs at just inside the simple support on the net-section, 
also at element 208. Thus at element 208 both inplane shear failure and fiber failure 
occur. This is the same location as the fiber failure location of the [ ± 45/0/90] 2S and 
[±45/0 2 ] 2S plates with a hole in the two half-wave configuration. Compared to the 
[±45/0/90] 2S and [±45/0 2 ] 2S laminates, the fiber failure and inplane shear failure 
loads of the two half-wave configuration are higher for the [±45/0 6 ]s laminate. Con- 
sidering the damage accumulation scenario, Fig. 5-13, it can be seen that like the 
[ ± 45/0/90]zs and [ ± 45/0 2 ] 2S plates, Figs. 5-7 and 5-10, respectively, transverse fail- 
ure occurs in the central region of the quarter plate, away from the net-section. The 
area just inside the simple support is damaged by transverse failure and inplane 
shear failure. This is similar to the previous two plates. The highest value of failure 
parameter Rn is 0.61, and the highest value of ft 2 j is 0.28 at the load level of 
P/EA = 0.00161, a level slightly greater than the first fiber failure load. 
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• Fiber failure 
“ Trans, failure 

▲ Shear failure 



x 


y 


(Number in the figure represents the no. of failed layers.) 

Load: -^—=0.00161 
EA 


Fig. 5-13 Predicted failure locations and failure modes of the [ ± 45/0e]s plate with a hole In 
the two half-wave configuration. 
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5.1.4 The [±45] 4S Plate with a Hole 


The load vs. end shortening relation for the [±45]4s plate with hole is shown in 
Fig. 5-14. In the one half-wave configuration inplane shear failure occurs at a very 
low load level. Table 5-2 shows that the first inplane shear failure occurs at the net- 
section hole edge, element 197. For the other plates failure begins about at least 3.5 
times the buckling load. For this plate, however, failure begins at less than 2 times 
the buckling load. For this plate the transverse failure load is considerably greater 
than the inplane shear failure load. This has not been the case with the other plates. 
The first transverse failure occurs in the central region of the quarter plate, away from 
the net-section, in element 140. First fiber failure occurs in element 129, near the 
clamped end and near the simple support. This is the first plate to experience fiber 
failure in this location. 

The damage accumulation scenario, Fig. 5-15, shows that at a load slightly above 
the inplane shear failure load the plate is damaged severely by inplane shear failure. 
In the area outside the simple support and near the net-section all 16 layers have 
failed. Table 5-2 showed that the first inplane shear failure occurred at the net- 
section hole edge. But in Fig. 5-15 failure marks for inplane shear do not appear 
there. This means that even though the first inplane shear failure occurred at the 
net-section hole edge, the location of severe inplane shear moved to the simply 
supported area near the net-section. That is, even though the value of R\ 2 at the net- 
section hole edge is greater than 1.0 for a single layer, the value at the net-section 
simple support is greater than 1.0 for many more layers, i.e., 16 layers. The highest 
values of the interlaminar shear failure parameters are Rn = 0.19 and R zi = 0.22 at the 
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Fig. 5-14 Failures and load-end shortening relations for the [±45]4s plate with a hole. 
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- Trans, failure 
▲ Shear failure 


x 



y (Number in the figure represents the no. of failed layers.) 

Load : = 0.00485 

EA 


Fig. 5-15 Predicted failure locations and failure modes of the [ ± 453m plate with a hole in the 
one half-wave configuration. 
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load level slightly above the first inplane shear failure load. It is expected that failure 
begins outside the simple support near the net-section and propagates into the plate 
along the net-section. 

Referring to Fig. 5-14 and Table 5-3, as the plate is loaded from zero load, inplane 
shear failure appears first in the two half-wave configuration for this case also. The 
first inplane shear failure is predicted to occur at the net-section hole edge, element 
158. Compared to the other three plates, that is a slightly different location. The in- 
plane shear failure load in the two half-wave configuration is higher than the inplane 
shear failure load in the one half-wave configuration. Transverse failure begins at a 
considerably lower load for the two half-wave configuration than for the one half-wave 
configuration. It occurs at the hole edge half way between the net-section hole edge 
and x axis, element 79. Fiber failure occurs near the net-section simple support in 
element 207. Figure 5-16 shows the damage accumulation at a load level between the 
transverse failure load and the fiber failure load, i.e., at P/EA = 0.00394. It shows se- 
vere inplane shear failure along the simple support slightly away from the net- 
section. Sixteen layers have failed. It also shows transverse failure along the hole 
edge partway between the net-section and x axis. Because of the large number of 
layers failed it is expected that inplane shear failure dominates in the two half-wave 
configuration as it did in the one half-wave configuration, Fig. 5-15. The highest value 
of the interlaminar shear failure parameters are /? 13 = 0.35 and R n = 0.42 at the load 
level of P/EA = 0.00394. The final failure is expected to begin at the area outside the 
simple support slightly away from the net-section, and propagate into the plate to- 
ward the net-section and hole edge. 
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(Number in the figure represents the no. of failed layers.) 


Load: 0.00394 

EA 


Fig. 5-16 Predicted failure locations and failure modes of the [ ±45]4s plate with a hole In the 
two half-wave configuration. 
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5.2 Failure of Solid Plates 


5.2.1 The Solid [±45/0/90] 2S Plate 


As an example of the characteristics of the solid plates, the deflected shapes of 
the solid quasi-isotropic plate for one half-wave and two half-wave configurations are 
shown in Fig. 5-17 and Fig. 5-18, respectively. The load levels used in the figure are 
about a factor of 6 and a factor of 4 times the first and the second buckling loads, re- 
spectively. Comparing the deflected shape of the one half-wave configuration for the 
solid plate with the deflected shape of the one half-wave configuration of a plate with 
a hole, Fig. 5-3, no reverse curvature appears in the solid plate. It is felt that the 
stiffness at the center of the plate, for the case of no hole, is responsible for pre- 
venting this reverse curvature. For the cases of the two half-wave configuration, the 
plate with a hole, Fig. 5-4, and the plate without a hole are quite similar. 

Figure 5-19 shows the load-end shortening relations of the solid [±45/0/90] 2S 
plate. As before, the failure load of each failure mode is marked on the load-end 
shortening relation for the one half-wave and two half-wave configurations. It can be 
seen that for both one half-wave and two half-wave configurations all three failure 
modes occur at almost the same load level of about 6 times the lowest buckling load. 
Considering that for the plate without a hole all of three failure loads are close each 
other, the failure process is expected to be sudden and catastrophic. Referring to 
Table 5-4 and Fig. 5-2, for the one half-wave configuration as the plate is loaded from 
zero load, the first transverse failure occurs in element 50 in the center of a quarter 
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Fig. 5-17 Shape of [ ± 45/0/90]k solid plate deflected In the one half-wave configuration at 
P/EA = 0.00289. 
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Fig. 5-18 Shape of [ ± 45/0/90] 3S solid plate deflected in the two half-wave configuration at 
P/EA= 0.00269. 
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Fig. 5-19 Failures and load-end shortening relations for the solid [±45/0/90] 2S plate. 
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plate. The first inplane shear failure then occurs at the corner of the clamped edge 
and simple support, element 62. Fiber failure occurs just inside the simple support 
toward the clamped end, element 69. This is a location where curvatures begin to 
develop in the x direction. The fiber failure loads for one half-wave and two half-wave 
configurations are greater than the fiber failure loads of corresponding configurations 
of the quasi-isotropic plate with a hole. It is interesting to note that the first fiber fail- 
ure load for a plate without a hole deformed in the two half-wave configuration is the 
same as the first fiber failure load of the plate with a hole deformed in the one half- 
wave configuration, namely P/EA = 0.00256. The failure locations, however, are dif- 
ferent. 

Figure 5-20 shows the damage accumulation scenario in the solid plate at a load 
level slightly above the first fiber failure load. The extension area near the net- 
section outside the simple support experiences severe inplane shear stress. From 
the figure it can be seen that half of the layers of the laminate are damaged there. 
Recalling that the quasi-isotropic plate with a hole shows severe failure at the net- 
section hole edge and along the net-section, Fig. 5-6, it is expected that for the one 
half-wave configuration the failure pattern of the solid plate will be different from the 
failure pattern for the plate with a hole. Regarding interlaminar stress, for the solid 
plate in the one half-wave configuration the highest values of interlaminar shear fail- 
ure parameters at the load level of P/EA = 0.00289 are F? 13 = 0.33 and R 2J = 0.37. These 
values show that the interlaminar shear stresses are low in the solid plate which is 
deformed in the one half-wave configuration. 

For the two half-wave configuration failure first occurs as transverse failure at 
roughly the center of the quarter plate, specifically in element 48. Inplane shear fail- 
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Load : TT“= 0.00289 
EA 


Fig. 5-20 Predicted failure locations and failure modes of the solid [ ± 45/0/90]j S plate In the 
one half-wave configuration. 
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ure then follows in the element 76, outside the simple support and slightly away from 
the net-section. Fiber failure then occurs at the net-section just inside the simple 
support. As shown in the damage scenario in Fig. 5-21, the extension area outside the 
simple support slightly away from the net-section experiences inplane shear stress. 
The location of the first fiber failure in this solid plate in the two half-wave configura- 
tion is the same location as first fiber failure of plate with a hole in the two half-wave 
configuration, namely at the simple support "net-section", net-section being in quotes 
since it does not really apply to the plate without a hole. For both the one half-wave 
and two half-wave configurations of the solid plate it is expected that the final failure 
will occur in the extension area due to severe inplane shear. The center of the plate 
remains intact for both the one half-wave and two half-wave configurations. The 
interlaminar shear stresses are low for the two half-wave configuration. The highest 
values of R, 3 and R a are 0.29 and 0.34 at the load level slightly above the first fiber 
failure load, respectively. 


5.2.2 The Solid [±45/0 2 ] 2S Plate 


In Fig. 5-22 the load-end shortening relations for the solid [ + 45/0 2 ] 2S plate for 
both one half-wave and two half-wave configurations are shown. For both the one 
half-wave and the two half-wave configurations failure occurs in the sequence of 
transverse failure, inplane shear failure, and then fiber failure. Between the two con- 
figurations the failure loads of each failure mode are about the same. However, for 
the one half-wave configuration failure loads are slightly higher. In general, trans- 
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Fig. 5-21 Predicted failure locations and failure modes of the solid [ ± 45/0/90] 2S plate in the 
two half-wave configuration. 
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Fig. 5-22 Failures and load-end shortening relations for the solid [ ± 45/0 plate. 
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verse failure occurs at about 4 times the buckling load, inplane shear failure occurs 
at about 6 times the buckling load, and fiber failure occurs at about 6 times the 
buckling load. 

For the one half-wave configuration transverse failure first occurs in element 48, 
near the center of the quarter-plate planform. First shear failure occurs at the net- 
section outside the simple support, element 73. First fiber failure occurs near the 
clamped edge but along the simple support, element 69. First fiber failure occurs at 
the same location as the solid quasi-isotropic plate in the one half-wave configura- 
tion. Figure 5-23, the damage accumulation scenario, shows that the solid 
[±45/0 2 ]k plate in the one half-wave configuration experiences inplane shear failure 
in the extension area near the net-section at the load level slightly greater than the 
first fiber failure load. Inside the simple support rather widespread inplane shear 
failure appears but it is less severe than that in the extension area near the net- 
section, i.e., less layers are failed. Transverse failure is widespread in the center of 
the quarter plate. The interlaminar shear failure parameters at the load level of 
P/EA = 0.00229, slightly above the first fiber failure load, are = 0.39 and R 23 — 0.37. 

For the two half-wave configuration first fiber failure occurs at the same location 
as with the quasi-isotropic solid plate deformed in two half-waves, i.e., in element 64, 
just inside the simple support on the net-section. This first fiber failure follows 
transverse failure in element 58 and inplane shear failure in element 76, the latter 
location being outside the simple support. Figure 5-24 shows that at a load level 
slightly above the first fiber failure the area on both sides of the simple support and 
slightly away from the net-section has significant damage accumulation due to in- 
plane shear failure. The interlaminar shear failure parameters at the load level 
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y (Number in the figure represents the no. of failed layers.) 

Load: -= 7 -= 0.00229 
EA 


Fig. 5-23 Predicted failure locations and failure modes of the solid [ ± 45/02 ]*s plate in the one 
half-wave configuration. 
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Load: -= 7 -= 0.00223 
EA 


Fig. 5-24 Predicted failure locations and failure modes of the solid [ ± 45/02] 2S plate In the two 
half-wave configuration. 
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P/EA = 0.0023, slightly above the first fiber failure load, are Rn = 0.30 and R 23 = 0.32. 
From Fig. 5-24 it can be deduced that the final failure of this plate in two half-waves 
will occur slightly away from the net-section and along the simple support. 


5.2.3 The Solid [±45/0 6 ] s Plate 

The load-end shortening relations for the [ ± 45/Osjs plate are shown in Fig. 5-25. 
The failure load for each mode is very close to being the same for both the one half- 
wave and two half-wave configurations. In general, transverse failure occurs at about 
4 times the buckling load and inplane shear failure and fiber failure occur at about 6 
times the buckling load for both configurations. For the one half-wave configuration 
transverse failure first occurs at the net-section away from the simple support, in el- 
ement 46. Inplane shear failure occurs along the simple support in element 68, be- 
tween the net-section and the clamped edge. First fiber failure also occurs at exactly 
the same location. This is obviously a location of severe stress. Figure 5-26 shows 
the damage accumulation at a load just above the first fiber failure load in the one 
half-wave configuration. It shows that the transverse failure has developed near the 
net-section. The inplane shear failure and fiber failure continue to develop in the 
same location where the first failure for each mode is predicted. The interlaminar 
shear failure parameters at just above the fiber failure load are ft 13 = 0.35 and 

Rzi — 0.26. 

For the two half-wave configuration transverse failure first occurs toward the 
center of the quarter plate, in element 58. Inplane shear failure first occurs along the 
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Fig. 5-25 Failures and load-end shortening relations for the solid [ ± 45/0 b] s plate. 
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Load: -= 7 -- 0.00169 
EA 


Fig. 5-26 Predicted failure locations and failure modes of the solid [ ± 45/0*] s plate In the one 
half-wave configuration. 
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simple support half-way between the net-section and the clamped edge, element 67. 
First fiber failure occurs at the net-section near the simple support, element 64. 
Damage accumulation for the two half-wave configuration is shown in Fig. 5-27. It 
shows that at a load slightly greater than the first fiber failure load the transverse 
failure occurs from the center of the quarter-plate to the simply-supported edge. 
Where there is transverse failure, about a quarter of the layers are failed. Figure 5-27 
shows that inplane shear failure also develops along the simple support in the two 
half-wave configuration. It is expected that the final failure will occur along the simple 
support for [ ± 45/0,] s solid plate in the two half-wave configuration. The failure could 
begin at the net-section. The interlaminar shear failure parameters at the load level 
above the first fiber failure load are = 0.49 and = 0.29. It should be noted that 
the solid [ ± 45/0 e ]s plate, in either configuration, is the strongest plate studied so far. 


5.2.4 The Solid [±45] 4S Plate 


Figure 5-28 shows the load-end shortening relations for the solid [ ± 45]« plate 
for both the one half-wave and two half-wave configurations. For a given failure mode 
the failure load is practically independent of the deformed configuration. The inplane 
shear failure occurs at about 2.5 times the buckling load, transverse failure occurs 
at about 4 times the buckling load, and fiber failure occurs at about 5 times the 
buckling load. For the one half-wave configuration, as the load is increased, inplane 
shear failure first occurs outside the simple support on the net-section, element 73. 
With increasing load transverse failure occurs near the center of the quarter plate, in 
element 50. First fiber failure occurs closer to the clamped support and just away 
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Load: -= 7 -= 0.00177 
EA 


Fig. 5-27 Predicted failure locations and failure modes of the solid [ ± 45 /Qb] s plate in the two 
half-wave configuration. 
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Fig. 5-28 Failures and load-end shortening relations for the solid [ ± 45] 4 s plate. 
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from the simple support, in element 60. This location of first fiber failure being nearer 
the clamped support than with the other three solid plates is quite similar to the case 
of the [ ± 45]« plate with a hole. The [ ± 45]« plate with a hole experienced fiber 
failure nearer the clamped end than any other plate with a hole. Like the plate with 
a hole, it is expected that inplane shear failure will dominate the failure in this solid 
plate. This is supported by Fig. 5-29. Figure 5-29 shows the damage accumulation 
scenario at a load just above the first transverse failure load in the one half-wave 
configuration. It can be seen that the simply supported area near the net-section is 
severely damaged by inplane shear failure, all 16 layers having failed. Concentrated 
inplane shear failure can be observed outside the simple support and just inside the 
simple support. 

For the two half-wave configuration the situation is similar. Figure 5-30 shows 
the damage accumulation in the two half-wave configuration at slightly above the first 
transverse failure load. Transverse failure is minor. However, significant inplane 
shear failure is predicted near the simple support near the net-section. It can be seen 
that all of 16 layers are failed in the area outside the simple support and just inside 
the simple support. Comparing with the one half-wave configuration, the location of 
inplane shear failure is moved slightly away from the net-section. The interlaminar 
shear failure parameters are very small. For the one half-wave configuration the 
highest values of ft 13 and are 0.17 and 0.20, respectively, at the load level slightly 
above the first transverse failure load. For the two half-wave configuration R n is 0.18, 
and R 23 is 0.22. 
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Fig. 5*29 Predicted failure locations and failure modes of the solid [ ± 45] w plate in the one 
half-wave configuration. 
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Fig. 5-30 Predicted failure locations and failure modes of the solid [ ± 45] 4S plate In the two 
half-wave configuration. 
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5.3 Summary 


One of the most important conclusions from the analysis of failure in this chapter 
is that it appears as though the failure load may be independent of whether the plate 
deforms in the one half-wave configuration or the two half-wave configuration. An 
examination of the results presented shows this to be the case. 

In this chapter no interlaminar shear failure is predicted for either the one half- 
wave or two half-wave configuration in any plate. Bringing forward the results from 
Ch. 4 and the results from this chapter, another important conclusion can be made. 
In Fig. 4-12 it was shown that the interlaminar stress component r XI was much larger 
at the net-section hole edge than at the net-section simple support region. This figure 
was for the case of the two half-wave configuration and the quasi-isotropic laminate 
but the findings are general. As a result, it can be concluded that if the plate is de- 
formed in the two half-wave configuration and if failure is due to a high interlaminar 
shear stress component r„, failure will begin at the hole edge, not elsewhere. This 
conclusion should be useful when observing the location of failure in the exper- 
iments. It is tempting to say that if failure begins at any location other than the hole 
edge, interlaminar failure due to t„ is not involved. However, the plate could fail at 
the hole edge due to the other mechanisms, e.g., fiber failure when the plate is de- 
formed in the one half-wave configuration as predicted in this chapter (see Table 5-2 
and Fig. 5-6). 

Finally, the failure load of the [ ± 45/0 e ]s plate is expected to be the highest, the 
failure load of the [ ± 45/0 2 ] 2S plate next, and the failure load of the [ ± 45/0/90] 2S plate 
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next. The failure load of the [ ± 45]4s plate is expected to be the lowest. This ordering 
of failure loads is true for both solid plates and plates with holes. Only the 
[ ± 45/0/90] 2S plate with hole deforming in the one half-wave configuration is expected 
to fail at the net-section hole edge. For the other plates failure is predicted to occur 
near the simple support. The postbuckling failure loads of solid plates are expected 
to be higher than the postbuckling failure loads of the plates with holes with 
D/W = 0.3. 


Attention now shifts to the experimental aspects of this study. Correlation of the 
predictions of this and past chapters with experimental observations and measure- 
ments is critical to providing credibility for the findings. 
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6.0 Description of Experiment 


6.1 Test Fixture 


Considerable effort has been put into developing a fixture to load composite 
plates in the manner discussed in the preceding chapters. The boundary conditions 
and other conditions used in the analysis, and described in Fig. 1-1, were actually 
idealizations of the conditions that could be most conveniently built into an exper- 
imental set-up. It is more difficult to enforce a specified stress resultant distribution 
on the loaded edge than to enforce a specified displacement distribution. Hence the 
uniformly displaced edge is used. In addition, it is easier to clamp this displaced edge 
than to allow it to be simply supported. Similarly, as has been mentioned, it is im- 
possible to not have the plate width exceed the supported width. The extra width 
outside the supports is necessary to accommodate the plate width decrease when 
the plate buckles. As has been shown, the extra width must be included in the 
analysis. The extensions of the plate outside the supports provide extra stiffness and 
are a site of high inplane stresses. 
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The detail drawings of the fixture are shown from Fig. 6-1 to Fig. 6-9, and a pho- 
tograph of the fixture in the displacement-controlled Tinius-Olsen load frame is shown 
in Fig. 6-10. When placed in the fixture the two opposite simply-supported edges of 
the composite plate specimen are vertical, the displaced clamped edge is at the top, 
and the stationary clamped edge is at the bottom. 

The fully assembled fixture with the composite plate specimen is depicted in Fig. 6-1. 
The fixture consists of a top assembly, a bottom assembly, a stiffback steel bar, and 
a steel ball. The majorities of both the top and bottom assemblies are made of alu- 
minum and the top assembly moves vertically relative to the bottom assembly 
through the four linear bearings. 

The top assembly, shown in Fig. 6-2, includes a top plate, four hardened linear 
bearing shafts, of steel, for alignment of the top assembly with the bottom assembly, 
and the supports for clamping the top edge of the composite plate specimen. The 
steel shafts are 0.75 in. in diameter. Other dimensions are shown on the figure. 

The bottom assembly, illustrated in Fig. 6-3, consists of a bottom frame (which 
consists of a bottom plate, side gussets, and side plates) four linear bearings, four 
knife edges to effect the vertical simply supported edges, brackets for holding the two 
linear variable differential transducers (LVDT's), and clamped supports for the bottom 
edge of the composite plate specimen. More details of the bottom assembly are 
shown in Fig. 6-4, an illustration of the bottom frame. 
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Fig. 6-1 Full assembly of the fixture with plate specimen. 




Top plate 
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Fig. 6-2 Top assembly of the fixture. 




Pinned knife edge 
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bottom plate 
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Fig, 6-4 Bottom 








Included in the top and bottom assembly are machined hardened steel inserts. 
The hardened steel inserts are actually what the top and bottom machined edges of 
the composite plate specimen bear against in the top and bottom assemblies. These 
inserts are flat to within close tolerances (less than 0.0005 in.) so as uniform as pos- 
sible edge displacements are imposed on the top and bottom edges of the composite 
plate specimen, the bottom edge displacement, of course, being zero. In addition to 
maintaining flatness, the hardened inserts prevent the sharp edges of the composite 
plates from deforming the softer aluminum at high load levels. 

With the bottom assembly, though not shown in detail, but illustrated to some 
degree in Fig. 6-5, one of the clamped supports is afixed to the bottom plate with pins 
and Allen bolts. The other clamped support contained slotted holes so it can move 
in a horizontal plane on the bottom plate. This sliding clamp support is forced against 
the composite plate, which is, in turn, forced against the fixed clamp support. The 
sliding clamp support is forced against the composite plate by screws which is 
threaded into a supporting block attached to the bottom plate. The supporting block 
is shown in Fig. 6-6. In essence, the sliding clamp support squeezs the composite 
plate, in the thickness direction, against the fixed clamp support to enforce the 
clamped boundary condition. Once the plate is sufficiently "squeezed", both clamp 
supports are secured to the bottom plate with Allen bolts. 

Referring to Fig. 6-3, knife-edges are fastened to the side plates, the side plates 
also serving to hold the linear bearings that guide the hardened steel shafts from the 
top assembly. Though the horizontal forces, either left-right or fore-aft, are expected 
to be small, the side plates are reinforced with side gussets. Brackets for the LVDT's 
are fastened to the side gussets. The scheme for forcing the composite plate against 
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dimensions in inches 



the vertical knife-edges is similiar to the scheme for clamping it on the bottom edge. 
A schematic of a knife-edge is shown in Fig. 6-7. On both the left and right sides of 
the composite plate one knife-edge is fixed in place while the other is movable and 
can slide against the composite plate, in the thickness direction. The sliding knife 
edge is forced against the plate by another supporting block and Allen screw as- 
sembly, this supporting block component being shown in Fig. 6-8. 

The fixture is designed so the bottom fixed clamp support and the two fixed knife 
edge supports form a boundaries in a plane against which the composite plate is 
forced. That the plate is forced against three boundaries all lying in a plane mini- 
mizes eccentricities or misalignments of the plate with respect to the loading. 

The clamp supports on the top edge, supports that are part of the top assembly, 
are both sliding , much like the one bottom clamp support. One sliding support is 
positioned, by hand, against the composite plate which is already forced against the 
three other planar boundaries. This sliding support is secured with Allen screws. 
The other sliding boundary is forced against the composite plate, as at the bottom 
edge, by a screw which is threaded into the supporting block attached to the top 
plate. 

The LVDT bracket, shown in Fig. 6-9, clamps the round barrel of the LVDT and 
prevents motion of the LVDT relative to the bottom assembly. The clamping action 
is such that the LVDT can easily be moved vertically to balance the LVDT voltage 
output. The LVDT's are a vital part of the experiment. The average of the left and 
right LVDT's is used as a measure of the displacement of the loaded edge of the 
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composite plate. A comparison of the left and right displacement readings is a direct 
measure of how uniformly the edge is being displaced. 

The load from the load frame is transmitted to the fixture through a spherical 
steel ball on the top assembly. A thick stiffback steel bar between the ball and the 
top plate of the top assembly provides additional stiffness and helps enforce the 
condition of uniform inplane displacement. The steel ball provides a decoupling of 
any misalignments of the fixture relative to the Tinius-Olsen load frame, though any 
misalignments were expected to be slight. The linear bearings restrains all but ver- 
tical motions of the top assembly relative to the bottom assembly. Replacing a 
specimen is relatively easy. The pinned clamp support on the bottom edge and the 
pinned knife edge do not have to be loosened, nor does the one movable clamp 
support on the top edge. Only the movable components needed to be loosened to 
accomodate a new specimen. 

The fixture is designed to accomodate an 11 in. by 11 in. square plate with 0.08 
in. thickness. One-half inch of the plate on each edge is used for the various boundary 
conditions. Thus the analysis region is 10 in. by 10 in.. The 0.5 in. on the left and right 
edges of the plate, outside the knife-edges, are simply free. The 0.5 in. at the top and 
bottom are under clamps. The 10 in. by 10 in. dimensions were somewhat dectated 
by the size of the load frame available. Also, the dimensions coincide with the di- 
mensions used by the other investigator [31] and hence results could be compared. 
The loaded edges of the laminates are machined precisely to achieve uniform contact 
with the loading surface of the fixture, namely the hardened steel inserts. In Fig. 6-10 
the fixure with a composite plate is shown within the Tinius-Olsen load frame. The 
LVDT's are in place in their brackets. Wires from strain gauges (to be discussed next) 
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Fig. 6-10 Close-up of the fixture In load frame. 
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are visible. As will be discussed, the composite plate specimen is painted white on 
one side to enhance moire interferometry. 


6.2 Test Set-up 


As has been discussed, LVDT's are used to measure the inplane movement of 
the top end of the composite plate. The LVDT's used had sensitivity of 195 volts per 
inch and were manufactured by Schaeitz. 

Strain gauges are bonded at a number of points on the composite plate specimen 
to measure local response. The strain gauge response could be used to examine 
plate behavior, and also evaluate the validity of the test by identifying lack of sym- 
metry, unusal strain response, etc. 

The locations where the strains are measured are depicted in Fig. 6-11. Strains are 
measured at six locations on one side of the plate, gauges 1-6, and two locations on 
the other side, gauges 7 and 8. Gauges 7 and 8, along with gauges 2 and 5, respec- 
tively, provide back-to-back strain measurement. Gauges 1 and 3 are placed to 
measure left-right uniformity of the response at the clamped end, gauges 2 and 7 are 
used to measure bending at the clamped end. Gauges 4 and 6 are also used to 
measure left-right uniformity. Gauges 5 and 8 are used to measure bending effects 
at the hole edge. 
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dimension in inches 



View toward the black side 


Fig. 6*11 Strain gauging locations on the plate specimen. 
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After strain gauging, the side of the composite plate which has only two strain 
gauges is painted white for use with the shadow moire measurements. The shadow 
moire [57] is a technique used to observe the overall out-of-plane displacement pat- 
tern of plates. For shadow moire measurements a clean white surface gives better 
results. Thus the side with less strain gauges is chosen for painting white. For im- 
plementing the shadow moire technique a reference grid and light source are 
needed. Here the reference grid is a thin photographic film with 200 lines per inch 
mounted on a one-half inch thick piece of plexiglass. The plexiglass used has close 
tolerence on both its thickness and flatness. The lines of the grid are parallel with the 
simply supported edges (vertical). The shadow moire fringe pattern is recorded on 
video tape by means of a TV camera. Still photographs of the fringe patterns are also 
taken during the test. For the quantitative measurement of out-of-plane displacement 
at a specific point another LVDT is used. This LVDT is placed perpendicular to the 
surface of the plate. The simplified top view of the arrangement of the LVDT to 
measure the out-of-plane displacement and the shadow moire set-up is shown in Fig. 
6-12. Finally, acoustic emission signals are recorded. An acoustic emission sensor 
is attached on the side of the composite plate with six strain gauges. The acoustic 
emission sensor needs to be attached securely so that it does not fall off during the 
test. Since the composite plate bends during loading, it is rather difficult to maintain 
a good contact between the sensor and the plate. When a mode shape changes dur- 
ing postbuckling, there is a sudden motion of the plate. If the sensor is not attached 
securely, it can fall off. 

The voltages corresponding to composite plate end shortening, out-of-plane dis- 
placement, strain, and load are collected by a data acquisition system at regular time 
intervals. These voltages are converted to engineering units by being multiplied by 
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Fig. 6-12 Top view of arrangement of LVDT for out-of-plane displacement measurement and 
shadow moire. 
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the appropriate scale factors, when the data are analyzed. A schematic drawing of 
the overall layout of the test equipment is shown in Fig. 6-13. A photograph of the 
overall set-up is shown in Fig. 6-14. The data acquisition system for strains and loads 
is in the foreground on the left in the upper photograph. The data acquisition system 
for acoustic emission is in the background on the left. The console on the right is the 
control unit for the Tinius-Olsen load frame. The VTR equipment, which is placed in 
front of the of the load frame during testing, is shown in the lower photograph. The 
next chapter discusses the results obtained with this experimental set-up and the 
comparison of the experimental results with numerical predictions. 
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Load Frame 




Fig. 6-13 Equipment layout for experimental test, 
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(a) Test equipment -load frame, data acquisition system, A.E. equipment. 



(b) Video recording equipment. 


Fig. 6-14 Photograph for experimental set-up. 
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7.0 Experimental Results and Comparison with 


Numerical Results 


As shown in Table 7-1, two plates for each laminate were tested to failure. With 
these plates the even numbered plates were tested first, and the odd numbered 
plates were tested second. For each test, plate end-shortening, the out-of-plane de- 
flection of a point at net-section hole edge, strains at 8 locations on the plate, and the 
applied load were measured and recorded. Acoustic emission test and moire- 
interferometry measurements were also made. These measurements are discussed 
in this chapter, and they are compared with the predictions discussed earlier. 

In order to check the uniformity of the loading on the edge of the plate, for each 
of the plates tested the inplane displacements on the left and right side of the plate, 
namely, the responses of LVDT 1 and LVDT 2 in Fig. 6-11, were compared. The dis- 
placements used in the load vs. end-shortening relations which are compared with 
numerical results are the averaged values of these left and right LVDT responses. In 
the load vs. out-of-plane deflection relations, since the out-of-plane deflection was 
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measured at the hole edge of the net-section, comparison with numerical results can 
be made only for one half-wave configuration. Because the net-section is the nodal 
line of deformed shape, for the two half-wave configuration the theoretical out-of- 
plane deflection at this point is zero. For both types of load vs. displacement re- 
lations discussed, the numerically predicted failure loads are marked on the relation. 

For each laminate shadow moire patterns are compared with numerically calcu- 
lated out-of-plane deflection contour drawings. Finally, acoustic emission results are 
plotted in the format of end-displacement vs. accumulated counts. In this manner the 
acoustic emission data can be compared easily with the load-end shortening relation 
by drawing the two relations on one figure. Numerically predicted failure loads are 
also marked on this figure. In the following discussions the comparison of exper- 
imental results with numerical prediction will be made on a laminate by laminate 
basis. 


7.1 The [±45I0I90] 2S Plates 


7.1.1 Load vs. Deflection Relations 

Experimental load vs. end shortening relations for the [ + 45/0/90] 2 s plates, plate 
1 and plate 2, and a comparison of these relations with numerical predictions is 
shown in Fig. 7-1. Failure information, both experimental and numerically predicted, 
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Fig. 7-1 Comparison of load vs. end shortening relations and failure loads for the 
[ ± 45/0/903,8 Plates. 
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are shown, as is the level of the predicted buckling load. The predicted buckling load 
for this plate is, as was indicated in Table 2-3, P„/EA = 0.000411 , or = 258 lb/in.. 
The experimental buckling load was obtained from the experimental load vs. end 
shortening relations by drawing straight lines on the prebuckling and postbuckling 
path and finding the intersection of these two straight lines. The experimental 
buckling load of plate 1 was P Cf /EA = 0.000375 , or N? = 235 lb/in., and that of plate 2 
was p cr /EA = 0.000364 , or N" = 228 lb/in.. These loads are summarized in Table 8-1 
along with similar information for plates to be discussed. Chapter 8 will furthur dis- 
cuss the findings presented in this chapter. Referring to the figure, after the buckling 
load the experimental relations followed the initial postbuckling path characterized 
by one half-wave in both the x and y directions. At the load level of P/EA = 0.0015 ( 
N x = 940 lb/in. ) the mode shape in the x direction (the loading direction) suddenly 
changed to two half-waves, while the mode shape in the y direction remained in one 
half-wave. The experimental relations for both plate 1 and plate 2 showed mode 
transition from the one half-wave to the two half-wave configuration at similar load 
levels. As the load level approached the mode shape transition load, the mode 
transition began to develop from the bottom stationary clamped support. Suddenly, 
the mode shape transition occurred with a rather violent motion. With the first test for 
plate 2 the moire reference grating, which was placed in front of the plate without any 
support, was knocked out of the alignment by the violence of mode shape transition. 
The test had to be stopped and the plate unloaded. Upon repositioning of the moire 
reference grating with a support and reloading the plate, the mode shape transition 
occurred again at the same load level. The mode shape transition for the plate was 
quite suprising. The problem of predicting the mode shape transition is beyond the 
scope of this study. In this study the out-of-plane deformation shape before the mode 
shape transition is assumed to be in one half-wave in both the x and y directions, and 
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after the mode shape transition two half-waves in the x direction and one half-wave 
in the y direction. This assumption will be examined and justified in the following 
discussions by comparing the experimental results with numerical results in terms 
of load-end shortening relations, strains, and out-of-plane displacements. 

Referring to Fig. 7-1, the experimental results indicate that the stiffness after the 
mode shape transition was lower than the stiffness before the mode shape transition. 
The data for plate 2 showed slightly softer behavior as compared with plate 1, and 
as compared with numerical predictions. It was found that this was due to paint on 
the loaded edge of the plate. When the paint was sprayed onto the plate to enhance 
the shadow moire, the loaded edge was also sprayed. After a few tests it was realized 
that the excessive paint on the loaded edge affected the stiffness observed, and that 
when the paint was sprayed, masking the edges of a plate was necessary. 

Until the mode shape change occurred, the numerical prediction for one half- 
wave configuration agreed quite well with the experimental results for plate 1. Just 
after the mode shape change the experimental results for plate 1 agreed with the 
numerical prediction for two half-waves. At higher load levels the numerical predic- 
tion began to deviate from the experimental results. The experimental relation 
showed softer behavior than the numerical prediction. It is felt that the plate became 
soft, because the plate began to fail. As the plate began to fail, it began to lose load 
bearing capability. However, from the relations in Fig. 7-1 it can be seen that until the 
final failure occurred, the load bearing capacity of the plate did not decrease sub- 
stantially. Final failure occurred catastrophically for each plate at the load level 
marked by the horizontal bars. Plate 1 failed at a load of P/EA = 0.00234 ( N x = 1468 
Ib/in. ) and plate 2 failed at P/EA = 0.00271 { N„ = 1700 lb/in.). The numerically pre- 
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dieted failure loads are marked on the numerical load-end shortening curves. For the 
two half-wave configuration the numerically predicted inplane shear failure and fiber 
failure loads occur at P/EA = 0.00244. These failure loads are summarized in Table 
8-2 along with similar information for plates to be discussed. It can be seen that the 
experimental failure load of plate 1 was close to the numerically predicted inplane 
shear and fiber failure loads, while the failure load of plate 2 was higher than that of 
plate 1. The average of the failure loads of plate 1 and plate 2 is P/EA = 0.00253. With 
the numerically predicted fiber failure load considered as the failure load, the devi- 
ation of the numerical prediction from the average of the experimental result is 3 %. 

Figure 7-2 shows the degree of uniformity of the loading for plate 1 and plate 2. 
By comparing the displacements measured by LVDT 1 and LVDT 2, the left-right uni- 
formity of loading was checked. It can be seen that until the mode shape change oc- 
curred, both plates 1 and 2 were loaded uniformly, the left and right displacements 
being identical. After the mode shape change there was a slight difference between 
the left and right displacements. With plate 1 the side of the plate with LVDT 2 dis- 
placed more, and with plate 2 the side of the plate with LVDT 1 displaced more. This 
shows that the fixture was not biased. The discrepancy between the two side dis- 
placements was worse for plate 2 than plate 1. However, considering the magnitude 
of the displacements, the discrepancy was not serious. 

The experimentally measured and numerically predicted out-of-plane deflections 
of the plates at the net-section hole edge are shown in Fig. 7-3. The jump from the 
one half-wave configuration to the two half-wave configuration is very evident. Until 
the jump occurred, the experimental load vs. out-of-plane deflection relation agreed 
well with the numerically predicted relation. The almost flat part of the load vs. out- 
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Fig. 7-2 Comparison of displacements measured from LVDT 1 and LVDT 2 for plate 1 and 
plate 2. 
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Fig. 7-3 Comparison of load vs. deflection relations and failure loads for [ i 45/0/90]; 
plates. 
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of-plane deflection curve representing the mode shape change shows that the change 
occured in a very short period of time. Theoretically, the out-of-plane displacement 
at the net-section hole edge for two half-waves is zero. This is because the net- 
section hole edge is a node line for the two half-wave configuration. The exper- 
imental out-of-plane displacement after the mode shape change was close to zero, 
but it was not identically zero. There are a few reasons why it was not zero. First, the 
node line of the two half-wave configuration was not perfectly parallel to y axis. As 
shown in Fig. 7-2, it was impossible to load the plate in a perfectly uniform manner 
throughout the range of load. Second, with the sudden movements associated with 
the mode shape change, there is a risk that the LVDTs' may have shifted slightly 
when the change occurred. 


7.1.2 Load vs. Strain Relations 


Experimentally measured strains at the 8 locations on the plate 1 are shown in 
Fig. 7-4. The strain gauge locations were illustrated in Fig. 6-11, and the strain gauges 
are grouped in Fig. 7-4 by location on the plates. Gauges 1, 2, 3, and 7 are located at 
the clamped movable edge and gauges 4, 5, 6, and 8 are at the net-section. Exam- 
ination of each group of gauges provides information on the symmetry of the loading, 
when buckling began, local bending effect, etc. In Fig. 7-4 it can be seen that the 
strains after the mode shape change increased faster with load than before the mode 
shape change. This reflects the sharper curvatures associated with the two half-wave 
configuration. Strains 1 and 3 were very close. This indicates good but not perfect 
left-right symmetry of loading, boundary conditions, and response near the movable 
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Fig. 7-4 Experimentally measured strains of plate 1. 
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clamped support. These strains look identical despite the fact that the inplane dis- 
placements on the left and right side of the plate had slight differences, i.e., see Fig. 
7-2. Strain gauge 2 and strain gauge 7 were back-to-back gauges. Up to the buckling 
load both strain gauges registered compressive strains. After buckling the response 
of strain gauge 2 remained in compressive and continued to increase. However, the 
response of strain gauge 7 changed sign and became positive. Before explaining 
more on this, another back-to-back pair of strain gauges, strain gauge 5 and strain 
gauge 8, is examined. Until buckling occurred, they also both showed compressive 
strains. After buckling strain 8 remained in compression, but strain gauge 5 became 
tensile. The high level of compressive strain in strain gauge 8 reflects bending and 
strain concentration effects at the net-section hole edge. That strain 5 was positive 
indicates that the plate was deflected so that strain gauge 5 was on the convex sur- 
face of the plate, as in Fig. 7-5. Strain gauge 2 and strain gauge 5 were on the same 
side of the plate, and strain gauge 7 and strain gauge 8 are on the other side of the 
plate. It is interesting that strain 2 and strain 5 had signs opposite each other, and 
strain 7 and strain 8 had signs opposite each other. This can be explained by Fig. 7-5. 
Since the loaded edge was clamped, the curvature near the clamped edge was op- 
posite in sign to the curvature at the center of the plate. This made the signs of strain 
2 and strain 5, both being on the same side of the plate, opposite one another. The 
fact that the strains at location 5 and 8 dropped so dramatically after the mode shape 
change was due to the fact that these gauges were on the two half-wave configuration 
nodal line. With the two half-wave configuration there was not much strain for these 
gauges to register, whereas in the one half-wave configuration there were significant 
bending strains at the locations of gauges 5 and 8. 
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Fig. 7-5 Schematic drawing to explain the sign change of strains. 
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Like strain gauges 1 and 3, strain gauges 4 and 6 were not influenced much by 
bending, since the plate was held flat at those gauge locations by the simple sup- 
ports. It should be mentioned that gauges 4 and 6 generally registered higher strains 
than gauges 1 and 3, even though they were essentially in line with each other and 
near the simple supports, the supports keeping the plate flat. The higher response 
of gauges 4 and 6 reflects the fact that the presence of the hole caused more load to 
be transferred through the material near the net-section simply supported area than 
near the clamped end simple support area. Finally, it can be seen that before the 
mode shape change strains 4 and 6 were almost the same, indicating symmetry of 
response. After the mode shape change the symmetry of left-right strains was not 
as good as before the mode shape change. 

The strains of plate 2, shown in Fig. 7-6, show results similar to the results of 
plate 1, except that the left-right symmetry of the response after the mode shape 
change was not as good as it was for plate 1. 

In the following figures comparisons between the experimentally measured 
strains and numerically predicted strains are shown. This, of course, is a critical 
comparison as it indicates whether or not the numerical analyses discussed in pre- 
vious chapters accurately represent experiments. Figure 7-7 shows the comparison 
for strains 2 and 7 of plate 1. Overall, the figure shows very good agreement between 
experimental results and numerical predictions. Strain gauge 2 registered 
compressive strain and followed the one half-wave numerical prediction until the 
mode shape change. After the mode shape change the strain gauge response fol- 
lowed the numerical two half-wave prediction quite closely. Comparisons for strain 
gauge 7 are similar. Until the mode shape change strain gauge 7 followed the one 
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Fig. 7-6 Experimentally measured strains of plate 2. 
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Fig. 7-7 Comparison of strains near the clamped support of plate 1. 
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half-wave prediction, then it followed the two half-wave prediction, the response be- 
ing tension for the majority of the load range. The comparison for strain gauge 5 and 
strain gauge 8 are shown in Fig. 7-8. Comparisons between numerical prediction and 
experimental measurements for strain gauges 5 and 8 are remarkably close. Initially 
both strain gauges registered compression. After buckling strain gauge 5 registered 
tension, while strain gauge 8 continued to register compression. The gauge readings 
and the numerical predictions are almost in perfect correlation. After the mode shape 
change the gauges followed the two half-wave prediction. The comparisons for strain 
gauges near the simple support, strain gauges 1, 3, 4, and 6, are shown in Fig. 7-9. 
Before changing modes correlations between experimental results and numerical 
predictions are good. After the mode change the experimentally measured strains 
were higher than numerical predictions. After the mode shape change strain gauge 
4 registered considerably higher strain than the numerical predictions. In general, 
though, the strain gauges followed the one half-wave prediction until the mode shape 
changed. Then they followed the two half-wave prediction. The comparisons of the 
strains for plate 2 are shown in Fig. 7-10 through Fig. 7-12. The trends are similar to 
plate 1. For strains 2 and 7 and 5 and 8 excellent correlation between observations 
and predictions can be observed. Again, strain gauge 4 registered considerably 
higher strain than the numerical prediction. 


7.1.3 Shadow Moire Results 


A shadow moire fringe pattern represents lines of equal out-of-plane displace- 
ment. In fig. 7-13 the fringe pattern at a load level before the mode shape change is 
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Fig. 7-8 Comparison of strains at the net-section hole edge of plate 1. 
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Fig. 7-9 Comparison of strains near the simple support of plate 1. 
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Fig. 7-10 Comparison of strains near the clamped support of plate 2. 


Experimental Results and Comparison with Numerical Results 


206 










P/EA 






Fig. 7-12 Comparison of strains near the simple support of plate 2. 
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Fig. 7-13 Predicted and measured shadow moire fringe pattern for a [ ± 45/0/90] 2S plate with 
a one half-wave configuration. 
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compared with numerically predicted out-of-plane displacement contour plot for the 
one half-wave configuration at the same load. Note the voltmeter to the left in the 
lower photograph and the indication that 1 volt = 60,000 lbs. Given that the width of 
the plate was 11 in., this is interpreted to mean 1 volt corresponded to 
A/„ = 5450 Ib/in. . In the figure the voltage reading, 0.12 volts, corresponded to 
N„ = 655 Ib/in. , and a nondimensionalized value of P/EA = 0.00104. This load level 
was well above the buckling load (see Table 2-2 and Fig. 7-3). The out-of-plane de- 
formation shown in Fig. 7-13 is characterized by being predominantly a one half-wave 
configuration. This photo is from the testing of plate 1. Since the view is toward the 
white side of the plate (see Figs. 6-11 and 7-5), the plate in the photo is deflected into 
the page, the positive z direction. The experimentally measured out-of-plane de- 
flection at the net-section hole edge at this load was w/h = 2.7, The numerically pre- 
dicted out-of-plane deflection is w/h = 2.8. These comparisons are shown in the Fig. 
7-13. As the plate was loaded, the mode switched from predominantly a one half- 
wave configuration to a two half-wave configuration. The fringe pattern after the 
mode shape change is compared with the numerically predicted out-of-plane dis- 
placement contour plot in Fig. 7-14. The voltage reading, 0.25 volts, corresponded to 
N* = 1364 Ib/in. , and normalized value of P/EA = 0.00217. Because in the experiments 
the LVDT measuring out-of-plane displacements was placed on the node line, its re- 
sponse was minimal (see Fig. 7-3). Hence it was not possible to correlate displace- 
ments at a point for the two half-wave configuration as it was for the one half-wave 
configuration. 
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Fig. 7-14 Predicted and measured shadow moire fringe pattern for a [ ± 45/0/90]} S plate with 
a two half-wave configuration. 
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7.1.4 Failure and Post-Mortem Analysis 


A photo of a failed [ + 45/0/90] 2S plate is shown in Fig. 7-15. This is plate 1. The 
failure zone can be seen in the photo. It is clear that the failure propagated from the 
simply supported edge inward. The failure occurred near the node line. The left edge 
of the plate just above the node line failed, and the right edge of the plate just below 
the node line failed. For the 2 plates tested this failure pattern was the same. Com- 
pared with the numerically predicted failure locations in Fig. 5-7, the experimental 
results agreed well with the numerical predictions. In the numerical prediction no 
interlaminar shear failure was predicted. Failure was governed by inplane stresses 
at the net-section simple support. This was borne out by experimental observations. 

Figure 7-16 shows the contour plots of the failure parameters defined in eqs.(5.1) 
- (5.5) at selected layers and for each of the five failure modes at the load level where 
the first fiber failure is predicted in two half-wave configuration. These figures are 
drawn with the a posteriori knowledge that the plates failed in the two half-wave 
configuration. The layer shown for each failure mode is the layer exhibiting the 
largest value of the failure parameter for that mode. Figure 7-16(a) shows the high 
fiber direction stresses near the net-section simple support region in the outer bot- 
tom layer, bottom referring to the negative-most z position. These stresses were re- 
sponsible for failure. Figure 7-16(b) shows the transverse failure zone in the central 
region of the quarter plate. Figure 7-1 6(c) shows the inplane shear failure zone along 
the simple support away from the net-section. These could have contributed to fail- 
ure. From Fig. 7-1 6(d) and (e) it can be seen that the failure parameters related to 
interlaminar shear stresses are very low, the maximum value being about 0.35. Fig- 
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Fig. 7-15 Failed [ ± 45/0/90] JS plate in the fixture. 
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(b) Contour plot of R t on 16-th layer from bottom ( 45° layer). 



(c) Contour plot of R n on 16-th layer from bottom ( 45° layer). 

Fig. 7-16 Contour of failure parameters of the [±45/0/90]* 4 plates at P/EA = 0.00255 (contin- 
ued). 
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(d) Contour plot of R n on 10-th layer from bottom ( 0° layer). 



(e) Contour plot of R n on 6-th layer from bottom ( 90° layer). 


Fig. 7*16 Contour of failure parameters of the [±45/0/90]js plates at P/EA — 0.00255 (con- 
cluded). 
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ure 7-16(d) and (e) also shows that the interlaminar shear stresses increase very 
rapidly on the hole boundary. The conclusion to be reached from Figs. 7-16(a)-(e) is 
that inplane not interlaminar stresses were responsible for failure. Interestingly 
enough, referring to Fig. 7-1, if the fiber failure load is considered as the failure load 
of the laminate, then for these quasi-isotropic laminates the numerically predicted 
failure load is about the same, independent of whether the plate is deformed into the 
one half-wave configuration or the two half-wave configuration. 


7.1.5 Acoustic Emission Results 

Finally, acoustic emission test results are shown. The results are presented in 
terms of accumulated counts. The definition of count can be understood from Fig. 
7-17 [58]. Acoustic emission or stress wave emission is the phenomenon of transient 
wave generation due to a rapid release of strain energy caused by a structural alter- 
ation in a solid material [59]. The upper portion of Fig. 7-17 shows a typical example 
of an acoustic emission signal. Count is the number of waves which have a higher 
amplitude than a predetermined threshold amplitude. The acoustic emission test 
results for plate 1 and plate 2 are shown in Fig. 7-18 and Fig. 7-19. In an acoustic 
emission test keeping a good contact between the test specimen and acoustic emis- 
sion transducer is important. Unfortunately, it is not easy in the testing of a large plate 
specimen. In the testing of plate 1 the good contact between the test specimen and 
transducer was lost. Figure 7-18 shows that acoustic emission signals were picked 
up only for the audible mode shape change and final failure. The data for damage 
accumulation between the mode transition and final failure was not collected cor- 
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Fig. 7-17 Definition of terms related to acoustic emission test [58]. 
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Fig. 7-18 Acoustic emission activity and load-end shortening relation 
[±45/0/90] 2S . 
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rectly. Figure 7-19 shows good acoustic emission results. As well as the mode shape 
change and final failure, the damage accumulation can be observed qualitatively from 
the acoustic emission result. Immediately after the mode transition there was little 
acoustic emission. However, with increasing load, acoustic emission increased. This 
correlates with the numerical prediction that transverse failure occurred at a load 
level somewhat higher than the mode transition. However, since the acoustic emis- 
sion results are only qualitative in the application here, the correlation can only be 
suggested. Shortly before failure the acoustic emission count increased dramatically. 
The high acoustic emission count can also be correlated with the numerical predic- 
tion that inplane shear failure and fiber failure was predicted to occur at similar load 
levels. This implies that the final failure would occur catastrophically. Practically 
speaking, the horizontal axis of Figs. 7-18 and 7-19, i.e., the end shortening, has the 
same meaning as time. Therefore, the solid circle on the load-end shortening relation 
representing the final failure falls on the portion of the rapid increase in acoustic 
emission at the final stage of count-end shortening relation. 


7.2 The [+4510 J 2S Plates 


7.2.1 Load vs. Deflection Relations 

Experimental load vs. end shortening relations for the [ ± 45/0 2 ] 2S plates, plates 
3 and 4, and the comparison of these relations with numerical prediction are shown 
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Fig. 7-20 Comparison of load vs. end shortening relations and failure loads for the 
[±45/02 ] 2 s plates. 
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in Fig. 7-20. It can be seen that the results from plate 3 agree well with the numerical 
predictions for both the prebuckling and the postbuckling behavior. The predicted 
buckling load for this laminate was PJEA = 0.000321 ,orWj r = 294 lb/in.. as shown 
in Table 2-3. This load level is marked on Fig. 7-20. The experimental buckling loads 
obtained in the same way as for [ ± 45/0/90]„ plates were P Cf /EA = 0.000300 , or 
Nr = 275 lb/in. for plate 3, and P C ,/EA = 0.000309 , or Wf = 283 lb/in. for plate 4. Until 
the mode transition occurred at P/EA = 0.00138 ( N x = 1265 lb/in. ) in plate 3, the load 
vs. end shortening relation agreed well with the one half-wave configuration numer- 
ical prediction. Just after the mode transition, the load vs. end shortening relation 
agreed well with the two half-wave numerical prediction. At the final stages before 
failure the experimental relation became softer than the numerical prediction. This 
again is felt to occur because the plate began to fail and loose load bearing capability. 
In Fig. 7-20 plate 4 shows softer behavior than plate 3 and the numerical predictions. 
The reason is the same as it was with plate 2. Since plates 2 and 4 were tested first, 
the same mistake, excessive paint on the loaded edge, made this particular plate soft. 
As can be seen from Fig. 7-20, there was a difference in mode transition load between 
plates 3 and 4. The mode transition load of plate 4 was P/EA = 0.001 53 ( N x = 1402 
lb/in. ), while the mode transition load of plate 3 was P/EA = 0.00138, as said earlier. 
However, the failure loads of the two plates were nearly identical, P/EA = 00201 ( 

N x = 1842 lb/in.) for plate 3 and P/EA = 0.00199 ( N„= 1824 lb/in.) for plate 4. Numer- 
ically predicted failure loads for the one half-wave and the two half-wave configura- 
tions are marked on the load vs. end shortening relations. It can be seen that the 
experimental failure loads were very close to the predicted fiber failure load. Again 
it can be seen that the fiber failure load for this plate did not depend greatly on the 
deformed configuration. The numerically predicted fiber failure load is 
P/EA = 0.00206, the experimental failure load of plate 3 was P/EA = 0.00201, and ex- 


Experimental Results and Comparison with Numerical Results 


222 



perimental failure load of plate 4 was P/EA = 0.00199. With the numerically predicted 
fiber failure load considered as the failure load of the laminate, the difference be- 
tween the numerically predicted failure load and the average of the experimental 
failure loads is 3 %. 

Figure 7-21 shows uniformity of end loading for plates 3 and 4. With plate 3 the 
uniformity of loading was excellent before the mode shape change. After the mode 
shape change the side with LVDT 2 was displaced inplane slightly more than the side 
with LVDT 1. With plate 4, throughout the test the side with LVDT 1 is displaced 
slightly more than the side with LVDT 2. Compared with the quasi-isotropic plates, 
after the mode transition the uniformity of loading of plate 3 was about the same as 
plate 1, and the uniformity of loading of plate 4 was about the same as plate 2, after 
the mode transition. 


Load vs. out-of-plane deflection relations are shown in Fig. 7-22. It is noted that 
plate 4 deflected in the negative z direction. Correlations between experimental re- 
sults and the one half-wave configuration are good until mode transition occurs. At 
the load level of the mode transition the experimental relation showed some devi- 
ation from the numerical prediction. That is because the deformed shape just before 
the mode transition was not predominantly one half-wave. As the load level ap- 
proached the mode transition load, the deformed shape of the plate began to change 
from predominantly one half-wave to a mixture of one half-wave and two half-waves. 
As the load was increased, the deformed shape of the plate then changed rapidly to 
two half-waves. Numerically predicted failure loads for the two half-wave configura- 
tion are marked on the two half-wave relation. This relation corresponds to an iden- 
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Fig. 7-21 Comparison of displacements measured from LVDT 1 and LVDT 2 for plate 3 and 
plate 4. 
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Fig. 7-22 Comparison of load vs. deflection relations and failure loads for the [ ± 45/0 j] 2S 
plates. 
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tically zero deflection at the net-section hole edge. Again note the closeness of the 
failure loads for the two repeat tests. 


7.2.2 Load vs. Strain Relations 


Experimentally measured strains for plate 3 are shown in Fig. 7-23. Again, these 
grouping of strain gauges allow assessment of the overall as well as the detailed re- 
sponse of the plate. As with the quasi-isotropic plates, Fig. 7-4, strain gauge 1 and 
strain gauge 3 registered negative strains before and after the mode transition. 
However, the left-right symmetry of strains was not as good as it was with the 
quasi-isotropic plates. As was discussed for quasi-isotropic plates, by the fact that 
strain gauge 5 registered positive strain in the load range between the buckling load 
and the mode shape change load, it can be seen that before mode transition plate 3 
was deformed so that strain gauge 5 was on the convex side of the deflected plate. 

Experimentally measured strains for plate 4 are shown in Fig. 7-24. Strain gauge 
5 in plate 4 registered negative strain all through the test. From this it can be seen 
that plate 4 deformed so that strain gauge 5 was on the concave side. That is, as said 
before, the deflected direction of plate 4 was opposite plate 3. Strain gauge 1 and 
strain gauge 3 registered negative strains, as with the previous plates. However, it 
can be seen that the responses of strain gauge 1 and strain gauge 3 of plate 4 after 
the mode transition were different than the responses of those gauges for other 
plates after mode transition. With plate 4 after the mode shape change the magni- 
tudes of strains 1 and 3 were reduced, while in the other three plates the magnitudes 
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Fig. 7-23 Experimentally measured strains of plate 3. 
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Fig. 7-24 Experimentally measured strains of plate 4. 
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of these strains increased. Figure 7-25 explains the mechanism related to the above 
strain variation. Figure 7-25(a) shows the cases of plate 1, 2, and 3. In these cases 
the locations of strain gauges 1 and 3 were on the locally concave side of the plate. 
At the locations of strain gauges 1 and 3 the out-of-plane deflections were very small. 
Therefore, the strains were determined predominantly by inplane compression. 
However, the curvature due to the out-of-plane deflection did contribute to these 
strains. From Fig. 7-25(a) it can be seen that the curvature at strain gauges 1 and 3 
would increase, in magnitude, after the mode shape change. Therefore, strain 
gauges 1 and 3 became more compressive after the mode shape change. When the 
plate was deflected in opposite direction, as with plate 4 and as shown in Fig. 7-25(b), 
strain gauges 1 and 3 were on the locally convex side of the plate. After the mode 
transition the local curvature increased so as to produce tensile strain. Therefore, 
the negative strains at the locations of strain gauge 1 and 3 were relieved. In Fig. 
7-24, for plate 4 strains in strain gauges 4 and 6 were identical before the mode shape 
change. After the mode shape change these strains showed a big difference. 


Comparison of the experimentally measured and numerically predicted strains 
are shown in Fig. 7-26 - Fig. 7-28 for plate 3. In Fig. 7-26 the comparison of strains 2 
and 7 show good correlation. Before the mode shape change the measured strains 
agree well with the one half-wave prediction, and after the mode shape change the 
measured strains agree well with the two half-wave prediction. The comparison of 
strains at the net-section hole edge, strains 5 and 8, show excellent correlation in Fig. 
7-27. The comparisons for strain gauge 1, 3, 4, and 6 are shown in Fig. 7-28. It shows 
that the strains initially tended to follow the one half-wave prediction. After the mode 
shape change the strains generally followed the trends of the two half-wave predic- 
tion, though strain gauge 4 tended to be larger than the numerical prediction. 
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Fig. 7-25 Curvature change at the locations of sgl and sg3 after mode shape change. 
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Fig. 7-26 Comparison of strains near the clamped support of plate 3. 
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Fig. 7-27 Comparison of strains at the net-section hole edge of plate 3. 
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The comparison of strains for plate 4 is shown in Fig. 7-29 through Fig. 7-31. Re- 
call that plate 4 deflected in the direction opposite to plate 3. In Fig. 7-29 the com- 
parison of strains for strain gauge 2 and 7 show reasonable correlation. However, it 
can be seen in Fig. 7-30 that the strains at the net-section hole edge, strains 5 and 
8, show excellent correlation with numerical prediction. Relative to Fig. 7-28 for plate 
3, the load vs. strain relations for the one half-wave and two half-wave predictions for 
strain gauge 1 and strain gauge 3 are interchanged. Experimentally, the jump in the 
data for these gauges was opposite Fig. 7-28. Before the mode change the exper- 
imental data of gauges 1 and 3 followed the one half-wave prediction. After the mode 
change the data followed the two half-wave prediction. 


7.2.3 Shadow Moire Results 

The comparisons of shadow moire fringe pattern and numerically predicted con- 
tour lines of out-of-plane displacement are shown in Fig. 7-32 for the one half-wave 
configuration and in Fig. 7-33 for the two half-wave configuration. The load level for 
the one half-wave configuration in Fig. 7-32 indicated at the volt meter, 0.151 v, cor- 
responded to N„ = 824 lb/in., and for the two half-wave configuration, 0.315 v corre- 
sponds to N, = 1715 lb/in.. For the one half-wave configuration experimental 
out-of-plane deflection at this load was w/h = 3.0, and the numerical out-of-plane de- 
flection was w/h = 3.1. 
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Fig. 7*29 Comparison of strains near the clamped support of plate 4. 
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(a) strain 1 


(b) strain 3 




Fig. 7-31 Comparison of strains near the simple support of plate 4. 
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Fig. 7-32 Predicted and measured shadow moire fringe pattern for a [ ± 45/0 j], s plate with a 
one half-wave configuration. 
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Fig. 7-33 Predicted and measured shadow moire fringe pattern for a [ ± 45/0,]** plate with a 
two half-wave configuration. 
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7.2.4 Failure and Post-Mortem Analysis 


The photograph for the failure of plate 3 is shown in Fig. 7-34. It shows that the 
failure occurred at the node line of the two half-wave configuration. It can be seen the 
failure occurred at the simply supported edge, not the hole edge. The failure for plate 
4 also occurred at the simply supported edge. The contour plots of failure parame- 
ters at the load level of first fiber failure in the two half-wave configuration are shown 
in Fig. 7-35. From 7-35 (a) it can be seen that fiber failure is predicted at the net- 
section simple support. Figure 7-35(b) shows that the transverse failure region is 
wide spread in the center of the quarter-plate. From Fig. 7-35(c) it can be seen that 
inplane shear failure is predicted along the simple support and away from the net- 
section. From Fig. 7-35(d) and (e) it can be seen that the interlaminar shear failure 
parameters are very small over all of the plate. However, it can be seen that R 13 and 
Ri 3 are increasing rapidly at the net-section hole edge. 


7.2.5 Acoustic Emission Results 

The acoustic emission results are shown in Fig. 7-36 for plate 3 and in Fig. 7-37 
for plate 4. These results also show clearly the audible mode shape change and the 
final failure event. Figures 7-36 and 7-37 show some acoustic emission activity oc- 
curred before the mode shape change. The activity correlates well with numerically 
predicted transverse failure just before the mode shape change. These figures also 
show increasing acoustic emission activity near the load level at which the inplane 
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Fig. 7-34 Failed [ ± 45 / 02 ],! plate In the fixture. 
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(b) Contour plot of R t on 14-th layer from bottom ( 0° layer). 



<c) Contour plot of /?« on 16-th layer from bottom ( 45° layer). 

Fig. 7-35 Contour of failure parameters of the [ ± 45/0 plate at P/EA = 0.00206 (continued). 
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(d) Contour plot of Rn on 8-th layer from bottom ( 0° layer). 



(e) Contour plot of R n on 6-th layer from bottom ( -45° layer). 


Fig. 7-35 Contour of failure parameters of the [ ± 45/Oj]** plat® at P/EA- 0.00206 (concluded). 
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Fig. 7-36 Acoustic emission activity and load-end shortening relation for plate 3, [ ± 45/0j] js . 
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Fig. 7-37 Acoustic emission activity and load-end shortening relation for plate 4, [ ± 45/0j] js - 
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shear failure is predicted. Qualitatively, the the numerical prediction for this plate 
agrees with the acoustic emission test results. 


7.3 The [±4510^5 Plates 


7.3.1 Load vs. Deflection Relations 

Figure 7-38 shows experimental load vs. end shortening relations for [ ± 45/0 6 ] 2S 
plates, plate 5 and plate 6. With this laminate a mode shape change did not occur. 
Until the plate failed, the mode shape remained in the one half-wave configuration in 
both the loading and transverse directions. The figure shows the excellent repeat- 
ability of the load-end shortening relations for the plates. The numerically predicted 
buckling load shown in Table 2-3, P Cf /EA = 0.000245 , or Aft' = 318 lb/in., is marked 
on Fig. 7-38. Experimental buckling loads were P„/EA = 0.000241 , or N? = 301 lb/in. 
for plate 5, and P cr /EA = 0.000227 , or A/" — 283 lb/in. for plate 6. It can be seen that 
there was a difference in the experimental failure loads between plate 5 and plate 6. 
Plate 5 failed when P/EA = 0.00143 ( A/,= 1786 lb/in. ), and plate 6 failed when 
P/EA = 0.00162 ( N x = 2023 lb/in.). Figure 7-38 shows that the numerically predicted 
load-end shortening relation agrees very well with experimental results. The numer- 
ically predicted failure loads are marked on the numerical load-end shortening re- 
lation curve. The experimental average failure load of plate 5 and plate 6 was 
P/EA = 0.00154. The numerically predicted fiber failure load listed in Table 5-2 was 
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Flfl. 7-38 Comparison of load vs. end shortening relations and failure loads for the 
[ ± 45/0e]s plates. 
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P/EA = 0.00163. With the fiber failure load considered as the final failure load, the 
numerical prediction is in error relative to the experimental result by 6 %. 

Figure 7-39 shows the uniformity of loading for plate 5 and plate 6. It can be seen 
that the plate 6 was loaded with left-right uniformity to failure. However, with plate 5 
the uniformity of the loading was not maintained for P/EA> 0.001. After the tests it 
was found that plate 6 failed with about equal damage on both the left and right sides. 
However, plate 5 failed on only one side. From Fig. 7-39 it can be deduced that the 
failure of plate 6 progressed equally on both sides, and the failure of plate 5 began 
on the side of LVDT 2 and the damage on the side of LVDT 2 increased as the load 
increased. As that side of the plate was damaged more and more, the displacement 
increased more and more relative to the other side. The difference between the LVDT 
1 and LVDT 2 at the maximum load was 13 % of the average end displacement. 

The comparison of experimental and numerically predicted load vs. out-of-plane 
deflection relations is shown in Fig. 7-40. It should be noted that both plates deflected 
in the negative z direction. As with the load vs. end shortening relation of Fig. 7-38, 
plate 5 and plate 6 followed exactly the same load vs. out-of-plane deflection relation. 
In addition, the relation showed good agreement with numerical prediction. From Fig. 
7-40 it can be seen that plate 5 failed with a maximum deflection level of w/h = 4.55 
at the center of the plate, and plate 6 failed with a maximum deflection level of 
w/h = 5.0. This is a rather large out-of-plane deflection. The previous four plates had 
relatively small maximum deflection when they failed, since those plates were de- 
formed in the two half-wave configuration. 
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Fig. 7-39 Comparison of displacements measured from LVDT 1 and LVDT 2 for plate 5 and 
plate 6. 
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Fig, 7-40 Comparison of load vs. deflection relations and failure loads for the [ ± 45/0e]s 
plates. 
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7.3.2 Load vs. Strain Relations 


Figure 7-41 shows the experimentally measured strains for plate 5. It can be 
seen that the response of strain 1 and strain 3 agreed with the measurement of LVDT 
1 and LVDT 2 shown in Fig. 7-39, the displacement of LVDT 2 being greater than LVDT 
1. By the fact that strain gauge 5 registered negative strain until plate failed, it can 
be realized that the plate was bowed so that the strain gauge 5 was on the concave 
side of the deflected plate. The experimentally measured strains for plate 6 are 
shown in Fig. 7-42. From the response of strain gauges 1, 3, 4, and 6 it can again be 
seen that the response of plate 6 was more symmetrical than plate 5. Strain gauge 
5 of plate 6 also registered compressive strain until the plate failed. Therefore, it can 
be understood that plate 6 also deflected in the same direction as plate 5. 

The comparison of numerically predicted strains and experimentally measured 
strains are shown in Fig. 7-43 - Fig. 7-45 for plate 5 and in Fig. 7-46 - Fig. 7-48 for plate 
6. Generally, plate 6 shows better experimental-numerical correlation than plate 5. 
In light of the lack of left-right symmetry of the response with plate 5, this was ex- 
pected. 


7.3.3 Shadow Moire Results 

Figure 7-49 shows the moire fringe pattern at the load level of P/EA = 0.0008, or 
W x = 993 lb/in.. The pattern is compared with the numerically predicted out-of-plane 
displacement contour drawing. At this load level the experimental out-of-plane de- 
flection at the net-section hole edge was w/h = 3.17, and numerically predicted out- 
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Fig. 7*41 Experimentally measured strains of plate 5. 
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Fig. 7-42 Experimentally measured strains of plate 6. 
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Fig. 7-43 Comparison of strains near the clamped support of plate 5. 
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Fig. 7-45 Comparison of strains near the simple support of plate 5. 
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Fig. 7-46 Comparison of strains near the clamped support of plate 6. 
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Fig. 7-47 Comparison of strains at the net-section hole edge of plate 6. 
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Fig. 7-48 Comparison of strains near the simple support of plate 6. 
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Fig. 7-49 Predicted and measured shadow moire fringe pattern for a [ ± 45/Ojs plate at a 
lower load level. 
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Fig. 7-50 Predicted and measured shadow moire fringe pattern for a [ ± 45/0b] s plate at a 
higher load level. 
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of-plane deflection is w/h = 3.21. Figure 7-50 shows the moire fringe pattern at the 
load level of P/EA = 0.0013, or N x = 1636 lb/in., a higher load than the load in Fig. 7-49, 
and it's comparison with the numerically predicted out-of-plane displacement contour 
drawing. The experimental out-of-plane displacement at this load level was 
w/h = 4.20, and the numerically predicted out-of-plane displacement at this load level 
was w/h = 4.55. It can be seen that the fringe pattern at the higher load was quite 
different than the fringe pattern at the lower load. At the corners of simple support 
and clamped support displacements of opposite sign appeared. The numerical pre- 
diction also predicts this reverse displacement. Recalling the 3-dimensional surface 
drawing of the out-of-plane displacement at the high load level for quasi-isotropic 
plate in Fig. 5-4, this fringe pattern can be better understood. 


7.3.4 Failure and Post-Mortem Analysis 

A photograph of failed plate 5 in the fixture is shown in Fig. 7-51. It can be seen 
that failure occurred along the simple support on one side. Plate 6 also failed along 
the simple supports but on both the left and right sides. Compared with the numer- 
ically predicted failure locations in Fig. 5-12, the experimental failure locations and 
failure modes are completely different. In a [±45/0,] s laminate the inner 12 layers 
are all 0° layers. The fiber direction in the 0° layer is parallel to the simple support 
line. It can be easily understood that the inner 12 0° layers are weak in resisting the 
interlaminar shear stress t that is, t Z3 . The failure along the simple support could 
well be due to the interlaminar shear stress t m . However, the numerically 
predictedresults did not predict the mode of failure. The numerically predicted stress 
state in the laminate can be understood better by the contour drawing of failure pa- 
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Fig. 7-51 Failed [ ± 45/Ob]* plate in the fixture. 
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line of simple support 



line of clamped support 


(a) Contour plot of Ri on 2-nd layer from bottom ( —45° layer). 



(b) Contour plot of R t on 14-th layer from bottom ( 0° layer). 



(c) Contour plot of Rn on 1-st layer from bottom ( 45° layer). 

Fig. 7-52 Contour of failure parameters of the [ ± 45/0s] s plates at P/EA = 0.001 69 (continued). 
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(d) Contour plot of R, 3 on 9-th layer from bottom ( 0° layer) 



(e) Contour plot of Rn on 9-th layer from bottom ( 0° layer). 


Fig. 7-52 Contour of failure parameters of the [ ± 45/0*]* plates at P/EA =0.00169 (concluded). 
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rameters at the load level of fiber failure, as shown in Fig. 7-52. Figure 7-52(a) shows 
that ft, is high near the simply supported edge away from the net-section. Figure 
7-52(b) shows that ft 2 is high on the net-section toward the simple support. It can be 
seen that the maximum stress level is about twice the transverse failure load. Figure 
7-52(c) shows that ft« is high along the simple support and near the net-section hole 
edge. From Fig. 7-52(d) and (e) it can be seen that the interlaminar shear failure pa- 
rameters are small all over the plate. However, it is interesting that the interlaminar 
shear stress t 2J is relatively high along the simple support compared to other lo- 
cations. From this result it is surmised that since the inner 12 0° layers are clustered, 
the strength for of thick clustered layer was much lower than the strength for t 23 
of a single layer used in the numerical prediction. This point will be discussed more 
shortly. 


7.3.5 Acoustic Emission Results 


Acoustic emission test results are shown in Fig. 7-53 and Fig. 7-54 for plate 5 and 
plate 6, respectively. From the acoustic emission test results it can be seen that once 
the failure began, the acoustic emission count kept increasing at a constant rate. In 
previous discussions it was shown that acoustic emission count of the quasi-isotropic 
plate dramatically increased only at the final stage, and that the acoustic emission for 
[±45/0 2 ] 2S plate sharply increased shortly before the final failure. Considering the 
acoustic emission results of Fig. 7-53 in conjunction with the failure along the simple 
support, it can be imagined that the [ ± 45/0 6 ] s plate failed in a continuous rather than 
in a sudden fashion along the simple support. Even though the load vs. end shorten- 
ing relation did not show a large loss of load bearing capability, it can be imagined 
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Fig. 7-53 Acoustic emission activity and load-end shortening relation for plate 5, [ + 45/0 b] 
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Fig. 7*54 Acoustic emission activity and load-end shortening relation for plate 6, [ ± 45/0g] s . 
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that the plate was already damaged severely in transverse direction before final fail- 
ure, and it was much much weaker in the transverse direction than the original plate. 
Therefore it could be expected that the two previous laminates were less damaged 
than [ ± 45 / 0 6 ] 2 s laminate just before the final failure. 


7.4 The [±45] 4S Plates 


7.4.1 Load vs. Deflection Relations 


Experimental load vs. end shortening relations for the [ ± 45]« plates, plate 7 and 
plate 8, and their comparison with numerical results are shown in Fig. 7-55. Both 
plates exhibited a mode shape change, at about the same load, but the change was 
not abrupt. The numerically predicted buckling load, P cr /E A = 0.00098 , or A/?' = 240 
lb/in., is marked on Fig. 7-55. Experimental buckling loads were P er /EA = 0.00088 , 
or N? = 215 lb/in. for plates 7 and 8. Numerically predicted failure loads are marked 
on the numerical load-end shortening curves for the one half-wave and two half wave 
configurations. With the [ + 45]^ plate the predicted fiber failure load is much higher 
than the experimental failure load. Hence it does not appear in the load range viewed 
in Fig. 7-55. Up to the load level of mode shape change the two experimental re- 
lations for plate 7 and plate 8 followed exactly the same path. After the mode shape 
change the two relations showed a slight difference. Mode transition for plate 7 oc- 
curred at P/EA = 0.00300 ( N, = 734 lb/in. ) and for plate 8 at P/EA = 0.00294 ( 
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Fig. 7-55 Comparison of load vs. end shortening relations and failure loads for [ ± 45] 4 $ 
plates. 
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= 719 Ib/in. ). Plate 7 failed at P/EA = 0.00365 ( A/, = 892 Ib/in. ), and plate 8 failed 
at P/EA = 0.00350 ( N x = 856 Ib/in. ). Correlation with numerical result shows good 
agreement up to P/EA = 0.00189 ( N„ = 461 Ib/in. ) where the inplane shear failure is 
predicted. This load level was far lower than the mode transition load. It can be ob- 
served that as the plate began to fail, it lost load bearing capability. This resulted in 
an increased difference between the experimental results and the numerical predic- 
tion. At higher loads the numerical prediction does not agree with the experimental 
results at all. It may be worth mentioning that the mode shape change of this lami- 
nate occured very slowly, contrary to the previous laminates. It is deduced that due 
to the accumulated damage in the plates at the load level of mode shape change, the 
plate behaved differently than the previous laminates. 

Figure 7-56 illustrates the good uniformity of loading on the edge of plate 7 and 
plate 8, plate 8 being better as regards uniformity. Figure 7-57 shows the good 
agreement in the numerical and experimental load vs. out-of-plane deflection before 
the mode shape change. From Fig. 7-57 it can be seen that the plates 7 and 8 failed 
soon after the mode shape change into the two half-wave configuration. Also, the 
lack of a sharp mode transition is seen in this relation. It should be noted that plate 
8 deflected in negative z direction. 

7.4.2 Load vs. Strain Relations 


Figures 7-58 - 7-65 illustrate the strain response of plates 7 and 8. In Fig. 7-58 the 
data for strain gauge 2 of plate 7 is missing. In the testing of plate 7 strain gauge 2 
was not connected well. Thus the data could not be collected. By the fact that strain 
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Fig. 7-57 Comparison of load vs. deflection relations and failure loads for the [ ± 45]*$ plates. 
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Fig. 7-58 Experimentally measured strains of plate 7. 
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Fig. 7-59 Experimentally measured strains of plate 8. 
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Fig. 7-61 Comparison of strains at the net-section hole edge of plate 7. 
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Fig. 7-62 Comparison of strains near the simple support of plate 7. 
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gauge 5 registered positive strain in the deep postbuckling range before the mode 
shape change, it can be seen that the plate 7 was deflected so that strain gauge 5 
was on the convex side of the plate, as in Fig. 7-25(a). Strains 1 and 3, and 4 and 6 
showed good symmetry of response of plate 7 before the mode shape change. Fig- 
ure 7-59 shows the experimental strains for plate 8. By the fact that strain gauge 5 
registered negative strain, it can be seen that the plate 8 was deflected so that the 
strain gauge 5 was on the concave side of the plate. This was opposite the response 
of plate 7. Gauges 1 and 3 showed reasonable agreement, while gauges 4 and 6 
showed good agreement. 

The correlations between the experimentally measured strains and numerically 
predicted strains are shown in Fig. 7-60 - Fig. 7-62 for plate 7 and in Fig. 7-63 - Fig. 
7-65 for plate 8. The correlation between predictions and experimental results for the 
strains before the mode shape change is not bad. Before the mode shape change the 
experimental strains seemed to follow the one half-wave solution. However, after the 
mode shape change the experimental strains did not agree well with numerically 
predicted strains. It is noticable especially in the comparison of strain 7 in Fig. 7-60 
and strains 2 and 7 in Fig. 7-63. 


7.4.3 Shadow Moire Results 


The comparison of the out-of-plane deflection contour drawing and shadow moire 
fringe pattern for the one half-wave configuration is shown in Fig. 7-66. This corre- 
lation is good. For the two half-wave configuration the comparison of out-of-plane 
deflection is shown in Fig. 7-67. 
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= 0.00205 


Flfl. 7*66 Predicted and measured shadow moire fringe pattern for a [ ± 45]4s plate in a one 
half-wave configuration. 


Experimental Results and Comparison with Numerical Results 


2B3 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 




E^=0.00335 


Fig. 7-67 Predicted and measured shadow moire fringe pattern for a [ ± 45]** plate In a two 
half-wave configuration. 
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7.4.4 Failure and Post-Mortem Analysis 


The failed status of plate 7 is shown in Fig. 7-68. It can be seen that failure 
occured along a 45° line along one simple support. In plate 8 failure occured in the 
same fashion but along both simple supports. Even though the load-end shortening 
relation was quite different than the numerical prediction, the failure location seemed 
to agree with numerical prediction of Fig. 5-16. The contour drawing of failure pa- 
rameters at a load level slightly above the matrix failure load is shown in Fig. 7-69. 
It can be seen that R, is low overall, R 2 is high at the center of the quarter-plate, and 
Ri 2 is very high along the simple support slightly away from the net-section. Figure 
7-69(d) and (e) show the interlaminar shear stresses are very small all over the plate. 


7.4.5 Acoustic Emission Results 

Acoustic emission test results for plate 7 and plate 8 are shown in Fig. 7-70 and 7-71, 
respectively. Unfortunately, the contact between the specimen and acoustic emission 
transducer was bad. The results for plate 7 and plate 8 were poor and essentially 
meaningless. 


7.5 Summary of Findings and Further Results 


A large amount of information obtained from the experiments has been pre- 
sented which allows for the characterization of the postbuckling response of the 
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Fig. 7-68 Failed [ ± 45]*$ plate In the fixture. 
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(b) Contour plot of R t on 16-th layer from bottom ( 45° layer). 
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(c) Contour plot of on 1-st layer from bottom ( 45° layer). 


Fig. 7-69 Contour of failure parameters of the [ ± 45]« plates at P/EA = 0.00395 (continued). 
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(d) Contour plot of on 8-th layer from bottom ( 45° layer). 



(e) Contour plot of /?« on 8-th layer from bottom ( 45° layer). 


Fig. 7-69 Contour of failure parameters of the [ + 45]4s plates at P/EA = 0.00395 (concluded). 
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Fig. 7-70 Acoustic emission activity and load-end shortening relation for plate 7, [ ± 45] 4S . 
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Fig. 7-71 Acoustic emission activity and load-end shortening relation for plate 8, [ + 45] 4$ . 
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specific composite plates studied. In general, the numerically predicted response 
and the observed response agree quite well. Though the mode shape change was 
not predicted, nor was there an attempt to predict it, it was shown that before the 
mode shape change the response of the plates followed the one half-wave config- 
uration prediction, and after the mode shape change the response of the plates fol- 
lowed the two half-wave configuration prediction. The numerical prediction of failure 
loads for f + 45/0/901,.= and £ + 45/0,1,.= laminates agrees very well with the exper- 
imental results. It was shown that the failure of these laminates was controlled by 
inplane stresses. 

With the [±45]« laminate the postbuckling response did not agree with the nu- 
merical prediction. It was observed that the deviation of the numerical prediction from 
the experimental results began at the load level where the inplane shear failure was 
predicted. It was mentioned in chapter 7 that due to the inplane shear failure the 
[±45]« laminate loses load bearing capability. This resulted in the disagreement 
between the numerical prediction and experimental results in postbuckling response. 
Therefore, predicting the failure load with the present method is not accurate. It is 
recommended that progressive failure analysis be used for the [±45]^ laminate. 
Other researchers have had difficulty in correlating the observed response of ±45 
laminates with numerical prediction based on linear elastic analysis. The nonlinear 
softening stress-strain behavior associated with the inplane shear response is gen- 
erally felt to be the cause of this. Progressive failure would add a softening effect 
also. It is felt both issues must be addressed to fully understand the effect of lami- 
nates with no 0° layers. However, it was shown that the prediction of failure mode and 
failure location agreed with experimental results. 
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With the [±45/0 6 ] s laminate predicted postbuckling deformation response 
agreed very well with numerical results. However, predicted failure location and fail- 
ure mode did not agree with experimental results. In the numerical prediction no 
interlaminar shear failure was predicted. However, it appears that the [ + 45/0 6 ] s 
plates failed due to excessive interlaminar shear stress r« along the simple supports. 
It is felt that the interlaminar shear strength of the clustered layers, the inner 12 lay- 
ers of the [ ± 45/0e]s laminate, could have been much lower than the interlaminar 
shear strength of a single layer. To explore this idea further, based on the exper- 
imental observation of the [ ± 45/0 s ] s laminate, the interlaminar shear strength r 23 for 
the clustered inner 12 layers will be estimated. The damage accumulation scenario 
will then be re-examined for this laminate. 

The average experimental failure load of the plates 5 and 6 is N„ = 1918 lb/in.. In 
the numerical failure analysis the maximum value of failure parameter R t3 was 0.25 
when N, = 1918 lb/in. if the interlaminar shear strength S 23 = 8 ksi was used, the 
value shown in Table 5-1 . Therefore, it can be deduced that if S 23 of the clustered inner 
12 layers of the [ ± 45/0 6 ] s laminate was 2 ksi, the interlaminar shear failure parame- 
ter would be R 33 = 1.0. Figure 7-72 shows the re-examined damage accumulation 
scenario of the [ ± 45/0$] s laminate deformed into the one half-wave configuration 
with estimated interlaminar shear strength of S 23 = 2 ksi. The other strength values 
are the same as in Table 5-1. In Fig. 7-72 it can be seen that interlaminar shear failure 
due to t 23 is predicted along the simple support. Away from the net-section but along 
the simple support over one-half of the inner 12 layers have failed. Figure 7-73 shows 
the damage accumulation scenario for the [±45/0 e ] s plate with a hole deformed into 
the two half-wave configuration. It can be seen that severe interlaminar shear failure 
due to T 23 is predicted for the two half-wave configuration also. 
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Fig. 7-72 Predicted failure locations and failure modes of [ ± 45/0b]s plate with hole In the one 
half-wave configuration with revised value of S». 


Experimental Results and Comparison with Numerical Results 


• Fiber failure 
— Trans, failure 
▲ Shear failure 



(Number in the figure represents the no. of failed layers.) 


Load: -f_ = 0 .00161 
EA 


Fig. 7-73 Predicted failure locations and failure modes of [ ± 45/0e] s plate with hole in the two 
half-wave configuration with revised value of S a . 
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Based on the comparison between the experimental results and numerical 
predictions for the plates with holes, the following can be said for the plates without 
holes. It is expected that the postbuckling deformation response and failure of the 
solid [±45/0/90] 2S and [ ± 45/0 2 ] 2S laminates can be correctly predicted as in chap- 
ters 3 and 5. The postbuckling deformation response of the solid [ ± 45/0 $ ] s laminate 
can be predicted correctly as in chapter 3. However, for this particular laminate fail- 
ure analysis using the lower interlaminar shear strength value, S 23 = 2 ksi, will prob- 
ably give a more realistic prediction. The reevaluated failure prediction for the one 
half-wave configuration is shown in Fig. 7-74. It can be seen that the solid [ ± 45/0 6 ] s 
plate in the one half-wave configuration would experience failure due to interlaminar 
shear stress r 23 along the simple support. In the same region a high fiber stress is 
also predicted, along with a high inplane shear stress. Figure 7-75 shows the failure 
prediction of the solid [ ± 45/0 6 ] 3 plate deformed into the two half-wave configuration. 
In Fig. 7-75 it can be seen that interlaminar shear failure due to x 23 is not as severe. 
The solid [ ± 45/0 e ] s laminate in the two half-wave configuration will fail due to in- 
plane stresses. 
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Fig. 7-74 Predicted faliure locations and failure modes of solid [ ± 45/0b] s plate in the one 
half-wave configuration with revised value of S w . 
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Fig. 7-75 Predicted failure locations and failure modes of solid [ ± 45/0b] s plate In the two 
half-wave configuration with revised value of S^. 
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8.0 Concluding Comments and Recommendations 


8.1 Concluding Comments 


The numerical buckling loads obtained in chapter 2 and the experimental 
buckling loads discussed in chapter 7 are summarized in Table 8-1. Numerical and 
experimental failure loads of each of the laminates discussed in chapters 5 and 7 are 
summarized in Table 8-2. Table 8-2 shows a failure prediction resulting from the use 
of a revised value of S M for the [ ± 45/0$] s laminate. Since known failure results were 
used to obtain the revised prediction of Table 8-2, the "prediction" is actually mean- 
ingless. What is important is the fact that the location of failure associated with the 
"prediction" is accurately predicted to be along the simple support. 

From the results in Table 8-2 it can be concluded that the understanding of the 
failure in the postbuckling range has certainly been enhanced by this study. Failure 
of two laminates of practical importance, namely the [±45/0/90]^ and the 
[±45/0 2 ] m laminates, appears to be completely predictable. Highly orthotropic lami- 
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Table 8-1 Experimentally determined and numerically predicted buckling loads of each 
laminate. 
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plate 7 0.000880 215 

[±45] 4t 0.000880 215 0.000980 240 

plate 8 0.000880 215 

























Table 8-2 Experimentally determined and numerically predicted failure loads of each laminate. 
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. interlaminar shear failure 2 : fiber failure, one half-wave (Table 5-2) 
: Using the revised interlaminar shear strength to take into account the 
layer clustering effect along the simple support 



nates such as [± 45/0 8 ]s are not likely to be used in the near future. In fact, the 
[ + 45/0 2 ] 2 s laminate is quite orthotropic by current standards. Thus, not being able to 
a priori predict the failure of the [ + 45/0 8 ]s laminate may be of little practical conse- 
quence. Also, not using any 0° layers is uncommon and not being able to predict the 
failure, or even the deformation response, of the [ ± 45]« laminate may also be of 
little practical importance. In practice the issue of interlaminar shear stresses on the 
control of failure may not really be an important issue, either. For an aspect ratio of 
1:1, at least for the value of D/W considered here, whether the [ + 45/0/90] 2S and 
[ + 45/0 2 ] 2S laminates are deformed into the one half-wave configuration or the two 
half-wave configuration, the failure load and failure mechanisms are very similar. 
Specifically, failure due to excessive inplane stresses is the controlling factor. 


Based on this study of the postbuckling failure analysis, the following specific 

conclusion can be reached: 

1. Despite the presence of multiple waves in the loading direction, inplane failure 
mechanisms control the failure of square [ ± 45/0/90] K and [ ± 45/0 2 ] 2S laminates 
with holes of D/W = 0.3. 

2. In the [± 45/0 6 ] s laminate interlaminar shear failure governs failure. From the 
experimental results it was observed that the lower interlaminar shear strength 
of clustered layers could be responsible for controlling interlaminar failure along 
the simple support. 

3. The [±45]« laminate response is influenced by softening behavior, a character- 
istic not modelled. This characteristic must be modelled to accurately predict 
both the displacement response and failure. 
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The following general conclusions can also be reached as a result of this study: 


1. Even though a plate may experience a mode shape change in the postbuckling 
range, the failure load of a plate with a mode shape change may not be too dif- 
ferent from the failure load of a plate without mode shape change. This makes 
the designer's task much easier. 

2. The maximum stress criterion is valid for the postbuckling failure analysis. 

3. Since a criterion for the mode shape change in the postbuckling range is not 
available, experimental verification of numerical prediction is important. 

4. The use of acoustic emission as an indicator of failure for structural level exper- 
iments has potential. 


8.2 Recommendations for Future Research 


Clearly the primary recommendation to be made regarding future research is to 
repeat this research for plates with other values of D/W and for plates with no holes. 
Since the [ ± 45/0»] s and [ ± 45]4s laminates may not be of much practical importance, 
it is recommended that besides the [ ± 45/0/90]* and [ ± 45/0?]* laminates only one 
other laminate be considered. The three laminates to be considered are 
[ ± 45/0/90]*, [±45/0?]*, and [±45/90?]*. These laminates would bracket the 
range of orthotropy in current use, though the [ ± 45/90?]*, like the [ ± 45]*, does not 
contain any layers with fibers in the loading direction and therefore may, in fact, not 
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be realistic. The issue of plates with no holes is quite important. These could be 
considered the baseline cases and experimental data need to be gathered. Plates 
with larger holes than those tested would result in even more load being transmitted 
along the simple support region at the net-section. It would be expected that inplane 
stresses would control failure in those cases. For plates with smaller holes it is not 
apparent what would happen. If the hole is small enough, the plate may behave as if 
it had no hole. These are important issues. All aspects should be explored. Also of 
importance is the issue of mode change. Factors which control the triggering of mode 
change should be explored. Candidate factors include lack of perfect flatness in the 
plates, slight misalignments in the load, and a lack of ideal boundary conditions. 
Nonuniformity of material properties, in a spatial sense, could be important too. 

Secondary issues recommended for future research includes consideration of 
other materials, e.g., IM7/8551-7a, a material with stiffer fibers and a softer matrix. 
Absolute size may be another factor. Repeating the experiments with plates twice or 
three times as large may reveal that material allowables need to be adjusted to re- 
flect a large volume of material. 

In summary, more work can and should be done, both experimentally and nu- 
merically. The ultimate goal should be contributing to guidelines for designers to take 
advantage of postbuckling strength. 
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